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On Whisker Growth Shapes} 


By 8S. AMELINCKX 
Laboratorium voor Kristalkunde, Rozier, 6, Gent, Belgium 


[Received January 17, 1958] 


BRENNER (1956) has observed helically twisted silver whiskers grown by 
thermal decomposition of the silver halide. From his photograph it is 
quite clear that the side faces of the whisker are formed by parallel facets, 
indicating that the whisker is most probably a single crystal. 

Frank (1956) suggested that this particular shape results from the 
circular movement of a catalytic particle at the tip of the whisker. The 
present author suggested at that time that a helical dislocation of which the 
end point describes a circular path during growth, could provide an 
adequate explanation as well. Since then, vapour grown whiskers of 
PbTe, presenting the same striking shape, have been produced at the 
General Electric Laboratories (Apker 1957, private communication). It is 
doubtful that in the latter case catalytic action is responsible. We feel 
therefore that it is not unreasonable to assume that a helical dislocation 
is at the basis of this particular growth shape. In this note we describe in 
some detail a possible mechanism which would yield kinked and helically 
wound whiskers. 

In view of the difficulty of representing the mechanism in space we 
illustrate it by means of two dimensional analogon, which can however 
immediately be transposed to three dimensions. 

Suppose that we have a perfectly straight whisker, growing at the tip, 
by means of the screw dislocation mechanism (Frank 1949, Sears 1955). 
The screw dislocation is stable at the centre of the whisker, as shown 
by Eshelby (1953). There is even an attraction of the screw towards the 
axis of the cylindrical crystal. The screw will leave the cylinder only if it 
is brought by some externally applied stress nearer than about half way to 
the surface of the cylinder. 

We further assume that the habit ‘faces’ of the “two dimensional 
whisker ’ are (01) and (11). When we start with a whisker having the 
shape shown in fig. 1 the screw will have no tendency to be displaced, 
unless by some climb due to condensing vacancies, as the crystals cool. 
The dislocation may then acquire helical shape (Amelinckx et. al., 1957) 
this however cannot be represented adequately in a plane. 

As long as the end facet is present, the displacement of the emergence 
point will lengthen the dislocation. As the facet containing the emergence 
point of the screw will however grow faster than the others, the end facet 
will practically disappear at some moment, we then have the situation 
shown in fig. 2. The top facet has now the size of the critical nucleus, 
which is very small for the supersaturation of interest here. 


+ Communicated by Professor Dr. W. Dekeyser. 
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We now can distinguish between two possibilities : the screw maintains. 
its direction or is able to move slightly. In the first case, it will stabilize 
the shape of the whisker tip. This will be such that the intersection 
point of all end facets coincides with the mean position of the emergence 
point of the screw. All end facets will grow alternatively by growth layers 
initiated at the screw as it moves on its circular path. The probability 
that this happens is largest at lower temperatures. 


Figs. 1-7 


Ne) ocean 


i 
(7) 


’ 5 
Successive stages of the developments of a (two dimensional) helical whisker 
The displacement of the dislocation is indicated by the arrows 


If the screw is however able to move at higher temperatures it will 
do so because it implies a decrease in length. This leads to the configura- 
tion shown in fig. 3 where the screw now ends normally to one facet. 
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The actual shape is determined by a balance between the line tension 
which tends to shorten the dislocation, the surface energy which tends to 
shorten the exposed edge, and the attraction of the screw towards the 
axis of the cylinder. 

The facet AB will now grow more rapidly than the other as it contains 
the screw ; a situation as fig. 4 results. At that point the screw can 
again shorten by moving over towards the centre of the facet BC. This 
facet can now serve as a base for the development of a little section of 
whisker, where again the same phenomenon will repeat. The sequence 
of events is represented schematically in figs. 1, 2, 3,4,5,6and7. Inthe 
two dimensional model the dislocation ‘ wiggles’ from left to right and 
adopts a wavy shape, along the ‘axis’ of the zig-zag shaped whisker. 
The movement of the screw is in great part conservative. Ifthe plane of 
the drawing is a glide plane it is entirely a glide motion. In the three 
dimensional case however a small amount of climb will be necessary. 

The driving force for the climb motion is in the first place the line 
tension and in the second place the small supersaturation of vacancies which 
one may expect immediately behind the growing tip. At the temperature 
of growth, the dislocation movement should not be too difficult. The 
vibration of the whiskers during growth may also favour displacement of 
the dislocation. 

The fact that climb is involved induces the screw dislocations to adopt 
a helical shape, so that even in the absence of the line tension, the 
emergence point would tend to describe some circular path as the crystal 
lengthens. 

The dislocation which is left in the coiled crystal, has helical shape and 
is of mixed character. Cottrell (1956) has demonstrated that such a 
situation is a stable one. 

It is clear that not all whiskers will present this particular shape. It 
is quite well possible that the situation shown in fig. 2 subsists through 
the whole period of growth of the whisker especially at the lower 
temperature. In that case the whisker remains perfectly straight, with a 
tip exhibiting end facets. It may also happen that as a consequence of 
some growth accident the screw leaves partly or completely the whisker. 
In the first case the screw will end in a side face, and at high supersatura- 
tion a side branch may now develop on the whisker; as Brenner has 
sometimes observed. At lower supersaturation the whisker may become 
a blade. A screw may also stabilize after one kink, e.g. in the situation 
shown in fig. 5. One will then end up with a kinked whisker. In 
fact the mechanism may stop accidently in any one of its phases, so that 
partially coiled and partially straight whiskers, may well occur. It is 
clear that this mechanism only applies to whiskers growing at the tip. 

The square ‘ tubes ’ of platinum described by Brenner (1956) probably 
originate in a normal whisker, where growth was inhibited in the centre 
of the terminal facet, e.g. by a foreign particle. A mechanism quite 
similar to the one described previously for the growth of ‘ hopper ’ 
crystals (Amelinckx 1953) has probably operated in this case. 
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ABSTRACT 


The equivalent orbital method of describing the electronic structure of 
valence crystals, introduced in an earlier discussion of diamond, is extended 
in various ways. The interactions of more distant bonds are now included 
so that a more detailed description of the electronic structure can be given. 
The role of the equations as interpolation formulae is illustrated and used to 
deduce the effective masses of holes for diamond, silicon and germanium. 
The application of the method to the conduction bands is also considered. 


§ 1. INTRODUCTION 


THE self-consistent field treatment of the electronic structure of a crystal 
is based on the use of a determinant as an approximate wave-function 
for the crystalasa whole. The elements in this determinant are the orbitals, 
functions of the position of one electron, each multiplied by a spin function. 
In order that a wave-function of this form should give the best energy, 
these orbitals must satisfy certain integro-differential equations. These 
equations alone do not define the orbitals uniquely; any unitary trans- 
formation of the orbitals is another solution and gives a determinantal 
wave-function of exactly the same form and value. This freedom is used 
in the equivalent orbital theory to define two different types of orbital. 
The crystal orbitals are defined by the additional condition that the 
equations should reduce to eigenvalue form and so make the orbitals 
independent of one another. The equivalent orbitals are defined by the 
additional condition that they should be invariant or permute into one 
another under the operations of the symmetry group of the crystal. Thus, 
in diamond, each equivalent orbital has a centre of inversion which leaves 
it invariant and every other operation permutes it into another equivalent 
orbital. It is not obvious that it is always possible to find such orbitals 
and the group theoretical conditions required have been discussed 
previously (Hall 1950). The remarkable fact is that these conditions are 
always satisfied in the equilibrium nuclear configuration of any molecule 
for which a classical valence diagram can be drawn, and that the symmetry 
of the bonds is the same as the symmetry of the equivalent orbitals. For 
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this reason we may think of the equivalent orbitals as giving a quantum 
mechanical realization of the chemical bonds. 

In diamond the transformation to equivalent orbitals is only possible 
within a limited range of CC internuclear distances. Thus, as the calcula- 
tions of Kimball (1935) show, there is a crossing of some of the valence 
and conduction bands when the CC distance is about twice its equilibrium 
value and another when it is about half this value. Outside this range 
the conditions for a transformation to equivalent orbitals of this type are 
not satisfied and there are no chemical bonds. It is probable too that 
outside these limits the use of a wave-function consisting of a single 
determinant no longer corresponds to the ground state. This phenomenon 
of a definite internuclear distance at which bonds are broken, is also 
found in many molecules. 


§ 2. THe ENERGY ZONES 


The equivalent orbitals, from their definition, can be transformed 
exactly into the crystal orbitals by a unitary matrix. This gives a new 
method of discussing the band structure of crystals. The energies of the 
crystal orbitals are the eigenvalues of the matrix whose elements are the 
energies of interaction of the equivalent orbitals. The remarkable 
constancy of bond properties between small molecules and crystals 
suggests that in practice the equivalent orbitals are not widely spread, 
so that all but a finite number of these interactions between a bond and 
its neighbours will be negligible. In these circumstances the energies 
of the crystal orbitals can be deduced in terms of the remaining equivalent— 
orbital interaction energies. In principle all of these interactions could 
be included and knowledge of them is exactly equivalent to knowledge 
of the crystal orbital energies. 

In an earlier paper (Hall 1952) this method was used to discuss the 
valence bands of diamond. Only two parameters were used a, the ioniza- 
tion energy of an equivalent orbital, and 6, the interaction energy between 
two neighbouring equivalent orbitals. It was found that two of the four 
bands were broad, with widths depending on b, but two remained degenerate 
with energies independent of wave-vector. This was clearly an over- 
simplification but there was no information available then to enable further 
parameters to be estimated. Much more information about diamond-type 
crystals is now available and a more detailed treatment can be given. 

Each bond in a diamond lattice has six nearest neighbours. These 
were discussed previously. There are also two kinds of next nearest 
neighbour. Six of these are parallel to the original bond and may be 
called trans bonds and twelve are non-parallel cis bonds. The corres- 
ponding interaction energies will be denoted by c and d respectively. 
The matrix of interactions e,, between the equivalent orbitals y, can now 
be set up. The eigenvalues of this matrix, which are the energies of the 
crystal orbitals, are then found by solving the determinantal equation 


|e,,;—#8,, | =0. Se ae Be nll) 
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Alternatively, the matrix 7',,, which reduces ¢,, to diagonal form, may be 
found by solving the system of homogeneous linear equations 


Deieritaly=O, fel,2,... . . . « (2,2) 
8 


This matrix is largely determined by the translational symmetry of the 
crystal. 

The unit cell of the diamond lattice is a parallelepiped spanned by three 
basic vectors a,, a,, a3, and containing two atoms. The two atoms in the 
(J, m, n) cell, where J, m, n are integers, are located at 


a=la,+ma,+ Naz 
a’ = (I+ f)a + (m+ t)as+ (n+ t)as 


The unit cell contains four bonds. That joining the two atoms has an 
equivalent orbital which will be labelled y (J, m,n, 0) and that joining 
a’ to a+a, as y(/,m,n,r). The translational symmetry suggests the 
substitutions 


Tegel je il s={Ln, 7,0) 

On ry lies =U, mM, 0,7), Ve ee iC aeS) 
where K labels the crystal orbitals and is a vector in the reciprocal lattice, 
viz. 

==) Di d-0Dn a UDa 2, 8 sake eo (24) 


and 
b,=6,;;. . . . . . . . : (2.5) 


The effect of these substitutions is to reduce the infinite system of eqns. 
(2.2) to four equations 


ee f= 0, 601, 2, 3. ren) 
§ 
In terms of the parameters a, b, c, d the symmetrical matrix «,, has elements 
Xo =a + 2c{cosf+cosg+cosh}, 
%4,4= 4+ 2c{cos f+ cos (g—f) + cos (h—f)}, 


doo = a+ 2c{cosg + cos (h—g) + cos(f—g)}; 

%33= a+ 2c{cosh + cos(f—h) +cos (g—h)}, 

Xo1 = Hy) = 2b cosh f+ 2d{cos (gy — $f) + cos (h—$f)}, (227) 
X99 = Xo = 2b cos $9 + 2d{cos (h — $9) + cos (f—39)}, 

Xog = Aq = 20 cos $h + 2d{cos (f— 4h) + cos (g— 3h) }, 


O45 =U, = 2b cos (g —f) + 2d {cos $(g +f) + cos (h— 
013 = %g = 2b cos $(h —f) + 2d {cos $(h +f) + cos (9 — 
hog = Ugo = 2b cos $(h —g) + 2d{cos £(h +g) + cos(f— 
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The elimination of the P, , from eqn. (2.6) leads to a determinantal equa- 
tion for the crystal orbital energies L,, viz. 


|&, = 2 0,, |e Ore vn es 


This equation has four roots and these give the values of Hy in the four 
valence bands. The corresponding crystal orbitals will be 


= > {Py eo *y(l, m,n, 0) +P, on Bey m,n, 1) 
Imn 


+ Py ,eik@+ia)y(1, m,n, 2) +P ge -@+i%)y(L, m,n, By} 5 ne ee 


In the previous treatment it was always possible to factorize the 
eqn. (2.8) but this is no longer possible for arbitrary values of K. When 
K is restricted to certain lines, fixed by the symmetry of the energy surfaces, 
factorization is still possible. Thus, along the line h=0, g=f the four 
energy bands are 


Ey 7, =a+ 264+ 2c+ 4(c+d)cosf+4(b+2d) cos$f, . (2.10) 

a—2b+2c+4(c—d)cos f (twice). ion ae A (2 

Similarly along the line h=g=/f, or the line g=h=0, the four bands are 
Ey, 7, p= at 264+ 2¢+ 2d + (4c + 2d) cos f+ 24/{[b+ce+d+ (d—c) cosf}? 

4+-3(0-4:2d)* costa fe. 2 50. Gn See oe A en ee 

a—2b+ 4c — 2d + (2c — 2d) cos f (twice). — (2.13) 


These formulae show that the two bands, which were degenerate before, 
now have a width proportional to (c—d). At (0,0,0) there are two 
energies, the lower at 

Eo,o,jp=%+6b+6c+12d . . . . . (2.14) 


is singly degenerate, while the upper at 


Hoo p= a= 204604 ee 
is triply degenerate. 


§ 3. DETERMINATION OF PARAMETERS 


This description of the valence bands of a diamond-type crystal now 
contains four parameters which could be determined in various ways. 
In principle it might be possible to solve the equation defining one equi- 
valent orbital since this one function defines all the matrix elements, but 
in practice this is too difficult. The previous paper used values of a and 
6 which were determined purely empirically from the ionization potentials 
of paraffin molecules. Unfortunately these potentials are not known 
with sufficient accuracy or in sufficient numbers to enable further para- 
meters to be determined. 


q 
; 
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There are many advantages in a purely empirical determination of the 
parameters. The greatest is that the labour of calculation is avoided. 
There are risks of overfitting or underfitting the experimental data but, 
_ if the experimental accuracy is known, these risks can be estimated by a 
x” significance test. One method of determining the parameters 
empirically might be to use the x-ray emission spectra of the crystals. 
The width of the spectrum, for example, is a measure of 
8b + 16d ek ae, a eG 
but, unfortunately, the width is not well defined experimentally because 
of the background. The development of this method will depend on a 
better understanding of the details of these spectra. 

Another method of determining the parameters is by using the formulae 
to fit the results obtained by other methods of calculation. This method 
was applied in the previous paper to the calculations of Kimball (1935) on 
diamond. The method has since been applied to metals (Hall 1953, 
Slater and Koster 1954). It is particularly useful as an interpolation 
formula to deduce the energy of any crystal orbital from the small number 
of energies which can be calculated directly by other methods. The 
calculations of Kimball are not accurate enough to determine c and d 
but those of Herman on diamond (1952) and germanium (1954) by the 
orthogonalized plane wave method, and of Jenkins on silicon (1956) by the 
variational cellular method, have been fitted in this way and the resulting 
values of the parameters are shown in table 1. 


Table 1. Parameters Determined by Fitting other 
Calculations (in Rydbergs) 


Diamond Silicon Germanium 


—1-547 —0-353 —Q-248 
—(Q-135 —0-058 —0-110 

0-044 0-011 0-020 
—0-032 —0-031 0-004 


Despite the labour of their calculations, these authors do not claim 
high accuracy in their results, primarily because of the uncertainty of the 
crystal potential. For this reason these values of the parameters may not 
be very accurate despite the closeness of the fitting. 


§ 4. SHAPE AT THE TOP OF THE VALENCE BANDS 
In diamond, silicon and germanium the top of the valence bands occurs 
when K = 0 at the energy given by (2.15). The form of the energy surfaces 
in this neighbourhood are of considerable importance in understanding 
many physical properties of the crystals. 
Three energy surfaces meet at the top of the valence band and in. 
consequence they vary in a way too complicated to be described by an 


434. G. G. Hall on the Hlectronic 


effective mass. To determine their form, the determinantal eqn. (2.8) 
is first transformed using the matrix 


1 1 1 1 


1 1 -—-1 -Il 
1/4 4.1 
Peel 1 —-l ba 
1 -l1 -l 1 


This is the matrix of eigenvectors appropriate to K=0. The cosines in 
the elements of «,, are then expanded up to the quadratic terms. Because 
of the transformation using (4.1) it is found that, to the same accuracy, 
one eigenvalue has separated from the others. This is the eigenvalue 
corresponding to the lower singly degenerate energy surface. The 
remaining equation is cubic and describes the upper three surfaces in the 
neighbourhood of K=0. 

It is convenient now to introduce rectangular coordinates in order to 
clarify the relation between K and the momentum. Unit vectors i, j, k 
are taken along the sides of the face-centred cubic lattice formed by the 
aatoms. If the distance between two neighbouring atoms is p, the basic 
lattice vectors can be expressed as 


ay qi + j), ] 
a,=q(it+k), } (4.2) 
a;=qj+k), | 
where 
qv/3=2p. (4.3) 
This implies that, if 
K=k,i+k,j+k,k, a bee ela 


its components are related by the equations 
f= dk, +ky), 
g=qk,+k,), er eee ee 
h=qk,+k,). 


The cubic equation for the energy surfaces near K=0 in terms of these 
new variables is 


Ak? +(B—A)k,?—e Dk,k, Dk, ke, 
Dkk, Ak? + (B—A)k,?—« Dk,k, “6 
Dk,k, Dk,k, Ak? +(B—A)k2—e 
where ~ aa (AOR 
b= he + hy? + ke? (4.7) 
and 
A= 29°(d—c), (4.8) 
B= D=¢(4b—2c—d), (4.9) 


——— =. >. — — 
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An equation of the form (4.6) has also been deduced for these surfaces by 
Herman (1954) using perturbation theory and with a different inter- 
pretation of A, B, and D. This equation requires to be modified when the 
original Hamiltonian is modified by including spin-orbit coupling terms. 
The effect of these terms has been considered by Elliott (1954) and he shows 
that at K=0 the triplet is split into a doublet, which is uppermost, and a 
singlet. The shape of the surfaces in the neighbourhood of the doublet 
is given by eqn. (21) in his paper which should be corrected to read 
e= 3(24 + B)k* + 3{(B—A)*k* + 3[D?— (B—A)?](k,2h,? + 
+ hth? + k2h.*) 2 Sy Dh ee erat oe en ee aa 
and, near the singlet, 
e=4(2A + B)k?. any eee ae) 
Thus the singlet energy surfaces are spherical and can be described 
adequately by a single effective mass but the doublet energy surfaces 
vary with angle in a more complicated fashion than can be described 
adequately by effective masses. It has been found experimentally that 
cyclotron resonance measurements can be explained by using a formula 
of this type, viz. 
e= —f?/2m[ak? + {Brkt + 7(k 2k 2+ kh 2+ k2k,7) $2], . (4.18) 
and experimental values of «, 8, y have been found for silicon and 
germanium. 
A theoretical estimate of «, 8, y can now be given using the eqns. (4.8), 
(4.9) and identifying terms in (4.11) and (4.13). When a, 0, c, d are given 
in Rydbergs and q in atomic units, the resulting equations for «, 8, y are 


a= —@q?(b/6—2c+d), | 
B= —@(b/6—d), a ee Te) 
y? = 3q*(d —c)(b— 2c — 4d). J 


These equations show that «, 8, y depend sensitively on the values of ¢ 
and d rather than on 0 and, unfortunately, they depend even more on the 
values of other smaller parameters which cannot yet be estimated. This 
sensitivity is clearly to be expected, since the higher order parameters 
describe the finer details of the surfaces, and the behaviour at a point is a 
fine detail. 

The results of using (4.14) in conjunction with the values of 6, c, d given 
in table 1 are shown in table 2 together with the experimental values 
(Lax, as quoted by Herman 1955). The calculation for diamond has 
been done in the same way as for silicon and germanium although the 
spin-orbit terms are much smaller and the separation of the doublet 
and singlet is very small. The calculations imply that the holes in the 
three bands of diamond all have effective masses of about 0-5 to 0:8 electron 
masses. The poor numerical agreement in table 2, particularly for 
germanium, is probably due to the inadequate accuracy of the calcula- 
tions from. which the parameters were determined. The eqn. (4.9) also 


436 G. G. Hall on the Electronic 


contains the relation, B=D, and because of this only two of the experi- 
mental quantities are independent. The relation between them is 
Sy2=(a—B)\(at5B)  .) 2°. 2 2 (ae 
and this relation is satisfied by the experimental values to within the 
experimental errors. 
§ 5. THE ConpuctTion Banps 
Slater (1956) has suggested that the conduction bands in these valence 
crystals should be considered using antibonding equivalent orbitals each of 
which has a node between the two atoms with which it is associated. 
Although this is probably possible, its justification is more difficult than 
that for the use of the bonding equivalent orbitals to describe the valence 


bands. 
Table 2. Calculated and Experimental Values of «, B, y 


Diamond cale. Si cale. Si expt. Ge cale. | Ge expt. 
o 1-61 1-63 4-] 1-54 13-0 
|p | O-11 0-56 1-4 0-63 8-7 
ly | 0-78 ()-92 By | | 1-19 11-4 


The conduction bands of a crystal may most easily be defined as the 
unoccupied bands obtained as higher solutions of the self-consistent 
equations for the crystal orbitals. The lowest of these bands should be 
good approximations to the orbitals occupied when extra electrons are 


added to the crystal. Since there are unlimited numbers of these ~ 


unoccupied orbitals available, it will always be possible to select those of 
the correct symmetry to transform into equivalent orbitals but there may 
well be others, not so selected, whose energies are lower. There is some 
justification, however, for believing that these antibonding equivalent 
orbitals should have low energies. Their form is determined by the form 
of the effective potential in the self-consistent equations and this is com- 
pletely defined by the occupied orbitals. Outside the atomic cores its 
most important part is the exchange potential. For these antibonding 
equivalent orbitals, this potential does provide a Coulomb attraction 
towards the nearest nuclei but it is rather weaker than for the occupied 
orbitals (Hall 1956). Furthermore, the unoccupied orbitals are higher 
solutions of the same equation and so will have more nodes and be more 
diffuse than the occupied orbitals. This reduced localization will mean 
that the interactions of higher order neighbours will be of greater 
importance and, hence, that the shapes of the bands will vary more from 
substance to substance. This general conclusion is well borne out both 
by experiment and by calculation. 

If the phases of the antibonding equivalent orbitals are adjusted so that 
the four which meet at any one atom all have the same sign, then it is 
easily shown that they can be treated exactly like the occupied equivalent 
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orbitals. The energies of the unoccupied crystal orbitals, for example, 
are given .by equations identical in form with those given above except 
that the parameters a, b, c, d must now refer to the interactions of the 
antibonding equivalent orbitals. Values for a, b, c, d have been found, as 
before, by a least squares fitting of the formulae to the energies calculated 
by the same authors and are given in table 3. 

Because of the uncertainties in the original calculations these results 
are not reliable. The values of ¢ and d show, too, that further parameters 
would have to be included before the fitting became satisfactory. A 
positive value of a implies that the orbitals are already in the continuum. 
The antibonding equivalent orbitals will probably, therefore, be little 
more localized than if they had been constructed from plane waves. 
Since the calculated bands for silicon are certainly too low, the true value 
of a for silicon is probably positive also. 


Table 3. Parameters Determined by Fitting Calculated 
Conduction Bands (in Rydbergs) 


Diamond Silicon Germanium 
a 0-106 —0-070 0-538 
b —0-119 —0-010 —0-038 
c —0-072 0-025 —()-025 
d 0-103 0-015 0-005 


§ 6. COMPARISON WITH THE LINEAR COMBINATIONS OF 
Atomic OrspitaLts MEetTHop 


The method of equivalent orbitals, as given in previous papers, has 
been criticized adversely by Slater and Koster (1954). Unfortunately, 
most of their remarks are based on the assumption that an equivalent 
orbital is defined as the sum of two overlapping atomic hybrids one from 
each atom of the bond. Such a definition would indeed lead to an over- 
simplified treatment of the electronic structure and would deserve most 
of their comments, but it is not part of this theory. The equivalent 
orbitals are defined as transformations of the best possible crystal orbitals 
and so must give the most accurate description of the crystal which can 
be given using a determinantal wave-function. The expansion of 
equivalent orbitals in terms of atomic orbitals is not yet known but would 
probably be much more complicated than Slater and Koster have assumed. 

The comparative simplicity of the equivalent orbital method makes it 
difficult to realize the number of difficult theoretical problems which it has 
solved or evaded. Some of these can be seen by comparing it with the 
linear combinations of atomic orbitals method, as applied by Slater and 
Koster to valence crystals. The use of atomic orbitals as a starting point 
might appear to make the method more concrete and easier to visualize 
than the equivalent orbital method, but, in fact, it leads to great difficulties. 
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In the first place, the best atomic orbitals to use for the solid are not known 
but it is clear from the example of the hydrogen molecule that the best 
orbitals for the atom need considerable modification, especially in scale, 
when a strong bond is formed. In the second place, the atomic orbitals 
on different atoms are not orthogonal and it is necessary either to include 
the overlap integrals or to transform to orthogonal orbitals. The inclusion 
of overlap would lead to very different equations from those used by 
Slater and Koster and to a difficult problem of finding independent 
parameters. The transformation to orthogonal orbitals is easier to carry 
out but has two disadvantages. It involves a delocalizing of the atomic 
orbitals and thus blurs the simple picture of atomic orbitals and makes 
necessary the inclusion of matrix elements between second and perhaps 
more distant neighbours. It also involves the danger of transforming 
away all the binding. In the hydrogen molecule, for example, the use of 
orthogonalized atomic orbitals fails to give any binding energy (Slater 
1951). The reason for this is that, in a valence crystal, though not in an 
ionic solid nor in the standard-excited-state treatment of a metal (Hall 
1953), the transformation to orthogonal atomic orbitals involves a 
genuine change in the value of the wave-function and, hence, in the form of 
the self-consistent potential. Some of these difficulties are masked in the 
Slater and Koster theory by the fact that they have not considered the 
self-consistent potential but treated it as independent of the orbital. In 
the accurate theory, the potential depends on the orbital so that the 
equations determining the coefficients of the atomic orbitals in the crystal 
orbital are cubic and not linear as they have assumed. A first approxima- 
tion to the solution of the cubic equations can be found by linearizing 
them, as for example in the standard-excited-state method, but the 
energies of the crystal orbitals are no longer the eigenvalues of the linear 
equations but require various corrections. 

The situation is, therefore, that the use of atomic orbitals involves an 
overlap problem, a problem of solving cubic equations and a problem of 
defining, even in principle, the best atomic orbitals. Also there is no 
reason, even then, to believe that the crystal orbitals can be expressed 
exactly as a linear combination of atomic orbitals, using a finite number 
on each atom. These problems are all solved in the equivalent orbital 
method. The orthogonalization of the equivalent orbitals leaves the total 
wave-function unchanged, the equations to be solved are rigorously 
linear, because of the definition as unitary transformations of the crystal 
orbitals, and the resulting crystal orbitals are exact. In practice, of 
course, the method involves approximations, due to the use of a finite 
number of parameters, but this is involved in any fitting method and 
refinements can be made as more accurate information becomes available. 
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ABSTRACT 


The twist due to a screw dislocation parallel to the axis of an isotropic 
cylinder of arbitrary cross section can be found from the solution of the 
ordinary torsion problem for the same cylinder. Some particular cases 
are worked out. The results are also valid for certain kinds of anisotropy. 


§ 1. INTRODUCTION 


Iv has been pointed out (Mann 1949, Eshelby 1953) that when a screw 
dislocation is introduced along the axis of a circular cylinder the cylinder 
suffers a uniform twist equal (in radians per unit length) to the Burgers 
vector of the dislocation divided by the area of the cross section of the 
cylinder. Recently twists have been observed and accurately measured in 
crystal ‘ whiskers ’ whose cross sections are not circles but, for example, 
regular hexagons or rectangles (Webb et al. 1957, Webb and Forgeng 1958, 
Dragsdorf and Webb 1958). 

In $2 we give a general method for finding the twist due to a screw 
dislocation in a cylinder of general cross section. The dislocation line is 
parallel to the cylinder axis, but its position in the cross section is arbitrary. 
The problem reduces to the solution of the ordinary torsion problem for the 
cylinder. Some particular cases are discussed in §3. These results refer 
to an elastically isotropic material, but in § 4 it is shown that they can also 
be applied with little or no modification to anisotropic crystals in a number 
of particular situations which are likely to occur in actual whiskers. 


§ 2. SOLUTION OF THE Twist PROBLEM 


Let & denote the whisker cross section and C the curve bounding it, and 
let a(x, y) be the twist produced by a screw dislocation parallel to the 
whisker axis and passing through the point (a, y). In order to calculate « 
we give the whisker an arbitrary additional twist «, by applying couples to 
its ends. 

To obtain the stresses associated with «, we have to find a function % 
satisfying 


V2) 0 ane is b= 4(x?+y?) on C. Sa eee tl 


+ Communicated by the Author. 


Oe 
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(See any book on the theory of elasticity : our notation agrees with that of 
Sokolnikoff 1946.) In terms of the auxiliary function 


Pim Yared Oe aee e ee ns () 
the non-vanishing stresses are 
Top Or (OU, . Ty —yor ies oo... (3) 


where , is the shear modulus. 
If the dislocation has a Burgers vector 6 these stresses exert a force 


F,=6t,,= — pba, OF /dx 
F,= —br,,= — pba, oF /dy 


on unit length of it. Thus ¥ is a potential whose gradient determines the 
force on the dislocation. (Compare the analogous role of the flexure 
function in a bent whisker (Eshelby 1953).) 

The work done by the applied couples (per unit length of whisker) when 
the dislocation at (x, y) moves right out of the whisker along any path s is 


| (F,dx+F,dy)= —bayp oF as= — bay F(a, y) 
since Y=0o0n C. But this work may be calculated directly: if D is the 
torsional rigidity of the whisker it is simply —«,«(x,y)D, the work done 
(per unit length of whisker) by the couple «,D as the twist due to the 
dislocation falls from its original value «(x, y) to the value zero when the 
‘dislocation reaches C. Thus ba, =«,«D and we have our basic result 


a(x, y) = F(x, y)ub/D. iay! le ance ie (4) 
Since 


we have 
/ 
x(a) = WE, 9) || it dyteght eke 2 (B) 


The integral of » over Ris $b. This implies that a uniform distribution of 
small dislocations of total strength b per unit area spread over the cross 
section produces a total twist $b, in agreement with a result in the con- 
tinuum theory of dislocations (see, for example, Nye 1952, Bilby et al. 
1958). 

This physical argument (which is equivalent to an application of 
‘Colonnetti’s theorem), is, perhaps, formally unsatisfactory in that we have 
in effect subtracted the infinite self-energy of the dislocation from the two 
sides of the energy balance — «,«D = —ba,¥, and so we sketch an analytical 
proof which also suggests a way to determine « when the ‘’-function is not 
available. It is adapted from the method used in a previous paper 
(Eshelby 1953). We imagine that we have obtained a solution of the elastic 
equations representing a dislocation in a whisker whose side-surfaces are 


P.M. 2H 
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stress-free, and in which the displacement w is required to be everywhere 
parallel to the whisker axis. In this state of anti-plane strain the whisker 
is not twisted but the stresses on its end-surfaces have a certain moment VW 
about the whisker axis. If the necessary end-couples are not supplied 
the whisker will evidently have a twist —M/D. 

Suppose that the dislocation is at r’ =(«’,y’). Then the displacement is. | 
harmonic, behaves like (b/27) tan [(x—«’)/(y—y’)] near r’, and, since the 
sides of the whisker are stress-free, its normal derivative dw/dn vanishes. 
on CO (compare (6) below). It is most convenient to use the conjugate 
function y related to w by 


—= = — S — = ; AP oe Ak: 6): 
Ox ay Tisai an Poel ft ( ) 


If s denotes arc-length along C, (6) shows that dy/0s = — dw/dn, and so on Cy 
has a constant value which we may take to be zero. Further, y is harmonic 
and behaves like (b/27)In|r—r’|nearr’. More compactly, x is defined by 
the requirements: 

Vy=bo(r—r’) int, y~=0 ond. Ger ee 


This implies that — 2zy/b is the Green’s function which solves Dirichlet’s 
problem for the contour C. That is to say, if f(r) is some function defined 
only for values of r on C then the expression 


BO) = Ff fs alee!) date Cae aS eee 


regarded as a function of r’ is harmonic and takes the value f(r’) on C. 
The end couple required to maintain anti-plane strain is given by 


M = | @ — YT xq) du dy = a (25 —u ge) dey a HOM 


or, in terms of y, 


M=-—ty | (grad r?) . (grad y) dx dy 
YR 


=| PV2y dx dy — 4 | polkas | sae 
R c On 

by Green’s theorem. According to (7) the first term is $ubr’2. Equation 
(8) implies that the second term is harmonic in R and reduces to — 4ybr’? on 
C; it is therefore identical with —ubi(x',y’)(eqn. (1)). Thus, dropping 
the primes, M(x, y) is identical with —p»bW (x,y). To find the twist when 
the end couples are removed we divide —M by the torsional rigidity D and 
so recover (4). 

The sign convention we have used for the Burgers vector gives positive: 
twist for positive b. The actual relation between Burgers vector and 
direction of twist is most simply stated thus. Suppose that there is a set. 
of crystal planes perpendicular to the axis of a cylindrical whisker. When 
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a screw dislocation is introduced they join up to form a spiral staircase. 
A creature ascending this staircase would follow a helical path with a certain 
sense. ‘The helix formed by a generator of the undislocated cylinder has 
the opposite sense, and, naturally, a much coarser pitch. 
When *' is not easily available it may be practicable to evaluate M by 
mapping the whisker cross section conformally on to the unit circle. 
Suppose that the dislocation is at the origin and that the relation 


Re ey TCO eee. oa ry leet (LE) 


maps the whisker contour C on to the circumference y of the unit circle 
|¢| <1 in the complex ¢-plane and in such a way that z=0 maps into 
€=0. (The complex variable (11) is, of course, not the same as the z 
which appears as a suffix in (3).) We shall use a suffixed ¢ to distinguish 
quantities taken in the ¢-plane from those in the z-plane. Let us associate 
with each point ¢ the value of x at the corresponding point in the z-plane. 
We obtain a function y.(2) = x[2(¢)] which is easily seen to be the y-function 
for a dislocation at the centre of the unit circle y, that is, 


X¢(¢) = (6/27) In |]. 
It follows from the theory of conformal mapping that 


Ox, Oxe 
pes an ds, 


if ds, ds, are corresponding elements of C and y. Since dy,/dn,=pb/27 
eqn. (10) with r’=0 may be written as 


| |= | BD 
4dr), 
If the transformation (11) can be expressed in the form 


flO= Sa, 


n=1 
then 
|M f= FI > DY din Un* exp {i(m —n)O} dd 
Air ym n 


since (=exp28 ony. Thus we have the simple result 


|M[=td laPo. 2. 0) 


and «=| |/D. 

The point z = 0 can be chosen at will within C, and so our assumption that 
z=f(£) maps z=0 into ¢=0 does not in principle limit the position of the 
dislocation in the whisker cross section. However, there will usually be 
one position of the origin for which f(¢) takes on a manageable form, and 
this may not coincide with the position of the dislocation. We should then 
be faced with the general problem of determining ; or Y’ at an arbitrary 
point within C using, for example, the analysis of Morris (1940) (see 
Sokolnikoff 1946) of which the above is a simple special case. 


2H 2 
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§ 3. Some EXAMPLES 
It is convenient to express the twist in the form 


Burgers vector 6 


twist ¢=K &_ ———— mm :- 
area of whisker cross section A 


We have, for example, «= 1 for any ellipse, « = 10/9 for an equilateral tri- 
angle, « = 1-048 for a square, « = 1-015 for a regular hexagon. In each case 
the dislocation is supposed to be at the centre of the cross section. The 
Y-functions for ellipse and equilateral triangle are simple expressions, 
(Sokolnikoff 1946, §36, §37) and for these and some other cross sections 
there is no difficulty in finding « for a dislocation at any point in the cross 
section. The value for the square follows from the discussion of the 
rectangle below. In the case of the hexagon « was calculated from (12) 
using the transformation 


= ef dt pits . / 6n4+1 
2(¢) = . a — 18)1/3 =C+ ee Ak 
with 
Fine. 2.8.14... {2+6(n—1)} 
aS n! 6”(6n+ 1) 


which transforms a regular hexagon of side a=2(1) = 2792 3-27 [T(4) }? 
into the unit circle (Kober 1957). For the torsional rigidity the value 
D=1-0359 pa* quoted by Polya and Szegé (1951) was used. 

In the case of a rectangle there are rapidly-converging series for Y(x, y) 
and D. (see, for example, Sokolnikoff 1946, §38). The following values 
refer to a dislocation at the centre of a rectangle of sides / and h <1. 


Ih 1 15 2 3 5 10 0 
2 1-048 1-026 0996-0948 0-858 0-800 8 
If we put 
(0, 0) = 4— (8/73) sech (zl/2h), 


i 
+ 
D/ pth? = 3 — (64/7°) {tanh (al/2h) + 3-5 tanh (371/2h)} 
the error in 


__ F(0,0)/h? 
D]ull® 


will not exceed $%. Ifl>3h we have 


k= 3/(1—0-630h/1) 
to better than 4%. 
The twist due to a dislocation not at the centre of the cross section of the 
rectangle is easily found from the series for ¥"(w,y). Asanexample, suppose 
that in a rectangle of sides 1 unit by 2 units the dislocation is displaced 


from the centre to a point equidistant from three sides. Then x is reduced 
from 0-996 to 0-850, or by about 12%. 


————— 
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For any whisker cross section C a well known analogy enables us to 
visualize roughly how the twist varies with the position of the dislocation. 
A hole having C for boundary is cut in a plate and a membrane is stretched 
across it. Ifa uniform pressure is applied to one side the deflection of the 
membrane at (a, y) is proportional to ¥'(a, y) and hence to « or x. There 
is also a membrane analogy for the function y of §2. Let the uniform 
pressure be removed and instead let the membrane be deflected by a single 
point-force applied at a point corresponding to the position of the disloea- 
tion. Then the deflection at (x, y) is proportional to y(x, y). The result 
following eqn. (5) says, in effect, that deflection by a multitude of point- 
forces and deflection by a pressure mean the same thing. 

The expression for the twist needs modification if the dislocation has a 
hollow core. Suppose that this takes the form of a circular cylinder of 
radius 7) small compared with the dimensions of the cross section. In the 
state of anti-plane strain described by y (§ 2), there is no stress across a 
curve y=const. Around the dislocation these curves will be very nearly 
circles, and so the elastic field will not be upset by the removal of the core, 
but | JZ | will be reduced by $ubr,?._ Thus (4) must be replaced by 


a(x, y) = pb (E(x, 4) — bre? }/D. 


Strictly, D as well as | 1/ | is decreased by the presence of the hole, but it is 
not so easy to calculate the correction. Ifthe dislocation happens to be at a 
centre of high symmetry of the boundary curve C the curves ‘Y= const. 
around the dislocation will be nearly circular, and D is reduced by ir7p7,4. 
Whether it has this simple form or not the correction to D will be negligible 
compared with the correction to | J |. 


§ 4. AntsoTROPIC WHISKERS 


If the whisker axis coincides with a four-fold or six-fold axis of the 
crystal the equations governing anti-plane strain and torsion are identical 
with those for an isotropic medium with »=c;, (Voigt 1929, Hshelby et al. 
1953) and our previous expressions for the twist can be used without 
modification. 

Again, suppose that the whisker axis coincides with the z-axis of the 
erystal and that the elastic coefficients C14, C15, C1, Coa, Co5> Cog AN Cy; Vanish. 
Then it is possible to have a state of anti-plane strain with the displacement 
parallel to the whisker axis. The governing equation is 


Grae OT, : Ow aw 
& SS th a Ore 2 > =C © 
- ae By 0 with ry, = C55 Aa Gye 24k oy 


We need to find a solution representing a screw dislocation and from it we 
have to calculate the end-couple required to maintain anti-plane strain. 
Suppose that the dislocation is at (x’, y’) and that the boundary of the 
whisker cross section is the curve 


C:2=F(y); 


forexample, «= + (a?—y?)!? when C isa circle. 
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Introduce the variables 
E=(Cqg/[C55)"?2,  n=Yy. 2 al) ¥ os) ae 


When C is described in the (x, y) plane the curve 


Ds €= (Cqq/C55)"? F(n) 


is described in the (£, 7) plane. It is easy to show that the substitution (13) 
reduces the anisotropic problem to the corresponding isotropic problem in 
the sense that we have to find a function w to satisfy 0?w/dé? + 0?w/dn2=0 
with the boundary conditions w~ (b/27)tan~ [(€—€')/(7—7’')] near 
(é’, n’) and dw/dv=0 on I’, regarding &, 7 as ordinary Cartesian coordinates. 
(visthe normal to I.) The integral (9) translates into 


iy os 0 
Mess | (ne -65) dé dn 


taken over the interior of [. This is the same as JM for a fictitious isotropic 
whisker bounded by [ and having a shear modulus p=c;;. On the other 
hand, it can be shown (Sokolnikoff 1946, §51) that the torsional rigidity 
of the anisotropic whisker is (c;;/c,,)!/? times the torsional rigidity of the 
fictitious whisker. Thus the twist in the real whisker is equal to the twist 
in the fictitious whisker mutiplied by the factor (c¢,,/c;;)"*. This factor 
can be eliminated from the expression for the coefficient « of §3 (the 
fictitious whisker has (c,,/c;;)1/? times the area of the real whisker) and we 
are left with the following simple result : 


Let the whisker axis coincide with the z-axis of a crystal in which 
the elastic constants C14, C15, Cig, Co4, Co5, Cog aNd Cy; vanish. Let the 
dislocation be at the point x= X, y= Y within the whisker contour 
x=F(y). Then the coefficient « is the same as « for a dislocation at 
& = (C44/¢;5)/2X, y= Y in an isotropic whisker bounded by the curve 
© = (Cqy/C55)"" F(y). 


In addition to the more obvious cases this also applies when the whisker 
axis is in the basal plane of a hexagonal crystal. (We use ‘ hexagonal ’ in 
the narrower sense, to denote a crystal which has only five independent 
elastic constants.) ‘The direction of the axis in the basal plane is immaterial. 
What we have called c,, and c;; would in this case be the quantities normally 
tabulated as c,, and cg, respectively. 

We may use our previous results for the isotropic case when the trans- 
formation (13) alters the shape of the whisker cross section but not its 
type. For example, an ellipse becomes an ellipse, and the rectangle 
x= + 3h, y= +31 becomes the rectangle x= + th (cq,/c55)"?, y= +41. On 
the other hand a rectangle whose sides are not parallel to the coordinate 
axes becomes a parallelogram, a figure we have not treated. 
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ABSTRACT 


A theoretical study has been made of the structure of strong ‘ hydro- 
magnetic ’ waves which are propagated, across a magnetic field, in a low 
density plasma where collisions can be neglected. Under these conditions the 
ions are accelerated in the direction of propagation and then brought to rest 
again. The thickness of the wave is determined by the characteristic distance 
/(me2/4ane?) and the wave velocity lies between the Alfvén speed and twice 
the Alfvén speed. 


§ 1. INTRODUCTION 


Astrom (1950a,b) has discussed the propagation of hydromagnetic 
waves through a low density plasma, where collisions are unimportant. 
The direction of propagation was taken at an arbitrary angle to the magnetic 
field, but the treatment referred to small-signal waves, i.e. it was assumed 
that the magnetic field was unaffected by the induced currents flowing 
in the plasma. Ferraro (1956) has described the propagation of strong 
collision-free waves which travel along the magnetic field; this type of 
wave is the collision-free analogue of the Alfvén wave (Alfvén 1942, 1950). 
The analysis presented in this paper is complementary to that of Ferrarof, 
and pertains to the structure of strong collision-free waves which are 
propagated across a magnetic field. 

The physical situation considered in this paper is that of a disturbance 
travelling into undisturbed plasma. Figure 1 is a schematic diagram of 
the disturbance (or wave) showing the trajectories of the electrons and 
ions ; it is convenient to use a coordinate system moving with the wave, 
so that the plasma appears to flow through the wave, from left to right. 
The energy associated with the wave is partly kinetic energy and partly 
magnetic energy. It is assumed that the energy which a particle acquires. 
in the wave is much greater than its original thermal energy, so that the. 
latter can be ignored. The increase in the magnetic field strength is due: 
to the transverse electron currents shown in fig. 1. A small difference 


+ Communicated by the Authors. 

{It was stated in Ferraro’s paper that the waves are ‘ hydromagnetic ’ if 
the frequency is small compared with the plasma frequency and the velocity 
is small compared with the speed of light. The correct limit on the frequency 
is that it should be small compared with the gyrofrequency of a positive ion in 
the guiding field. 
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between the electron and ion densities produces an electrostatic field in 
the direction of propagation. 

In the laboratory frame of reference the massive positive ions are 
accelerated by the electrostatic field, and then brought to rest again. 
Both the ion density and magnetic field increase within the wave, and 
then return to their initial values. The plasma returns to its initial state 
after the passage of the wave, except that each particle has been displaced 
in the direction of propagation. The thickness of the wave is determined 
by the characteristic distance 4/(mc/4ne?) and the wave velocity lies 
between the Alfvén speed and twice the Alfvén speed, the velocity increasing 
with the strength of the wave. 


Fig. 1 


ELECTROSTATIC rates 
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Schematic diagram showing the trajectories of the electrons and ions, viewed 
in a coordinate system moving with the wave. An induced electric 
field exists in the y-direction. 


§ 2. FUNDAMENTAL HQuaTIONS 


The coordinate system is chosen so that the magnetic field is in the 
z-direction. The case considered is that of a plane wave with no variations 
in the y-direction or the z-direction. Thus the equations of motion for 
the electrons and ions can be written as follows, where suffixes 1 and 2 refer 
to the electron and ion respectively, and wv and v are the velocity components 
in the a- and y-directions. 


Ou, 8 “| Bi aerate ss! (1) 
ey wera Pasig : , 


Coates Ot, era: E -4=| hea a hetie)' 
Br hae aa ae y i 
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back rece eee me Vee | ge en ae 
ot + MEG Mz ot c (3) 
OV, OU, ~~ ey B Ue Ds ys 
OL i a Oe Ma\|. * C (4) 
The field equations, neglecting the displacement currents, are 
4 
e =— a [142,04 + Melee] = & «= « «se 
EH, 4 
= oe oy eee ee a) 
cand 
aby __ 128, a 
Ox c dt 


The numerical values of » and k (the permeability and dielectric constant 
of a vacuum) are both unity, but they are retained in the equations. 

These equations are much simplified by the use of the fundamental 
property of a plasma, namely that |n, — | <n, (Tonks and Langmuir 1929). 
‘Thus both n, and n, will be replaced by n in all equations except eqn. (6); 
the difference between n, and n, produces the electrostatic field H, which 
plays an important role in this type of wave. It has been assumed that 
the ions are singly charged, so that e, = — ey. 

Only waves of constant profile will be considered in this paper, and it 
is convenient to use a coordinate system moving with the wave, as already 
mentioned. The undisturbed plasma has a density n) and velocity uw, in 
this coordinate system, and is situated in a uniform magnetic field By 
at large negative values of x (see fig. 1). As the charge separation is small, 
it follows that u,=wu,=u, to a high degree of approximation. 

The relative velocity between the electrons and ions will be denoted 
by v, i.e. v=(v,—v,) and the suffix z will be dropped from B,, since the 
magnetic field has no other component. Eliminating #,, between eqns. (1) 
and (3) and putting 0/dt=0 gives 


we = _ vB * (8) 
dx (m,+m,)c 
Equations (2) and (4) give 
le SEP atypia) ES =| Me () 
dx M Mz c 


and eqn. (5) reduces to 
dH — Amngtlyeav 
due a 
SINC — (Nea) Ue 
The electric field in the undisturbed plasma is wu )B)/c, when measured 
in the coordinate system moving with the wave. Equation (7) shows that 
0H ,/dx=0 when 0B,/0t=0, thus EL, has the constant value u)B,/e for 
all values of a. 
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§ 3. SOLUTION OF THE EQUATIONS 


The dimensionless quantities x,=2/d, u,=ulv4, «= His, Vy = Clee 
and Bb, =B/B, where 


2 AI ergo Fm wale = al Lat FD 
Arrnopeg?(m + My) mau ae + ima | 
By?(m, +My) +g) 
ene 


will be introduced so that eqns. (8), (9) and (10) can be written 


and 


du 
LS B 
Un ee Un n? 8 a) 
AV, 
NP Re a (a UnBy) (9 a) 
and 
dB, AV» 
dae Te (10 a) 
Eliminating v, between ( 8 a) and (10 a) gives 
UN 
ppNRC ee ee 
Integration gives 
Ve Beals oye = eke ae geees Be C9) 


where the constant of integration is chosen so that wu, =a when B,=1 
Eliminating x, between eqns. (8 a) and (9a), and substituting for w,, 
and du,/dB,, using (11) and (12) gives 


peat 1+ | \ B,2—-2B +C 
where C is a constant of integration. In the undisturbed plasma, 6, =1, 
V»,=0 so that C can be evaluated and the resulting expression for v,, is 


Raia (eee Boe ei (13) 


Substituting for w, and v, in eqn. (10 a), using (12) and (13) gives 
aB,, B,-1)V/[(2«+1+B,)(2«—-1—B,)] 
dx, a 7 2a? + 1 —B,? 
which shows that B,, has stationary values at B,=1 and 5, =(2«—1). 
The positive sign in (14) will be used for the present, and the reason 


for this choice will appear later. 
Inversion and integration of eqn. (14) gives 


(14) 


a 


ie n= — vil Qat1+B,)2e—1— B+ Fea 


(202—1—B,,) — +/{(a? — 1)(20+14+B,,)(2« —1—B,,)} 
«In| “(B,— 1) {et ue 
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The constant of integration has been chosen so that x is measured from the 
point where B,, has the stationary value (2a —1). 

Examination of (15) shows that B,>1 as x, — 00, which is the required 
boundary condition for the undisturbed plasma; this justifies the choice 
of the positive sign in (14) since the negative sign leads to the result B,>1 
asa,>+ 00. In order to determine B, for positive values of x,, however, 
the negative sign of (14) must be used, since the positive sign yields no. 
further information. The resulting curve is an even function of w,, and 
is plotted in fig. 2 (a) for different values of «. 

Figure 2 also contains plots of the electric field (fig. 2 ()), the relative 
velocity (fig. 2 (c)), and the ion density (fig. 2 (d)). The relative velocity v 
was calculated using (13), and the ion density nyu,/u was obtained from (12). 
The electric field was eliminated from the initial equations, but it is easily 
shown that 


B= — (Mae pe Se ee 
My+M,/) ¢ 
or, using dimensionless quantities, 
EL, = —0,B,, iy oioo eae ig an 
where the unit of electric field strength is 
(“2 a m1) VBo 
M,+m,) Cc 


The range of values of «, the parameter denoting the strength of the 
wave,isl<a<2. This can be shown as follows: substituting B, = (2«—1) 
into eqn. (12) shows that w,, =(2—«) ata, =0. Waves containing multiple 
streams of ions are not considered in this paper, so that w, is always positive 
and «<2. Equation (13) shows that v,>(B,—1){(a+1)(«—1)}?/« 
as B,>1, so thata>1. Thus 1<a<2, ie. v, <u) < 2v, for a wave of this. 
type. 


§ 4. EXAMINATION OF THE PLASMA PROPERTY |n,—7,|<7, 
The net charge responsible for the electrostatic field E, is described 
by eqn. (6), viz. 


dH, 4, 
a yp om te 2) ns Ree eee ae 


This quantity is greatest at the centre of the wave, as shown by fig. 2 (6); 
furthermore it increases with the value of x. Differentiation of eqn. (16) 
gives dH’,/dx, and at x =0 (where dB/dx=0) 


(=) _ _ (M@z—m,\ Bdv 
da }2=0 iI) Cb Oe oan is 


In the strongest waves, in which « is slightly less than two, L>uBle 
at «=0. Using this condition eqn. (9) shows that 


(2 _ __ &s(y + Mg) Uy By 
dx} x=0 UMMC 


since £,, =u )Bo/c. 
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The variation throughout the wave of (a) the magnetic field, (5) the electrostatic 
field, (c) the relative velocity between electrons and ions, and (d) the 


ion density. 
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Substituting this expression in (18) and combining the result with (6) 


gives 
My—Ny\ _ (My—M (2a — a) can Oem 
n Ms 47nym,c* 
since B=(2«—1) B, and nyuy=nu. As a>2 
Ny — Ny My —™M4 Vd oP 20 
( ” ) Fa ( Me baa 7) 


which can be written 


My—M\ 4 (M2— M1 pkvy? 
n Mz+mMy, oe A 
Thus the assumption that |n,—7,|<n, is justified if (v,/c)? <1. 


§ 5. DIscussION 


Three conditions must be satisfied if waves of the type illustrated in 
figs. 1 and 2 are to be excited. 


(a) (vy/c)?<1. It has been shown in §4 that this condition must hold 
if |ny—n,|<n; it is also the condition for non-relativistic equations to be 
valid. An alternative form is By?/47n)u<m,c?. 

(6) kT,<B,?/47n,u. The initial random energy is negligible if this 
condition holds, since the electrons are given a velocity of the order of v,. 
Combination of (a) and (b) gives kT, < B,?/4anyu<m,c?. 

(c) (d/v4)<(npvQ)-1 where v is the random velocity of the electrons 
and Q is the positive ion cross section for momentum transfer. This is the 
criterion for collisions to be unimportant, since the time which an ion 
(or electron) spends in the wave is ~ (d/v,), and the mean time between 
collisions is (nvQ)-. 

An interesting result of the present work is the emergence of the charac- 
teristic distance d(==5 x 10°n—?.em). It is identical with the collision-free 
‘skin-depth ’ of a plasma, for frequencies well below the plasma frequency 
(Spitzer 1956). 

A strong wave imparts an energy of 4(B,?/87n)u) to ions which it 
encounters. ‘This energy may be appreciable, for example if n»=10% 
ions/em® and B=1000 gauss the ions are accelerated to an energy of 
10% electron volts. Such waves may be excited if a plasma containing 
a magnetic field is rapidly compressed ; the energy given to the plasma at 
the boundary is conveyed to the interior, partly as magnetic energy, and 
partly as kinetic energy of organized motion. In certain cases, some of 
this energy will eventually be degraded into thermal energy, by coulomb 
collisions or collisions with neutral atoms. 
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Note added in proof.—The characteristic distance d is an important 
parameter in Rosenbluth’s recent theory of the pinch effect. 


RosENBLUTH, M., 1957, Magnetohydrodynamics, ed. R. K. M. Landshoff, p. 57, 
(Stanford). 
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Tue theoretical prediction by Lee and Yang (1956) that parity might not 
be conserved in B-decay has stimulated an interest in the study of the 
angular distribution of B-rays from polarized nuclei. This was further 
increased by the experiment at the National Bureau of Standards (Wu 
et al. 1957), which first confirmed this hypothesis. As a result we were 
led to consider possible B-emitters which could be polarized by the low 
temperature method. One of these, Pr seemed to be a favourable 
isotope since it is a pure B-emitter and if the B-intensity is recorded by the 
blackening of a photographic film (Grace ef al. 1957) the absence of 
y-radiation could be an advantage. To test whether polarization of 
143Pr was likely to be obtained we have carried out experiments on the 
neighbouring 19hr isotope !42Pr whose orientation can be detected and 
measured by observation of its y-radiation. 

In earlier experiments (e.g. Ambler et al. 1956, Cacho et al. 1955) isotopes 
of some rare earth elements have been oriented in both the double nitrate 
and ethyl sulphate salts. Paramagnetic resonance experiments on Pr?* 
(Baker and Bleaney 1955, Duffus 1954, Cooke and Duffus 1955) show that 
the lowest electronic level is a magnetic doublet in both salts and that the 
observed spectrum can be represented by the spin Hamiltonian, 


we = 9, PH S,+ ASI, FAS. + ADS 35 @r eee ae (1) 


with an effective spin S=4}. It is seen that the hyperfine structure is 
is completely anisotropic and therefore axial alignment should arise in 
zero magnetic field when the crystal is cooled to sufficiently low temper- 
atures. The quantity A= /(A,?+A,?) is not sharply defined but has a 
two dimensional Gaussian distribution of standard deviation A centred on 
zero, and has been ascribed to distortions of the local crystal field (Bleaney 
and Scovil 1952); its magnitude is difficult to measure accurately from the 
paramagnetic resonance spectrum and may be expected to depend 
sensitively on concentration and nature of the constituents of the crystal. 
The measurements on ™!Pr in lanthanum magnesium nitrate at 4°K 
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show that A =0-077em-, 9, = 1:55 and A is about 0-lem— (Cooke and 
Duffus 1955, Duffus 1953). For !!Pr in lanthanum ethyl sulphate at 
20°K, A=0-075 cm—, 9, = 1:52 and A is about 0-04em-! (Baker and 
Bleaney 1955). 

For any particular ion, the energy levels are given by 


W= +}V/{(9,8H,+Am)?+ A*}, 


where m is the nuclear magnetic quantum number. As A increases, the 
different m-states are crowded together so that the degree of alignment 
at any given temperature is reduced. The effect of A can be removed 
by the application of an external magnetic field such that 9, PH,>A, 
when the separation between the different m-states will approach the full 
value of $A. From the expression for the level positions the relative 
population distribution of the nuclear magnetic states is easily calculated 
as a function of temperature for given values of the parameters A, A and 
g,8H,. The population distribution for a given value of A is obtained 
by a Gaussian averaging over all A. The y-ray anisotropy can then be 
calculated as a function of temperature for any particular decay scheme. 

Since the values of A for ‘Pr are similar for the two crystal lattices we 
decided to incorporate '4*Pr into cerium magnesium nitrate. This choice 
was governed by the following considerations; firstly, temperatures 
down to a few millidegrees can be reached by adiabatic demagnetization o 
cerium magnesium nitrate (Daniels and Robinson 1953); and secondly, 
since g, for Ce ** in this salt is small it is possible to apply the axial magnetic 
field necessary for overcoming the effect of A without producing a pro- 
hibitive rise in temperature. 

The decay scheme of 14*Pr, which is formed by neutron capture in 
stable #41Pr, is shown in fig. 1. (Pohm et al. 1954). The unique first 
forbidden f-transition to the ground state of 14Nd shows that 1%Pr has 
a spin of 2~. The f-transition to the 1-57 Mev level in 1**Nd is classified 
by its comparative half-life as first forbidden with 1,=0 or 1 and hence 
the spin of the 1-57 Mev level is 1+, 2+ or 3+: the choice of 2+ was made by 
Pohm et al. on the basis of the general systematics of the energies of the 
first 2+ excited states of even—even nuclei. This spin sequence has been 
used for the interpretation of the present experiments. 

In the experiment a single crystal of cerium magnesium nitrate (about 
1g) was grown from solution containing Pr and this was mounted in 
the demagnetization cryostat: it was estimated that about 3% of the 
Ce?+ ions had been replaced by Pr?+. The crystal was then cooled by 
adiabatic demagnetization and the y-ray anisotropy «{=1—(0)/I (5)} 
was measured as a function of temperature in zero field. These measure- 
ments were also repeated with axial magnetic fields of 300 and 600 gauss 
and the results are shown in fig. 2. The amount of activity which can be 
incorporated in the crystal is determined by the heating up rate due to 
B-activity and in this case was chosen to be about 50 pc. Because only 
5% of the decays are accompanied by y-rays the counting rate 1s neces- 
sarily small, and gives rise to the large statistical errors. As can be seen the 
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The dependence on temperature (1/7 in degrees K~1) of the anisotropy of the 
1-57 Mev y-ray from Pr in zero field (@) and with axial magnetic 
fields of 300 gauss (O) and 600 gauss (X). For each of the three experi- 
ments, a typical statistical error is indicated, 
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anisotropy at a given temperature is increased by the application of the 
field. This change in y-ray anisotropy shows that the degree of nuclear 
alignment, ¢(J,2), has been increased. Since the nuclear orientation 
arises through a magnetic hfs coupling, any appreciable increase in Cie) 
through the application of a small magnetic field must be accompanied 
by a nuclear polarization (/,). Thus, although we have not performed 
an experiment specific to nuclear polarization, the increase in y-ray 
anisotropy associated with the applied magnetic field implies that a 
nuclear polarization has been produced. The figure also shows that the 
observed anisotropies in 300 and 600 gauss are closely the same. This 


indicates that the effect of A has been largely removed by these fields. 
The sign of the anisotropy rules out the possibility of dipole radiation and 
confirms that the 1-57 Mev level must have a spin greater than one. This 
is in accord with the assignment of 2+ made by Pohm ef al. The temper- 
atures in the applied field have been calculated from the susceptibility 
assuming g,,=0-03 for Ce?+ (Wheatley and Estle 1956). There is an 
apparent contradiction between the relatively large temperature rise on 
application of the magnetic field, and the small value of g,. There are 
two possible explanations ; either the g, value is larger than that given by 
Wheatley and Estle, or application of the field produces irreversible 
heating through some secondary processes. ‘The existence of this rise in 
temperature does not, however, lead to any uncertainty in the inter- 
pretation of the results which depend only on the knowledge that the 
effect of the A term has been removed by the magnetic field. 

From the anisotropy of the y-radiation in the magnetic field the hfs 
constant A and hence the nuclear magnetic moment can be determined. 
Since the angular momentum change (¢,) accompanying the f-transition 
can be 0 or 1, two possible sets of values are obtained. Taking 


e= +0°1404+ 0-015 at 1/7’ = 250 


we find 
Acm>* wn.m, 
1p=0 0-0027 + 0-0003 O-11+0-01 
ip= il 0-0042 + 0-0005 Ont + 0:02. 


With these values of A, the calculated nuclear polarization ({J,/1)) of 
142Pr in fields above 300 gauss and at 1/7’ = 250 are 0-44 (t, = 0) or 0-62 (t, = he 
By comparing the ¢-1/7' curves with and without field we find that A is 
about 0:004em— if i,=0 and 0-006cm™ if 1,=1. The curves shown 
have been calculated using these values of A and A. 

The known spins and moments of !4'Pr and !°Nd are in accord with 
shell model predictions and suggest that in Pr the odd proton and neutron 
are in d5/2 and f7/2 states respectively. The use of free proton and neutron 
moments leads to a value for the nuclear moment of !47Pr of about 1-5n.m. ; 
this is about ten times greater than that observed experimentally. How- 
ever, calculations by Caine (1957) show that the effect of configurational 
mixing is to reduce this value substantially and it may be possible, there - 
fore, to understand the small moment on this basis. 


212 
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From the present experiment, the following conclusions may be drawn. 


(i) A substantial nuclear polarization of ‘Pr has been obtained. 

(ii) The orientation of !42Pr in cerium magnesium nitrate is satisfactorily 
described by the spin Hamiltonian of eqn. (1). 

(iii) The nuclear moment of 14Pr is estimated to be 0-11 n.m. (¢g=0) or 
O-17n.n (tp=1). 

(iv) Since the nuclear moment of 14*Pr is likely to be larger than that of 
142Pr and similar to that of stable 441Pr (u4=3-9n.m.) we can expect to 
obtain a large nuclear polarization in this isotope by the present method. 
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ABSTRACT 


Values have been obtained for the latent heat and saturated vapour 
density of helium between 2-2°xK and 5:2°K, and the critical constants have 
been re-assessed. It is shown how the calculated vapour density may be 
used to examine the consistency of a liquid helium temperature scale. 


§ 1. [IvTRopUCTION 


IN a previous paper (Berman and Poulter 1952—referred to as I), some 
values were given for the amount of heat required to liberate a mole of 
helium vapour from a calorimeter containing liquid helium. The quantity 
actually determined, /,, in such measurements is related to the latent 
heat, L, as usually defined, by the equation 


Det Oe polod\” pattie teen eae) 


where p, and p, are the densities of saturated vapour and of liquid 
respectively. The factor (1—p,/p,) arises because a fraction, p,/p,, 
of the vapour produced occupies the space in the calorimeter vacated by 
the evaporated liquid. 

The measurements reported in I were made in order to calibrate a 
calorimeter in which the volume of gas liberated was used to measure 
enthalpy changes. It was later shown (Berman and Swenson 1954) that if 
a liquid helium temperature scale (the temperature-vapour pressure 
relation) is assumed, then the saturated vapour density can be derived, 
without a knowledge of the liquid density, from the relation 
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The vapour density can also be obtained from the virial equation, or 
measured directly, and a comparison of such values with those calculated 
from eqn. (2) can be used to examine the consistency of the assumed 
temperature scale. We have reported this work briefly (Berman and 
Mate 1957), and in the present paper we describe some experimental 
details and present the results in full. 
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§ 2. THe EXPERIMENTS 


The essentials of the apparatus were described in I, so that here we 
mention only the changes made for the present series of experiments. 
The relevant part is shown in fig. 1, in which 8 is the high pressure pot of 
a Simon liquefier. A is a small annular vessel into which helium can 
be condensed through B, which leads to a helium supply and a burette 
system. F is the dropping tube used in the enthalpy measurements, and 
H is the small heating coil needed for determinations of L,. E is a heat 
exchanger used to increase the efficiency of condensation into A. 


Fig. 1 


F B 


Sketch of calorimeter. 


It had been found that with helium boiling at atmospheric pressure, 
the evaporation from the pot S sometimes became so large that it emptied 
before readings could be completed. At higher pressures this occurred 
more frequently, and experiments in the temperature range from the 
boiling point to the critical temperature were usually abortive. It was 
concluded that the well known oscillations of a helium gas column were 
responsible for these effects, since the onset of high evaporation rates was 
quite sudden and the instability increased with pressure (see ter Haar 
1955). Several changes were made in the apparatus in attempts to 
eliminate the oscillations, which were eventually shown to be associated 
with the heat exchanger, K. The oscillations appeared to arise in the 
tubing B, connected to the tail but their main effect was to convey heat 
to 8, increasing the evaporation from the tail only slightly. Replacement 
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of the heat exchanger by two separate straight tubes to the pot and to the 
tail entirely eliminated oscillations at all pressures. Naturally, more 
helium was evaporated from S during condensation; nevertheless several 
latent heat determinations were still possible for each expansion. Later 
to reduce the background evaporation, the central tube, F (required for 
enthalpy measurements), was cut above S and sealed, preventing radiation 
from reaching the tail; the elimination of conduction along this tube 
halved the evaporation from the pot, more than compensating for the 


loss during condensation. 
§ 3. RESULTS 


All calculations have been made using the temperature scale 7’ 


proposed by Clement (see reference). SBE 
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3.1. Experimental Values of L, 


The experimental results are shown in fig. 2. They extend over a 
greater temperature range than the few measurements reported in I, 
and the values are of the order of 1% higher, the discrepancy being greater 
at higher temperatures. A discrepancy of 0-4% is accounted for by the 
neglect in I of the variation of mercury density with temperature. (The 
pressure of the collected gas was measured with mercury at room temper- 
ature.) Other systematic errors which might arise (e.g. calibration of 
the burette tubes, manometer scales, stop watch and electrical standards) 
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were thoroughly investigated during the present series of measurements, 
and were not found to account for any appreciable fraction of the remaining 
discrepancy. 

At the beginning of the present experiments only one completely 
satisfactory run was made with the apparatus in the same condition as 
for the earlier measurements. This gave a value of L, differing from 
the corresponding value in I only by the mercury correction discussed 
above. After the heat exchanger was eliminated higher values of L, 
were obtained, and at the same time several improvements in the behaviour 
of the apparatus were noticed: The evaporation from A returned to its 
background value much more rapidly after the heater current was switched 
off; also the final and initial backgrounds, which previously had nearly 
always differed, were now usually the same. A constriction in the heat 
exchanger could account for the first observation and the pressure built 
up in A could lead to low values of 1, as well as encourage oscillations. 
As the apparatus can now be used up to the critical pressure with no sign 
of oscillations, it seems safe to assume that these are absent at all lower 
pressures too, and we believe the present, higher, values of L, to be correct. 

It should be noted that the previous values seem to be in error systemati- 
cally, so that calibration of the apparatus for enthalpy measurements is 
not affected. 

The scatter of the experimental points is seen from fig. 2 to be of the 
order of 0-5° over most of the temperature range. At the lower tem- 
peratures (< 2-5°K) the scatter is greater due to the increase in volume of 
the background evaporation (the same mass of gas must be acommodated 
in the burette at a lower pressure) which, even for the largest burette tube 
we had, became a large fraction of the total evaporation measured. At 
the lower pressures involved, fluctuations in pressure also have more effect. 
As the general scatter was rather larger than expected we attempted to 
trace its cause, and looked for a dependence of L,, on the rate of heating or 
on the height of liquid helium in vessel A. In I measurements were made 
at different heating rates, but on examining the results more critically it 
was not considered certain that this had no effect. A variation with rate 
of heating would suggest that the background correction applied was not 
appropriate, as the extra flow of cold gas during the heating period alters 
the temperature gradient in tube B. By varying the temperature differ- 
ence between 8 and A and finding the change in background rate it was 
concluded that an insignificant part of the heat leak flows down B, and 
indeed no change in L,, was found in the present measurements when the 
rate of heating was varied by a factor of 5. If there is insufficient liquid 
in A to extract the heat efficiently, some heat from the heater will pass up 
the metal walls and reach the vapour instead of evaporate liquid. There 
was no change in L,, when the height of liquid was varied by a factor of 20. 

In the course of the work three heaters were used with resistances varying 
from 100 to 600 ohms, and the electrical wiring was altered. In addition, 
runs were made for different lengths of time for fixed rates of heating. 
Again, none of these variations was reflected in the values of L . 
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3.2. The Latent Heat 


Smoothed values of L, taken from fig. 2 are given in the table; these 
should be accurate to a few tenths of 1%. From L, and the T'55, Vapour 
pressure-temperature relation we have calculated p, by eqn. (2). p, 
derived in this way is denoted by p,, in the table. Below 4:4°x the liquid 
densities shown in the table are taken from Kerr (1957). The only data 
for p, at higher temperatures are from the measurements of Mathias et al. 
(1925) and Clement has tabulated values from them by taking account of 
the change in temperature scale. As these values are only 0-25%, higher 
than Kerr’s at 4:4°K, we use them to complete the table above this tem- 
perature without further adjustment. 

These values of Z,, p,; and p, have been used in eqn. (1) to calculate the 
latent heat, L, which is also shown in the table. Vapour densities derived 
from the equation of state are possibly slightly more accurate than our 
values below about 3°K as the non-ideality is small, while our values are 
dependent on a knowledge of dP/d7’. We have, however, used our values 
throughout since at the lower temperatures the error produced in L is 
negligible (even at 3-8°K p,/p, ~ 0-1, so that an error of 1% in p, will only 
affect L by 0-1%), and at higher temperatures, where p,/p, is large, the 
saturated vapour density can not be calculated from a virial equation with 
the present knowledge of the virial coefficients. 

Our latent heats differ from those tabulated as L.... by van Dijk and 


exp 
Durieux (1957) throughout the common temperature range. L,., was 
derived from the experimental data of Dana and Kamerlingh Onnes (1926) 
and Berman and Poulter, and the smoothed values given by Keesom (1942), 
At 2-4°x our values are about 0-5°% higher and the difference increases 
steadily to 17% at 4:2°xK. 
3.3. Critical Data 

In fig. 3 we have plotted the vapour density, p,,, the liquid density, p,, 
and the diameter, }(p,,+p,), using as before our calculated values of p,. 

The critical density obtained from fig. 3 is pg=0-0680 + 0-0003 g/em?®. 
This is in extremely good agreement with the value of 0-0675 + 0-0005 g/cm?® 
which Kerr suggests as most likely, from a discussion of several possible 
ways in which his values of p, can be used to obtain p,. 

Above 5°K the liquid density, as tabulated by Clement, decreases 
rapidly and the critical temperature is fixed fairly closely by the shape of 
this curve. Since we now have, in addition, values for the vapour density 
and can construct the diameter, it is possible to narrow the limits of un- 
certainty of 7. On the 7'5;, scale we find T= 5-194 + 0-002°K and the 
critical pressure P,=172-4+0-3emHg at 20°c. The critical coefficient 
(RT op¢)/Po= 320. Kamerlingh Onnes (1911) determined P,, directly as 
172-4 + 1-0cm Hg, which corresponds on the 7’;,,, scale to 5-194 + 0-:008°K. 


§4. THe Vapour DENSITY AND THE TEMPERATURE SCALE 
The density of a vapour can be represented by a virial equation of the 
form PV2=RTOUSB/V +OC/V2+D/Vi+...) . . . « (3) 
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where the virial coefficients B, C, D ... can be derived from measurements 
of PV as a function of V for fixed temperatures. Keller (1955) has 
measured such isotherms for helium at several temperatures. At 4°K 
the highest vapour pressure used corresponded to about 40°%, of the satur- 
ation pressure, and near 2°K the highest pressure was over 90% of 
saturation pressure. Within the limits of the experimental uncertainties 
Keller found that the measured densities could be represented by eqn. (3), 
using only the second virial coefficient B. We would not necessarily 
expect, however, that the saturated vapour can be represented in the 
same way when the temperature and consequently the vapour pressure 
is high enough for large departures from ideality to exist. Using B alone, 
the saturated vapour at 4-2°K would have a density 40%, greater than an 
ideal gas at the same temperature and pressure, and at 4:4°xK the calculated 
density becomes imaginary. Even at 3-0°K the departure from ideality is 
still 15%. As the saturated density decreases the vapour behaves more 
like an ideal gas, and the small departures should be adequately represented 
by the term in B. 
Fig. 3 


Density, grams/cm 


25 30 35 40 45 50 
ip 

Liquid density, vapour density and the diameter. 
We have calculated p,,, the saturated vapour density derived from 
eqn. (3), when only the first two terms are used. We have taken for B 
the values tabulated by Kilpatrick et al. (1955) which are based on Keller’s 
measurements. Any adjustments of 6 which might be necessary on 
account of changes in the temperature scale have a negligible influence on 
Py Figure 4 shows the difference, Ap=py7p— pj, between the Hens 
derived from eqns. (2) and (3). The rapid increase of Ap above 3-3°K 
merely indicates the inadequacy of using only two terms in eqn. (3) and 
shows that if the third term is responsible then C must be positive. We, 
would expect Ap to approach zero at lower temperatures in a regular way, 
and the irregular shape of the curve found here must reflect some errors in 
the parameters entering eqns. (2) and (3). The irregularities can not be 
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accounted for by the small uncertainties in the values of B, and except at 
the lowest temperatures, where the curve has been drawn as a broken line, 
they are greater than the uncertainties in L,. Fora measured pressure it 
is unlikely that 7’,,, is wrong by more than a few millidegrees, an error 
which could only account for about one tenth of 1% discrepancy. We 
therefore ascribe the shape of the Ap curve to errors in the values of dP/dT 
used to calculate p,,. If one point on the P-T' curve were known, we 
could find and integrate corrections to dP/dT to give corrections to 1’ at 
other temperatures. Failing this knowledge, we have assumed that at 
2-55°k—the centre of the region where Ap is small—7’,,,, is correct. We 
have then calculated corrections to 7';,,, which are shown in the table. 
We do not give values above 3-0°K as it is not known how much of Ap is 
due to the inadequacy of calculating p,,,, from B alone. The table does 


a A 
Fig. 4 


22 2:4 26 28 3:0 32 3-4 36 Se} 4-0 
7 SK 


Variation of Ap with temperature. 

show, however, that between 2-3 and 2-8°xK 7’,;,, can only be in error by a 
constant amount, apart from a variation of the order of 1/10 millidegree. 
Above this range the correction seems to start as a positive quantity, 
relative to any correction at 2-55°K, but we cannot make any deductions 
for temperatures above 3°K for the reason discussed. There is only a 
slight suggestion that below 2-3°xK the correction is in the opposite sense. 

Sydoriak and Roberts (1957) give a vapour pressure equation for ?He 
which fits experimental data between 0:45°xK and the critical point (3-3°K). 
They show a graph of the difference between T 55, and T’,,,, which bears a 


striking resemblance to the corrections we have suggested in the neighbour- 
hood of 2:5°K 
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We will not pursue further the relation of our measurements to the 
temperature scale in view of the proposals to make another attempt to 
produce an accurate scale. It is evident, however, that our results, when 
used in the manner we have outlined, will be useful in examining the con- 
sistency of any future scale. 
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ABSTRACT 


Zener’s theory of the volume expansion of self-strained elastic bodies is 
applied to screw and edge dislocations in isotropic media. The important 
physical parameters entering are the pressure dependence of the shear 
modulus, and, to a lesser extent, that of the bulk modulus. Numerical 
results for a number of materials are given. The volume expansion per 
atomic length of a dislocation line is found to be of the order of magnitude 
of one to two atomic volumes, thus confirming an earlier estimate by Stehle 
and Seeger. 


In 1942 Zener published a theory relating the volume expansion of a self- 
strained elastic body to the stored elastic energy. The present authors 
(Haasen and Lawson 1955, Seeger 1957a) subsequently have pointed out 
that Zener’s theory can be employed to calculate the decrease in density of 
a crystal upon introduction of dislocations. It is the purpose of the 
present note to give the results of such a calculation. We shall compare 
these with a formula recently proposed by Lomer (1957) and shall show that 
Lomer’s result can be obtained from Zener’s more general approach by 
making the following assumptions: the total elastic energy is dilatational 
energy ; there is no energy due to shear strains. 

According to Zener (1942) the average dilation of a self-strained isotropic 
body is given within the framework of second order elasticity theory by 


= 1 /dlogk 1 /dlogG 

Here Kk denotes the (isothermal) bulk modulus, G the shear modulus, V the 
specific volume. Wand W, are the average work done on the unit volume 
in establishing the dilation and shearing strains, respectively. It will be 
demonstrated below that for practical purposes W, and W , can be replaced 
by the average densities of the stored energies of dilation and shear strains. 
The physical basis of eqn. (1) lies in the fact that in a self-strained body W 4 
and W are proportional to K and G, respectively. These elastic constants 
depend on the specific volume. Therefore the free energy of the strained 
body can be lowered by a volume expansion. 

Treating a number of pertinent problems E. Mann (to be published) has 
shown that eqn. (1) gives the same result as Murnaghan’s (1951) theory of 
second order elasticity the application of which, however, is much more 
tedious than that of eqn. (1). . 


y+ Communicated by the Authors, 
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Equation (1) will now be applied to a straight dislocation of strength / 
which is, as usual, assumed to lie in the axis of an isotropic elastic cylinder 
(Poisson’s ratio v) with outer radius R and inner radius r,. In this case it 
is convenient to consider the volume expansion per unit length of the 
dislocation, which is given by 


1 /dK 1 /dG .-G 
Si) ax Ai yh . 
: K (i i) SiG (s z) sg aret gegen 


Tn going from eqn. (1) to eqn. (2) we have used the relation 


dp 
= — Y-— 
Kea) Ri ee Se oe OE) 


between bulk modulus, volume, and hydrostatic pressure p.H, and H, are 
the dilatation and shear contributions to the dislocation line energy. They 
are given by 

E,=9, 


G25) R gp eee! Si iOea eC), 


for a screw dislocation, and by 
1Gb?1—v—2v?, R 
been 4m (=e | : 
_ 26881 48 R [ ee ek) 
* 347 (1—v)? are | 
for an edge dislocationt. As is well known the total line energy of an edge 
dislocation is 
Gog oes ye fon Frid Wes Pa ASD) 
47 (1—v) Yo 
Inserting these expressions for the line energies into eqn. (2) we obtain for 
the volume expansion per unit length due to a screw dislocation 


H =H,+H,= 


tot 


62/06 anG R 
= — (| — — —)log— bt iets ed 
2G a(> z) mire a 
due to the dilatational strains of an edge dislocation 
62 1 /1—2v\2 (dK R 
= —-({—— } {—— — 1 J log_, geen Orcs) 
a= Fea T=r) (ap) 85, 6 
and due to the shearing strains of an edge dislocation 
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Total volume expansion Sv and total line energy #,,, of an edge dislocation 
are related by 


Lf l—v—20? 1 (dK l-v+v?2 (dG Z| ' 
an —__—__—_(__ — —} |#,,,. (10 
v= 5 pay zz 1) + ees (a ee) 


eee EEA 
+ For the sake of simplicity we have neglected boundary terms in eqn. (5) 
which for R/r)~ 104 are an order of magnitude smaller than those retained, 
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Table 1 
d log K/dp 
Author or d log G d log Cag dlog $(Cr1—C12) 
Metal and technique |@log 3(¢u+2¢)} dp dp dp 
dp 
Cu Birch (1937) 
polyery- | torsional = 2-71 = a 
stalline | vibration 
Bridgman (1935) 3.65 ae ks om 
static si 2 
Cu Lazarus (1949) 
single supersonic 2°84 -= 1-13 2-45 
crystal | pulses 
Ni Birch (1937) 
polyery- | torsional — 1-82 — — 
stalline | vibration. 
Bridgman (1935) 2.88 pbs = en 
static 
Bridgman (1929) 18 a = 
static : 
Al Birch (1937) 
polycry- | torsional — 7-6] = al 
stalline | vibration 
Bridgman (1935) 5.83 = — ae 
static mi 
Hughes and 
Maurette (1956) 5.9 9.8 J #) 
supersonic i a 
pulses 
Al Lazarus (1949) 
single supersonic 5:47 = 7-49 — 
crystal | pulses 
Fe peas and 
ate aurette (1956) bes 
ees supersonic Paes oie = < 
pulses 


NaCl Lazarus (1949) 
single supersonic 24-6 — 2-1 25:3 
crystal | pulses 


Experimental values for the pressure dependence of the elastic constants, 
in units of L0~® em?/kg=10~* bars-!. Bridgman’s (1935) data are corrected 
for the redetermination of the compressibility of steel (Bridgman 1940) as given 
by Birch (1942). Misprints in table V of Lazarus (1949) have been corrected. 
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In table | we have collected experimental data of the pressure dependence 
of the elastic constants of a number of f.c.c. metals, of «-Fe and of NaClf. 
Lazarus’ (1949) results on Cu single crystals are seen to differ appreciably 
from the polycrystalline data obtained by Birch (1937) and Bridgman 
(1935). Since we were unable to decide between the two sets of data we 
carried through the evaluation of eqn. (7) to eqn. (9) for both of them. 


co) kg/ yp we dG [screw Tyodge Tcdge | Tedge 
mm2 | mm? dp dp ; 
Lazarus 4:38 | 0-80 | 0-056 | 0-034 0-07 0-10 
ee 15440 | 4590 | 0-36 |--— es fae aA 
Birch, 5:64 | 1:24 | 0-106 | 0-047] 0-13 | 0-18 
Bridgman 
Ni 18000 | 7850 | 0-31] 5-10 | 1-42 | 0-11 | 0-069! 0-12 | 0-19 
Al 6780 | 2600 | 0:33 | 3-69 | 1-94 | 0-17 | 0-039 0-20 0-24 
Fe 16640 | 8060 | 0-29 | 4:23 | 3-40 | 0-30 | 0-060 | 0-32 | 0-38 
NaCl 2340 | 1500 | 0-24 | 5-77 | 2-05 | 0:317 | 0-25 0-30 | 0-55 


The physical constants used in the calculations are given in table 2. 
Asthe results of the numerical evaluation we list the dimensionless quantities 


f 6u b 
~ Qlog (R/r9)” 


Dlog R/r, has the meaning of the volume expansion due to a dislocation 
line of length 6, measured in atomic volumes Q (Q, = 63/21? for f.c.c. crystals ; 
Q = 463/33? for b.c.c. crystals; Q = 063/23? for the NaCl-structure). [sw 
refers to screw-dislocations (calculated from eqn. (7)), [g°* (calculated 
from eqn. (8)) and I’,°*2° (calculated from eqn. (9)) refer to the expansion due 
to the dilatational and shear parts of the strain-field of an edge dislocation, 
respectively. Finally we give the sum [t#e = [°° + De, Comparing 
the contributions of dilatational and shearing strains to the volume expan - 
sion of an edge dislocation the latter is seen to be dominant (particularly 
for the metals) and of the same magnitude as the total volume expansion of a 
screw dislocation. 

The dislocation densities in undeformed crystals correspond to 
log R/rys=8—10. The volume expansions of f.c.c. metals per dislocation 
length 6 are accordingly found to be of the magnitude one to two, in 


(11) 


+ For data on other materials see Birch (1942), Lazarus (1949) and Bridgman 
(1952, p. 386). 


P.M. 2*k 
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agreement with the results of Stehle and Seeger (1956) derived from a 
model based on inter-atomic forces. Stehle and Seeger’s results can be 
approximated by an expression of the form of eqn. (11), namely 


il a it 
screw — a eet eS ; : ; F . 5 12 
2 =) TENE CESS i} (12) 


Here « is the exponent in the repulsion potential between neighbouring 
atoms. A reasonable value for metals is%#=14. This gives MX’ =0-18. 

Lomer’s (1957) relation between volume expansion and stored elastic 
energy, involving Griineisen’s constant y, reads in our notation 


we ae 

O= Wo ona el BNP a Breas Wea 
where W,,,=W,+W, is the average density of the elastic energy. 
Assuming Hookeian forces between nearest neighbours the following 
relation between y and dK /dp can be derived (Leibfried 1955): 


dk 
= —-l. «te nek, oct ee eee ene 
ay iP (14) 
Lomer’s expression, which takes the form 

rb OS a 
ees zz = ') Wcites = Sey 

can now be compared with eqn. (1), which we rewrite as 

1 (dK dlogG dlogKk 

=—(—— -l —= = — P ve ph Le 
rer ) Fra ( dp dp )W. ae 


We see that Lomer’s formula yields either too large or too small values for ©, 
depending on whether dlog K/dp is larger or smaller than d log G/dp. 
From table 1 we see that the first possibility obtains in the case of Cu and 
Ni, whereas the second possibility applies to Aland Fe. In an isotropic 
medium with v=+4 shear strains contribute 7/9 of the line energy of an edge 
dislocation. If we consider an equal number of edge and screw dislocations 
we find that 13/15 of the stored energy is due to shear strains. In the 
application of eqn. (16) to dislocations the term due to W, is therefore in 
general not negligible. 

Wand W, have been defined as the reversible work done per unit volume 
in establishing the strains. The stored energy measured in a calorimetric 
experiment differs from the W’s by the heat absorbed from the temperature 
bath if the deformation is done isothermally. The heat Q absorbed per 
unit volume is proportional to the average dilatation © and given by 
(Darken and Gurry 1953, p. 183) 


Q=TekK @ . 2 Soe eee nei 


} Equation (14) is obtained by combining Leibfried’s equations (83.6) and 
(83.11). The latter equation has been corrected according to a private com- 
munication by G. Leibfried and W, Ludwig. ; 
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where 7’ denotes the absolute temperature, « the coefficient of thermal 
volume expansion and K” the isothermal bulk modulus. Using the 
_ thermodynamic equation 

KS Tx 

en =! of KS 
K C 
(KS =adiabatic bulk modulus, Cc, =specific heat at constant pressure per 
unit volume) and Griineisen’s relation 


(18) 


Dp 


y= Kga/c,, . ‘. ; A : i a (19) 
we obtain from eqn. (17) 
" KS—KT K'6 
gel <a (20) 


By combining eqn. (1) and eqn. (14) it is found that K"©/y is of the order 
of magnitude of W,,,. Since (KS—K")/K” is of the order of 10- (e.g. 
3 x 10? for Cu at room temperature, see Overton and Gaffney 1955) we see 
that W,,, can be taken equal to the stored elastic energy, as originally 
asserted by Zener (1942). A comparison between the theory and measure- 
ments of volume expansion and stored energy due to cold-work has been 
given elsewhere (Seeger 1957b). 
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ABSTRACT 


Experiments have been made on the small angle scattering of x-rays from 
polycrystalline copper and single grains of aluminium deformed by cold work 
or fatigued by pulsating tension. The results confirm in a direct and con- 
clusive way that almost all of the scattering is due to double Bragg reflections 
as proposed by Beeman et al. No surface scattering has been observed. 
Double Bragg reflections have been eliminated in small angle scattering 
experiments using cold neutrons: the results indicate the magnitude of the 
scattering from defects in the volume of the metal. The magnitude of the 
effect is consistent with scattering from dislocations. There was no observ- 
able effect from fatigued metals. Hence there is no direct evidence for the 
existence of cavities in metals as a result of small angle scattering experi- 
ments with either x-rays or neutrons. 


§ 1. INTRODUCTION 


The density changes associated with lattice imperfections which extend 
over several atomic dimensions can give rise to small angle scattering of 
x-rays or thermal neutrons. Hence, it should be possible to obtain 
information from small angle scattering experiments about dislocations 
and also about small voids which may be introduced into crystals by cold 
work, fatigue, quenching or radiation damage. 

Quantitative small angle scattering experiments on pure metals using 
x-rays (of wavelength 1-54 A) were first carried out by Blin and Guinier 
(1951, 1953) and Blin (1954, 1955) who examined the scattering from 
polycrystalline aluminium, nickel, copper and zinc in the range of scattering 
angles 55’ to 7°, i.e. in the range of s from 0-010 to 0-079 A+, where 
s=2sin 0/A~ 26/A; 20 is the scattering angle, and Ais the wavelength. They 
found substantial scattering from metals which had been plastically 
deformed in tension or by cold rolling. Blin and Guinier considered the 
observed intensities too large to be due to dislocations (Dexter 1953) 
and interpreted their results in terms of a concentration of 1015 to 1017 
per cm® of cavities of a size corresponding to about 300 to 7000 atomic 
volumes ; each of these cavities might have been formed by the aggregation 
of vacancies produced during deformation. However, Brammer (1955) 
found an immediate increase in scattering from a nickel foil stretched and 
examined at 77°K. Since vacancies are not thought to be mobile at this 
temperature, it is difficult to see how voids could be formed by clustering 
of vacancies. 
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More recently Blin (1957) has suggested that his results can be explained 
as being due to scattering from density changes associated with the first 
order elastic strains around edge dislocations. Seeger (1957 a, b) also 
believes that the observed scattering is due to dislocations, but considers 
that second order elastic effects predominate. Both these theories have 
been discussed recently by Atkinson and Hirsch (1958) who showed that 
the previous theoretical calculations either contained errors or were 
incomplete. and that the theoretically expected scattering from a reason- 
able arrangement and density of dislocations is in fact very much less 
- than that observed in x-ray experiments. 

Quite apart from these theoretical considerations it is difficult to inter- 
pret certain features of the x-ray scattering in terms of defects introduced 
during cold work. Thus, for example, the original results of Blin (1954) 
and recent more detailed studies by Smallman (private communication) 
and Wild ef al. (private communication) show that the scattering (at a 
fixed angle) at first increases with increasing degree of cold rolling, but 
that a maximum is reached after 10 to 20% reduction in thickness, beyond 
which the intensity decreases steeply and continuously with increasing 
deformation up to 95%, the maximum reduction used. In contrast to 
this variation of the intensity of scattering it is expected that the density 
of defects, and therefore the intensity scattered, should increase with 
increasing deformation. 

Recent experiments by Beeman ef al. (1957) have shown conclusively 
that a large proportion of the scattering from cold-worked metals can be 
explained in terms of a double Bragg reflection mechanism (Neynaber 1955). 
Thus, for example, Webb (1956) found that x-rays (FeK, radiation) 
scattered at 14° are partially polarized with their strong electric vectors 
in a direction radial to the undeviated beam. This polarization can only 
be explained by two successive large-angle Bragg reflections. The 
direction of the electric vector suggests that the process is one in which a 
ray reflected from one region of a grain suffers a further reflection by another 
region in which the reflecting planes are somewhat misaligned by a simple 
rotation about the direction of the once-reflected beam. ‘The particle 
size- and strain-broadening were considered to be negligible. However, 
the polarization decreases as the degree of deformation is increased. This 
can be explained by taking into account the particle size- and strain- 
broadening of the Bragg reflections which become important at the larger 
deformations. Defects associated with the surface have also been shown 
to cause scattering at small angles (Robinson and Smoluchowski 1956, 
Zhurkov and Slutsker 1957). 

On the basis of these experiments it appears therefore that the major 
part of the small angle x-ray scattering from metals deformed by cold 
rolling or in tension is due to double Bragg reflections and possibly also 
partly to the surface effect. 

More recently, small angle scattering has also beeu measured from 
metals which have been deformed by fatigue. Franks and Holden (1955) 
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have observed the effect of about 3-5 x 10° cycles of pulsating tension 
(between 0 and 9-5 kg mm~*) on the scattering from copper foil of thick- 
ness 60. An individual grain was examined with CuK, radiation, 
and x-ray photographs showed streaks through the central spot extending 
out to about 2°. These streaks were interpreted as being due to disc 
shaped cavities, about 1004 in diameter, lying on (111) planes. No 
measurement was made of the absolute intensity of the effect. 

Since 1955 the present author has studied the small angle scattering 
from deformed metals with the following aims in view. Firstly, it was 
thought desirable to demonstrate the double Bragg reflection mechanism 
by means of more direct experiments, and to determine whether this effect 
could also account for the scattering from metals deformed in fatigue. 
Secondly, a series of experiments were undertaken to study the small angle 
scattering with cold neutrons. By using neutrons of sufficiently long 
wavelength the possibility of double Bragg reflections was eliminated, and 
since the specimens examined were very large, surface effects were negligible. 
In this way it was possible to measure the true small angle scattering due 
to density changes associated with defects. A preliminary account of 
the neutron experiments has been published elsewhere (Atkinson and 
Lowde 1957). 

In the present paper the x-ray experiments will be described and some 
new neutron scattering results discussed. The neutron scattering 
experiments will be published in more detail elsewhere. 


§ 2. THe X-Ray EXPERIMENTS 


2.1. Equipment 


A point focusing monochromator and camera have been used in this 
work. Only a brief description of the instrument will be given here as 
full details are to be published elsewhere. A schematic diagram of the 
equipment is shown in fig. 1. Radiation from the copper target of the 
x-ray tube (operated at 20 ma and 30 kv) passes through a slit S, onto a 
single crystal of lithium fluoride plastically bent to conform to a toroidal 
surface so that the CuK , radiation is reflected from the (200) planes and 
brought to a point focus 23 em away. The specimen, placed midway 
between the crystal and the focus, is irradiated over a region of about 
1'8mmx7mm. Ifmmediately behind the specimen is a second slit S, 
to limit the effect of parasitic radiation from the LiF crystal. This slit is 
carefully adjusted so that it just does not touch the direct beam. With the 
slit settings used in this experiment the parasitic scattering is reduced 
in one direction ; the area indicated by the dotted lines in fig. land by the 
elongated region about the direct beam in figs. 2, 9, 11 and 14, Plates 14, 16, 
17 and 18 respectively, is obscured by parasitic radiation. This limits 
the smallest angle observable to 14°. The whole instrument is evacuated 
to eliminate air scattering. 
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The small angle scattering is recorded on a photographic film placed in 

a plane perpendicular to the direct beam. A round beam stop, of diameter 
6mm, intercepts the direct beam; this can be moved momentarily to 


obtain a picture of the direct beam. Measurements of the absolute 
intensity were made with a Geiger counter with a 1-2 mm square aperture. 


Fig. | 
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Schematic diagram of the x-ray point focusing monochromator and camera. 


The direct beam intensity, about 107 counts per second, was measured 
with a counter, using nickel foil absorbers; a correction was made for 
radiation of higher orders reflected from the LiF crystal. 


2.2. Deformation by Cold Work 


2.2.1. Cold rolled copper 


A photograph of the small angle scattering from heavily cold-rolled 
copper foil, 30 » thick, is shown in fig. 2, Pl. 14 (exposure 70 hours). The 
material used was of commercial purity. (Material from the same batch 
will henceforth be referred to as ‘commercial copper ’.) It is seen that 
streaks are produced, some of which do not pass through the trace of the 
direct beam. It will now be shown that all these streaks are due to double 
Bragg scattering. 

As the x-ray beam passes through the polycrystalline metal, the rays 
which are Bragg reflected must lie on cones centred about the direct 
beam (e.g. Coo) in fig. 3). Consider a particular singly reflected ray, 
e.g. OA: in passing through the rest of the specimen this ray itself may 
again be reflected (as described in § 1). The twice reflected ray will also 
lie on one of a series of cones (e.g. Cag 299 OF Coo, 111 IM fig. 3, correspond- 
ing to successive reflections from planes of the same and different 
spacings, respectively) only one of which, Cy99, 999 Will intersect the direct 
beam in the forward direction. If the crystal consists of many randomly 
orientated grains and if there is no strain broadening, the resultant scatter- 
ing will have circular symmetry about the direct beam, and will decrease 
steadily as 1/20 with increasing scattering angle 26. However, for angles 
20 > |20,4 — 28 ny|, & new term is added to the intensity, for it is then 
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possible for a single reflection described by the Bragg angle 6,,, to be 
followed by a second reflection at a different set of planes (,,, ny) and vice 
versa. Also, an abrupt increase in the intensity at small angles will occur at 
20 = (27—40,,,,), if there is a Bragg angle 6,,,,, near to 37. For copper and 
aluminium with CuK, radiation, the first increase should occur at 5° and 
34° respectively, given by 20= (269 —263;,;). On the other hand, the 
copper foil used in this experiment showed strong preferred orientation, 
and the small angle pattern expected from the double Bragg scattering 


Fig. 3 


Q INCIDENT BEAM 


(200) 
0 *\ CONES of TWICE 
SS ai) 
REFIEG TE LY CRATS 
\ YY C200, 111 ‘e 
200,200 
aot \ \ 
yk A PARTICULAR 
SS SINGLY 
C200 Ne a 
CONE OF : 
ONCE REFLECTED > —~ ‘ Der Ee Pepe ey 
(200) RAYS 
(200) , (111) 
SMALL ANGLE BEAMS 
(200) (200) 
\\ 


TRANSMITTED BEAM 


effect no longer has circular symmetry. The preferred orientation of the 
grains in the specimen is evident from figs. 4 and 5, Pl. 15. Figure 4 shows 
a transmission photograph taken in the point focusing instrument by 
placing a film parallel to the specimen and about 2em from it. This 
shows some of the first 111 and 200 reflections occurring under conditions 
the same as those in which fig. 2, Pl. 14 wastaken. (The central blackening 
. . . . . © 
is caused by the direct beam which is out of focus at this point. The 
size of the direct beam also causes some of the spreading of the reflections. ) 
Most of the other Bragg reflections occurring are shown in fig. 5; PEs 
In this case the film was bent in the form of a half cylinder. The x-ray 
beam (from a copper target with a nickel filter) passed along the axis of 
the cylinder onto a small piece of the copper foil. The normal to the 
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plane of the foil was parallel to the beam. The plane containing the 
X-ray beam and the rolling direction of the foil cuts the film along a line 
_ marked 0°; when this plane is rotated 90° about the beam it then intersects 
the film along the lines marked + 90°, ? 

A qualitative examination of figs. 4 and 5, Pl. 15, suggests that the small 
angle streaks which pass through the trace of the direct beam in fig. 2, 
Pl. 14, are associated with strong large angle Bragg reflections. To explain 
the details of the pattern, pole figures of the normals to the (111) and (200) 
planes have been determined photographically, and are shown in figs. 6 


Pole figure of the preferred orientations of the normals to 111 planes for cold 
rolled commercial copper. The plane of the foil is parallel to the plane 
of the paper. ‘ R.D.’ indicates the rolling direction. 


and 7. Using these pole figures, it is possible to predict the nature of the 
small angle scattering pattern expected ; this predicted pattern is shown 
in fig. 8, Pl. 14. The dotted lines on fig. 8 indicate the positions of other 
double reflections deduced from fig. 5, Pl. 15, but for which pole figures 
have not been determined. The streaks which do not pass through the 
origin are due to reflections first at (200) planes and then at (111) planes: 
the pole figures show that reflections in the reverse order, 111 and then 200, 
occur in a range of angles not covered by the photograph. Comparison 
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of fig. 2 with fig. 8, Pl. 14 shows that all the observed streaks are explained 
completely by the double Bragg reflection mechanism. 

The intensities of the scattering at some points of the pattern (tig. 2) 
have been measured with the Geiger counter; at @ the intensity is 160 
electrons/atom and at 6, 85 electrons/atom, i.e. of the same order of mag- 
nitude as found previously by Blin. The intensity of the scattering was 


Pole figure of the preferred orientations of the normals to 100 planes for cold 
rolled commercial copper. The plane of the foil is parallel to the plane 
of the paper. ‘R.D.’ indicates the rolling direction. 


unaffected by removing a surface layer 2:5 « thick from each side of the 
specimen by electropolishing. Removal of a further 2-5 4 from each 


side reduced the intensity slightly, but not by more than would be expected 
from the consequent reduction of the quantity of material in the beam. 


2.2.2. Aluminium deformed in tension 


Figure 9, Pl. 16 shows the scattering (64 hour exposure) from a large 
single grain of Al which has been extended ~ 18% by atensile load. After 
deformation, the large grain concerned was reduced to a thickness of 100 be 
by removing 100 ~ from each side by electropolishing. Using a Geiger 
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counter, the grain was orientated so that it would give a strong 200 reflection 
for CuK, radiation when placed in the point focusing instrument. Figure 
10, Pl. 16 shows that in fact two strong single 200 reflections occurred ; 
the spot a on this photograph produces the intense streak marked a in fig. 9: 
the second 200 spot (6) probably produces streak 6 (fig. 9). Three other 
strong reflections have been observed (outside the range of fig. 10): two 
are of the type 311, producing streaks c and d (fig. 9); the third (of type 
331) may also contribute to b and c. Spot e is due to the 311 reflection 
giving rise to streak d, followed by a 220 reflection. 


2.3. Deformation by Fatigue 

A small fatigue machine was constructed to apply to a thin foil a tensile 
load which alternated between a value of about zero and a fixed value, 
at a rate of 1200 cycles per minute (‘ pulsating tension’). The machine 
was automatically switched off when the specimen stretched more than a 
certain amount. 

In these experiments the full load was applied immediately and in each 
case the previously recrystallized specimen would stretch considerably 
during the first few cycles and then more slowly during the rest of 
its life. 


2.3.1. Copper 


Figure 11, Pl. 17 shows the small angle scattering (45 hour exposure), 
from a piece of pure copper foil (supplied by Johnson, Matthey & Co. Ltd.) 
30 » thick, which had broken after ~5 x 10° cycles of pulsating tension 
between 0 and about 7kgmm~*. The specimen stretched a total of at 
least 10°, during the fatiguing. The corresponding 111 and 200 reflections 
are evident in fig. 12, Pl. 17. Thisshows that the numerous grains in the 
specimen have a wide range of orientations. Hence the small angle 
pattern has streaks occurring at many different angles. Streaks marked 6 
in fig. 11 are due to double reflections of the type 200, 111. The intense 
scattering produced (e.g. 400 electrons/atom at a) shows that such a wide 
array of orientations increases the probability of the double scattering 
process. 

The intensity of scattering at a is approximately the same as that 
resulting from a foil of the same material which had been stretched by 
the same amount (about 10%) but not fatigued. Removal by electro- 
polishing of 5 from each side of either of these foils resulted in a negligible 
change in scattering. 

A piece of the commercial foil was fatigued under similar conditions 
to the above. It produced no observable small angle scattering. This 
can be understood by referring to fig. 13, Pl. 18, taken under conditions 
identical with those corresponding to fig. 12. It is seen that the few 11] 
or 200 reflections which occur are not sufficiently wide to produce scattering 
observable at angles greater than 13°. 
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2.3.2. Aluminium 

A piece of aluminium foil 250 . thick, with very large grain size (~ 1 em) 
was fatigued; during this process a selected grain stretched ~ 14%). 
The specimen did not break. This grain was then orientated to produce a 
strong 111 reflection. ‘The small angle scattering is shown in fig. 14, PI. 18 
(22 hour exposure): the streak is produced by the strong 111 reflection 
shown in fig. 15, Pl. 19. 

A photograph of the same grain orientated so that no strong III or 200 
reflections occurred showed very little small angle scattering. Because of 
the heavy distortion of the lattice it is difficult to find an orientation for 
which there is no scattering at all. 

The removal by electropolishing of 50 from each side of the specimen 
diminished only slightly the intensity of the streak in fig. 14. 


2.4. Discussion of X-Ray Scattering Results 


The experiments described above leave little doubt that nearly all 
the small angle x-ray scattering from metals deformed in tension, by cold 
rolling, and by fatigue stresses is due to double Bragg reflections. This 
mechanism, first suggested and investigated by Beeman, Brammer, 
Neynaber and Webb, has therefore been confirmed for several different 
types of deformation. It follows from the neutron scattering results (§ 3) 
that the x-ray scattering due to defects contributes not more than about 
1% to the total scattered intensity. It is therefore doubtful whether this 
scattering can ever be detected with x-rays. 

The present results also indicate that surface scattering is relatively 
unimportant. 

In view of these results it should be noted that there is now no evidence, 
from X-ray scattering experiments, for the possible existence of cavities 
or cracks as a result of cold work or fatigue. 


§ 3. THe Neutron Scatrerinc EXPERIMENTS 


In order to prevent entirely the occurrence of double Bragg scattering 
in a metal specimen, it is necessary to use radiation of wavelength greater 
than the Bragg cut-off wavelength, d,. for the metal used (A,=2dmax., 
where dmax is the maximum interplanar spacing: for copper A, =4174 
and for aluminium = 4-67 A). Long wavelength x-rays are so heavily 
absorbed that specimens only a few microns thick would have to be 
examined. However, for long wavelength neutrons, thick (~1 em) 
specimens can be used. 

In collaboration with R. D. Lowde at the A.E.R.E., Harwell, experiments 
have been carried out recently to measure the small angle scattering of 
cold neutrons from aluminium and copper deformed by tension, cold rolling 
and fatigue. The results of these experiments are now reported together 
with brief description of the experimental arrangement. A full deseription 
of the experiment will be given elsewhere. 
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The cold neutrons used came from the liquid hydrogen moderator 
(Butterworth et al. 1957) in the Harwell pile B.E.P.O. Filters of Be and Bi 
were placed between the source and the specimen to remove all the neutrons 
of wavelength A<6-5 A (Kgelstaff and Pease 1954). For A> 6-5 A, the 
intensity of the beam incident on the specimen is approximately 
proportional to A~>. | ; 

The flux of neutrons available even from this very intense source is 
still several orders of magnitude smaller than the x-ray flux used in the 
experiments described in § 2. Hence it is necessary to relax considerably 
the resolution of the apparatus so that, at the setting of the ‘ small-angle a: 
counter used, neutrons could be detected over a range of scattering angles 
between 3° and 11°. 


Fig. 16 


| Sse 2:54cm. 


NEUTRON 
BEAM 


GORPIER ALUMINIUM 


Fatigue test-pieces used in the neutron experiment. To the left of each test- 
piece is shown a finished specimen. 


The aluminium specimens were prepared from a length of rod of 99-99%, 
purity from the Northern Aluminium Company. The copper specimens 
were from a batch of metal of 99-995°% purity supplied by I.C.I. (Metals 
Division) Birmingham. This material was cut into suitably shaped test- 
pieces (the fatigue test-pieces are shown in fig. 16), and then recrystal - 
lized. One piece of each metal was kept as a standard ‘ recrystallized ’ 
specimen. The tensile and fatigue test-pieces were deformed using 
a 10 ton Losenhausen fatigue machine. Each specimen was then 
carefully machined from the gauge length of the test-piece. Finally, the 
copper specimens were electropolished ; 50 were removed from each 
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side of the tensile specimen (the only copper specimen which had been 
machined on surfaces normal to the neutron beam) and 12 from each 
side of the fatigue specimen. 

The final forms of the aluminium and copper specimens are shown in 
fig. 16. Details of the deformations are shown in table 1. 

‘In the actual experiment, a deformed specimen was compared with a 
recrystallized but otherwise identical one; the difference in the scattering 
from the deformed and recrystallized specimens should be due solely to 
defects produced in the volume of the metal by the deformation. The 
magnitude of this ‘ defect ’ scattering, expressed as a differential scattering 


= = 


Table 1. Small Angle Scattering of Neutrons from Deformed Aluminium 
and Copper 


Nuclear differential cross 


Metal State Details section for scattering by 
defects (millibarns/steradian/atom 


Recrystallized | at 550° c for 1 hour. 


Extended about 50%. 0-03 + 0-08 
Aluminium | Fatigued in push—pull to fracture after —0-05+0-08 
3-17 x 10° eycles at 
+ 2-3 kg/mm? 


Reerystallized at 700°c for 4 hours in a _ 
vacuum of about 10-° mm Hg. 


Extended about 50%. 1-5+0°8 
Cold rolled thickness reduced by 90%. 1:9+0°8 
Copper Fatigued in push-pull. The specimen —0-3+0-8 


fractured after being subjected 
successively to the following 
treatments : 
17 x 10° cycles at +4:7 kg/mm? 
1-2 x 10 eyeles at +5:-8 kg/mm? 
1-0 x 10° eyeles at +6°9 kg/mm? 
4-2 x 10° cycles at +7:6 kg/mm? 


cross section for neutrons in units of millibarns steradian—! atom—!*+ is 
shown in the last column of table 1. The standard deviation is also shown. 
The absolute magnitude of the results may be in error by +50% : however 
this does not effect the relative values of the quantities or the relative 
statistical errors. 

It is seen that the cross section for small angle scattering by cold-worked 
or fatigued aluminium is probably less than 0-08 mb str atm-!. A 
small effect is observed for cold-worked copper. The statistical probabi- 
lity that this effect is real, and is not due to chance, is 92° for the extended 
and 98% for the cold rolled specimens. The fatigued copper gives no 
observable scattering. 


+ A table relating the scattering expressed in neutron or x-ray units to 
[j1,, where I, is the scattering per cm? and J, is the intensity scattered from 
one atom under the same conditions, is given in Atkinson and Hirsch (1958), 
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For comparison, the x-ray results have been converted to average 
differential cross sections for neutrons, for the range of wavelengths and 
angles used in the neutron experiment. Table 2 shows average cross 
sections calculated in this way from the results of Blin (1954) and Smallman 
(private communication). The present x-ray experiments confirm these 
values, 

It is seen that not more than about 1% of the scattering observed in 
the x-ray experiments is due to true small angle scattering from defects in 
the volume of the metal. 


Table 2. Some X-Ray Scattering Results Converted to Equivalent 
Differential Cross Sections for Neutron Scattering 


Equivalent cross section 
Specimen Observer for neutron scattering 
(millibarns/steradian/atom) 


Copper 
Cold rolled 20% Smallman > 150 
90% Smallman SH) 
95°, Blin > 35 
Extended 5% Blin > 350 
Aluminium 
Extended 5% Blin Ss 70) 


As for the magnitude of ‘true’ scattering, the neutron experimental 
sensitivity seems to be just sufficient to detect an effect for copper. Al- 
though no scattering was observed from aluminium, this experiment 
enables an upper limit to be set for the effect. To determine whether 
these results are consistent with scattering from dislocations, average 
neutron scattering cross sections have been calculated using the theoretical 
treatment of Atkinson and Hirsch (1958), assuming reasonable dislocation 
densities of 10! and 101° em~? for heavily cold-worked copper and alumi- 
nium, respectively, and supposing that half the dislocations are in the 
edge orientation. The values obtained, 2 mb str~! atm™! for copper and 
0-04 mb str-! atm-! for aluminium, are in order of magnitude agreement 
with the experimental neutron scattering cross sections of table 1. 


§ 4. CONCLUSIONS 


The x-ray experiments show that most of the scattering at small angles 
is due to double Bragg reflections. No surface scattering effect has been 
detected. 

The neutron scattering experiments show that the true scattering from 
defects is very small, and that the scattering detected for rolled and ex- 
tended copper can be explained satisfactorily in terms of scattering from 
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a density of about 10! dislocations cm~*. ‘Thus there is no evidence for 
the existence of cavities in cold-worked or fatigued specimens of copper 


or aluminium. 
In view of the small neutron scattering intensities, and considering the 


large effect of the double Bragg scattering mechanism for x-rays, it seems 
doubtful whether the small angle scattering method can be usefully applied 
at present to the study of dislocations. 
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ABSTRACT 


A single crystal of nickel fluosilicate has been used for the low temperature 
alignment of cobalt-58. Angular distribution measurements of the 805 kev 
y-radiation are consistent with the 8-transition being pure Gamow-—Teller, 
the Fermi admixture being given by 


Cy*|My|2/{Cy2|My|?+ Car?|Mar|2}= — 0-003 + 0-005. 


§ 1. IyrRopucTION 


The recent experiments relating to the violation of parity conservation 
in B-decay (Wu ef al. 1957, Ambler, Hayward, Hoppes, Hudson and Wu 
1957) have emphasized the importance of obtaining complete information 
on the character of 8-transitions so as to provide a test of B-decay theory. 
In particular the presence of small admixtures of Fermi interaction in 
predominantly Gamow-—Teller transitions can lead to interference terms 
which may radically alter the spatial asymmetry of B-rays in the ‘ parity ’ 
type of experiment. At the same time, the presence of such admixtures 
can also be detected in experiments which do not rely on the breakdown 
of parity conservation, such as the study of the anisotropy of y-radiation 
from aligned nuclei in which the y-emission is preceded by a f-transition. 
58Co which might exhibit such an admixture is so far the only isotope that 
has been studied by both methods and in this case the results of the two 
are in disagreement. 

The decay scheme of *8Co is given in fig. 1 (see, for instance, Frauenfelder 
et al. 1956) and shows that most of the disintegrations proceed via 
K-capture or positron emission directly to the 805 kev level in Fe. The 
study of the weak branch through the 1-62 Mev state showed this state to 
be of spin 2 and the 815 kev y-ray to be mixed with 6(#2/J11)= 2-2 + 0:3. 
The spin assignment to the other states is made on the basis of nuclear 
alignment experiments by Daniels et al. (1952). Recently the ground 
state spin of 58Co has been confirmed by a paramagnetic resonance 
experiment by Dobrov and Jeffries (1957). Daniels et al. also showed 
the B-transition to be principally Gamow-Teller. A statistically more 
accurate result was obtained by Griffing and Wheatley (1956) who used 
the same technique. They showed that best agreement between the 
theoretical and observed y-ray angular distributions is obtained for 

aa Ox? | My [2 
~ Ogr"|Merl + Cx*|MeP 
+ Communicated by the Authors. 


=O-11 + 0-4 
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Cy», Cap ave the Fermi and Gamow-Teller coupling constants and My, 
May the corresponding nuclear matrix elements for the f-transition.. 


A subsequent analysis of the measurements of Daniels et al. showed that a 


value of \ of 0:1+0-1 would be consistent with their results: this figure. 


is referred to in the paper of Griffing and Wheatley. 
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The decay scheme of cobalt-58. 


In the ‘parity’ experiment of Ambler et al. the backward-forward 
asymmetry coefficient « was found to be + 0-325 + 0-047 for the positrons 
from polarized **Co. This value is in agreement with that predicted 
by the two state neutrino theory for a pure GT transition (+4) and is 
incompatible with a value of A greater than 0-03; this conclusion is 


supported by a similar experiment of Postma et al. (1957). Confirmation of 


this ‘parity ’ result is provided by the B-y circular polarization correlation 
measurements of Boehm and Wapstra (1957) who showed J to be less than 
0-06. Ambler, Hayward, Hoppes, and Hudson (1957) pointed out that 
the apparent, discrepancy between the results of the y-ray and the B-ray 
experiments might arise from the F-GT interference terms being largely 
imaginary but this hypothesis they showed to be untenable in the light of 


a fP-y correlation experiment on polarized ®8Co nuclei in which they 
found no evidence for such terms, 
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Since the establishment of a disagreement between these two types of 
experiment could be important in the development of future B-decay theory, 
it was decided to make a further study of \ by the y-ray method under 
experimental conditions which avoided some of the defects of the earlier 
alignment experiments. 


§2. DEsIa@N or EXPERIMENT 


The principle of the method for studying A through the y-ray angular 
distribution is as follows. The y-ray intensity at an angle @ with the axis 
of alignment is expressible in the form (see, for instance, Gray and Satchler 
1955) 


W(@)=1+ B,U,F,P,(cos 6)+B,U,F,P,(cos 6) 


B, is a temperature dependent term describing the degree of alignment, 
and F, an angular correlation coefficient for the y-ray transition. The 
character of the §-transition is contained in the terms U, which can be 
written as 

U,=dAU F+(1-A)U 


U,*, U,S™ being the values appropriate for Fermi and Gamow-Teller 
transitions. For Co all these terms are known except A and B,. Since 
in the simple case of completely anisotropic hyperfine structure (hfs) B,, B, 
each depend only on the Boltzmann parameter f(=A/2k7' where A is 
the hfs constant ; see below), it is sufficient to measure W(6) for two values 
of @ and to compare these with the theoretical values in order to determine A. 
It is usually convenient to choose to measure at the angles 6=0 and 47 
where both P,(cos @) and P,(cos @) are changing least rapidly. This was 
the method used in earlier alignment experiments, but these suffered from 
two defects: 


(a) Since the experiments were carried out on a number of crystals 
arranged with their axes parallel thermal equilibrium amongst the con- 


stituents was not ensured and therefore the B, must be replaced by B, 


where B, is some average over all the crystals. Such a change in the 
calculated B, could lead to erroneous values for A. 


(6) The Tutton salt crystals used in the experiments have two alignment 
axes inclined at about 70° to each other, neither of which is coincident with 
a crystallographic axis. For °*Co it is the axial direction in which the 
y-ray intensity is most sensitive to the value of \ and therefore, since in the 
biaxial Tutton salt this direction is not accessible for study, the measure- 
ments with these crystals are less sensitive to A than those on uniaxial 
crystals. In addition, the error introduced by uncertainty in the exact 
location of these alignment (tetragonal) axes makes the Tutton salt an 
unsatisfactory choice for this kind of measurement. A more favourable 
erystal would be one of the fluosilicate group which have a single axis of 
alignment coincident with one of the crystal axes. 


2L2 
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The fluosilicates M++SiF, . 6H,O form an isomorphous series in which M 
can be any divalent ion of the iron transition group. Experiments by 
Hill et al. (1958) using a magnetically diluted nickel salt of this series have 
shown that this can be cooled to temperatures of the order of 0-02°K 
by adiabatic demagnetization using H/T values of 25 kilo oersted/degree. 


W (%) 
ime) 1:05 HO 15 120 


w(O) 


) = O-t 


The normalized counting rate along the axis W(0) is plotted against that in the 
equatorial plane W($7) for the 805 kev y-ray from aligned ®8Co. The 
value 1 corresponds to the isotropic distribution. The solid curves show 
the expected dependence for different values of the £-transition mixing 
parameter - 


A =Cy"| Mpl?/Oy*| My? + Cgr*| Merl’). 
Only portions of the A=0-1 and 1 curves are shown, 


Bleaney and Ingram (1951) have studied the paramagnetic resonance 


spectrum for Co** in this lattice and have shown it to be describable by a 
spin Hamiltonian of the form 


H = ASI, + B(S,1,+8,1,)+98H-S 


where S=} and A=0-0184+0-0004 cm—, B=0-0047 + 0-0002 em-! for 
stable cobalt: the sign of A is taken to be positive (Abragam and Pryce 


——— 
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1951). Thus at temperatures of the order of 0:01°K considerable align - 
ment of the Co nuclei should occur in zero external magnetic field through 
the magnetic hfs (Bleaney 1951). The existence of the B term leads to a 
deviation from the uniform level spacing 14 and also to a mixing of the 
nuclear states. Both of these depend on the ratio B/A, a factor which is 
independent of the particular cobalt isotope. Thus for a given B/A the 
population functions B, can be evaluated in terms of 8 only. With 
these values a family of curves can be drawn showing the dependence 
of W(0) on W(37), each member of the family corresponding to a single 
value of A and each point on the curve corresponding to a different value 
of 8. Anumber of these curves are shown in fig. 2 and indicate the substan - 
tial differentiation between small values of A which can be obtained by this 
method. The sensitivity to A comes from the great difference in the 
y-ray intensities along the axis for the two cases: the minimum value 
of W(0) for the Fermi transition is 0-06 whereas for the Gamow-—Teller 
transition it is 0-822. It should be pointed out that the derivation of i 
from a measurement of W(0) and W(47) does not rely on any knowledge 
of the absolute temperature of the crystal but only on the validity of the 
spin Hamiltonian, the value of B/A and the sign of A. 


§ 3. EXPERIMENTAL 


Cobalt-58 was produced by the *®Ni(n, p)>®Co reaction using the fast 
neutron flux of the A.E.R.E. pile. Carrier free cobalt was then separated 
chemically from the nickel. The cobalt was largely °8Co, only about 
4% being due to Co. This low yield is in accord with that predicted by 
Hughes (1954). A single crystal of (10°% Ni, 90% Zn) Sif, .6H,O weighing 
3-1 g and containing about 70 pcuries of Co was grown from active 
solution. The surfaces of the crystal were clear and its axis could be located 
accurately. This was then held in a small Perspex carriage and mounted. 
with its axis horizontal on a rigid glass spill in the low temperature cryostat. 
Also mounted on the spill was a pill of Mn(NH,),(SO,). . 6H,O which acted 
on cooling both as a guard ring for reducing the heat influx to the crystal 
through its support and as a ‘ getter ’ for mopping up the last traces of the 
exchange gas. After cooling by adiabatic demagnetization from initial 
fields of up to 35 kilo oersted at 0-9°x, scintillation counters fitted with 
2in. cylindrical Nal(T1) phosphors were placed round the cryostat to record 
the counting rates on the axis of the crystal and in the equatorial plane. 
Each counter had an energy resolution of 8% for 660 kev y-radiation and 
subtended a half angle of 9° at the crystal. 

A pair of discriminators was set above and below the 805 kev photopeak. 
In this way account could be taken of the effect of the higher energy 
y-radiation from *8Co and from ®Co contamination on the counting rate 
in the 805 kev channel. A pulse height analysis showed that the whole 
spectrum could be interpreted in terms of these components. These 
counting rates together with the magnetic temperature were followed for 
about 5 minutes and then the crystal was warmed up to about 1°, the 
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temperature of the helium bath. By demagnetizing from different initial 
fields the whole temperature range was covered. 

The counting rates for any given value of the magnetic temperature were 
independent of whether this temperature was reached by direct demag- 
netization or by demagnetizing to a lower temperature and allowing the 
crystal to warm up. This shows that the temperature was sensibly homo- 
geneous within the crystal over the measuring period of 5 minutes. 
Frequent tests throughout the series of runs were carried out which showed 
that no systematic trends, due, for instance, to instability or to fatigue 
of the counters, were influencing these readings. These tests took the 
form of demagnetizations from small fields in which case no cooling was 
produced and the isotropic distribution was obtained. 

The counting rates for the channels on each counter normalized to the 
warm counting rates were then plotted one against the other and are 
shown in fig. 2. Each point is an average of a number of readings (about 
20/point). 


$4. RusuLtrs anp Discusston 


It is seen that the uncorrected experimental points agree closely with the 
A=Ocurve. The results were treated quantitatively in the following way. 
Each experimental point can be treated as a measurement of the value 
of X with a probable error and hence a weight arising from the counting 
statistics. In this way the average value of A is found to be + 0-003 with 
a probable error based on the scatter of these points of +0-004. This 
value is in agreement with that to be expected on the basis of the total 
number of counts (~ 10°) for the coordinates of each point (+ 0-003). 

This value for A now requires to be corrected for the following systematic 
effects : 


(i) The finite aperture of the counters (9°) gives rise to a correction to A 
of — 0-003 (see for instance Feingold and Frankel 1955). 


(ii) The presence of pulses due to the 0-815 and 1-62 Mey y-rays lying 
within the 0-805 Mev channel causes a correction which depends on the A for 
this weak low energy B-transition (A’). If A’ =0 the 1-62 Mev y-ray involves 
no correction while the 0-815 Mev y-ray causes a change in A of — 0-001. 
IfA’ is not zero but about 0-5 as the measurements of Griffing and Wheatley 
suggest then A has to be corrected by —0-004. For the present purposes 
we assume 1’ to be zero. 


(ii) The presence of Co in the source was responsible for about 0-13% 
of the total number of counts in the 0-805 mev channel. From the known 
values of A for Co and *8Co in this crystal the angular distribution of 
this radiation can be predicted and leads to a change in the value of A of 
— 0-001. 

(iv) The total thickness of the low temperature apparatus through 
which the y-radiation had to pass was equivalent to a vertical cylinder of 
brass 0-35 cm thick which caused scattering of the y-radiation. The 


ee = 
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‘discrimination level of the counter was such that it could not detect 
y-radiation which had been singly scattered through angles greater than 
30° with the corresponding degradation in energy. A calculation for a 
spherical scatterer of this thickness shows that scattering leads to a 
reduction in A of 0-001. It is considered that this closely represents the 
‘experimental conditions and that a full numerical integration would not 
lead to a significantly different result. This has therefore been used in 
‘determining the best value for X. 
In addition to these there are the following sources of uncertainty : 


(a) Possible misalignment of the crystal by 2° leads to an error in 
aot Fh. 
—0-002" 
(6) The value of B/A is only known to 5% and this gives rise to an error 


of +0-0015. 
Combining these corrections and errors we arrive at the result : 
A= — 0-003 + 0-005. 


This value is consistent with the complete absence of Fermi admixture 
and makes it unlikely that it exceeds 0-005. 


§ 5. CONCLUSIONS 
In the present experiment two possible sources of error which may 
have been present in the earlier **Co alignment experiments of Daniels 
et al. (1952) and of Griffing and Wheatley (1956) have been removed : 


(i) The use of a single crystal rather than a large number of crystals 
eliminates any error arising from temperature inhomogeneity in the 
crystal. 

(ii) The choice of a fluosilicate crystal with its single axis of alignment 
rather than the biaxial Tutton salt removes the necessity of determining 
the position of the tetragonal axes and increases the sensitivity of the 
y-ray method to the parameters of the f-transition. 

We conclude that the Fermi admixture is negligible, the value of y 
being —0-003+0-005. This result is consistent with those obtained in 
the ‘ parity ’ type of experiment. 

Note added in proof.—In subsequent experiments we have shown that 
the anisotropy of the 1-62 Mev y-radiation is closely the same as that of 
the 0-805 Mev radiation. This result confirms the hypothesis that 2’ is 
about zero and is incompatible with a value of \’ greater than 0-2. 
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ABSTRACT 


An expansion equivalent to Taylor’s series is derived for functions of two 
operators of the form #(4-+ B) where B is a shift operator with regard to A; 
the result is obtained as a special case of a theorem due to Kermack and 
McCrea. The development in powers of B involves finite differences of the 
function F(A) instead of derivatives as in the case of commuting variables. 
Exponential functions can be factorized into two terms, each depending on 
one operator only. 

The result allows a Gaussian potential energy to be expressed as a product 
of exponentials involving only creation and destruction operators. Explicit 
results are derived for the matrix elements of such a Gaussian function in the 
representations which diagonalize the Hamiltonians of the harmonic oscillator 
in one and two dimensions; they can be expressed in terms of Jacobi poly- 
nomials (finite hypergeometric series). 


§ 1. INTRODUCTION 


In many problems it is of importance to determine how a function F(A} 
of an operator A varies as A itself is altered. In the field of complex 
numbers this variation is determined by Taylor’s theorem 


Role IP a2r d 
et ne ee ee iG e ital 
1! dx i PAINTS is 33 exp( a) (”) ir) 


but this formula depends on the commutative law of multiplication. 
In consequence it cannot be used to express a function 


DG aa Td Ane) Duane we in ye vara (122) 


in powers of A where A and B are non-commuting operators and A is a 
complex number (a c-number), even if / has a well-defined meaning. 
Any formula analogous to (1.1) involves the whole algebra generated by 
A and B through the continued formation of commutators ; the resulting 
equations either involve repeated integrals (Feynman 1951) or very 
cumbersome commutator expressions (Kubo 1952) ; they gain considerably 
in simplicity in representations in which the operator A is diagonal 
(Schafroth 1951). 


eeeseeeiee 0 ee 
+ Communicated by the Author. Now at the Department of Theoretical 
Chemistry, University of Cambridge. 


F(x+h)=F(x) + 
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A formula for the expansion (1.2) has been derived by Kermack and 
McCrea (1931) which is applicable whenever A and B can be expressed as 
functions of two operators p and q which satisfy the simple commutation 
rule 

lg, pl=qp—pq=1. Ls yee GAP Ee ee 


This work is not well known even though one of their results 
exp (ap +bq)=exp (ap) exp (bq) exp (ab), . . . (1.4) 


where a and 0 are c-numbers, is widely used in field theory. 
One particularly important case of non-commuting operators arises 
when one of them is a shift operator with regard to the other such that 


[4 Bi-4AB- BA=ab . 


The Taylor expansion (1.2) for the case can easily be derived from the 
general formulae of Kermack and McCrea (1931). In § 2 their treatment 
will be briefly outlined and applied to the case of shift operators. 
Questions of convergence will be ignored as they have been fully discussed 
by Kermack and McCrea. 

In certain quantum-mechanical problems it is necessary to consider 
potential energy functions which are essentially quadratic at large distances 
from a given centre, but near this centre have curvatures that do not 
correspond to their asymptotic behaviour and may even be negative. 
Such functions are most easily expressed as the sum of two terms, the first 
of which is purely quadratic and the second Gaussian. It will be shown 
in § 3 that the matrix elements for such a Gaussian function in the represen- 
tation which diagonalizes the Hamiltonian of a harmonic oscillator can 
easily be calculated by repeated applications of the formulae derived in § 2. 


§ 2. THe ExpaNsION THEOREM INVOLVING SHIFT OPERATORS 


The most important of the several expansion formulae derived by 
Kermack and McCrea (1931) can be written formally 


Plg+Mi)1= Faq) + 1 HP) + % HR (q)+...=exp QH)FQ@) (2.1) 


where H is a symbolic operator acting on F(q) 


d id, & 


GI Bo 
Tee ee , 
ti dq e 2!dp dq? 3! dp? dq? ie (2.2) 
= 4 ld, @ ld, a3 
“dq 2! dp dq 3! dp? dq (2.3) 


he suffixes U and 2 attached to f indicate that in each term all the powers 
and derivatives of f(p) premultiply or postmultiply respectively the 


appropriate derivatives of F(q). This expansion is most easily derived if 
use is made of a theorem due to Dirac (1926, 1930) 


at (p)—f(p)q=4f/dp. oA Regen ee See 
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‘Total induction yields for positive integer m 
m — m m af m—1 m af m—2 
q"S(p)=f(p)q (1) Fa +o )aper et (2.5) 
The commutation rules (2.4) and (2.5) also apply if q is replaced by 


w=q+Af, 
so that for any positive integer power w” the first derivative becomes 


a(w") — — MFo yNr—M—L 
oA cas i de 


== n n—-1 n af n—2 ” af n—3 = nv 
@lz + (5) eu + 3 ane “Pee (1) 


The alternative form (2.3) for H is proved in an analogous way ; repeated 
application of (2.6) for a Taylor expansion in power of A leads to (2.1). 

The special commutation rules (1.5) for shift operators can be satisfied, 
in view of (2.4), by the choice 


(2.6) 


A=q, eed) eX 0D) ba eet tel 27) 
In this case we can evaluate the operator H in (2.2) 
H=B, (57 + a +...) = [exp (ad aq) —1}= 2 A,, (2.8) 
where A, is a finite difference operator acting on F(A) 
We AyaF(A +a)— F(A). 2. 2 oes = (2.9) 
The alternative expression (2.3) leads to 
Aho 13510 Nae eee RE Oe. Pees era? (col) 


The suffixes 1 and 2 in (2.8) and (2.10) have the same meaning as in (2.2) 
and (2.3). The complete Taylor expansion for (1.2) becomes with the use of 
(2.1), (2.8) and (2.10): 
F(A +AB)=exp [(AB,/a)A,]F(A) =exp [- (AB,/a)A_, |F(A). (2511) 
For the special case of an exponential function 
F(A)=exp (kA) 
we have in view of (2.9) 
AF (A)=exp [«(4 +a)]—exp («A)=[exp (xa)—1]#(A), (2.12) 
so that (2.11) factorizes: 
exp [x(A +AB)]=exp [(AB/a)(e — 1) ]e4# =e exp [(AB/a)(1—e™)], 
re) 


In the limit a+0, the divided differences A,f'/a become derivatives and the 
expansion (2.11) and (2.13) become the corresponding Taylor series (1.1) 


for commuting variables. 
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§ 3. Matrix ELEMENTS FOR A GAUSSIAN POTENTIAL ENERGY 


A very common representation of dynamical variables in quantum 
mechanics is that which diagonalizes the Hamiltonian of the harmonic 


oscillator 


H =3(p?+q"), pd-Wp=—1; NE OTE ae 


where the units are chosen in such a way that h=w=1. (The variables p 


and q in this section are not the same as in the first two sections ; the com- 
mutation rule (3.1) is not identical with (1.3).) The eigenvalues for this 
Hamiltonian are well known to be 


Ren=i.3, bo 2: ee 


If we try to evaluate the matrix elements for a Gaussian potential energy 
exp (— kg?) in this representation, we are faced with the difficulty that the 
matrices for g¥ become more and more complicated with increasing NV. 
In contrast, the matrices for the powers and exponential functions of the 
creation and destruction operators 


a*=4(q—tp), a=s(qt+up) ... . . (3.3) 
retain a simple structure. In fact, if we use Dirac’s bracket notation for 
the eigenstates of H, the matrix elements for a* and a become 


(nals) = (sla*|n) = 8, n, 1[3(m + 3) P? - + + (3.4) 


where 6,,,,, is the Kronecker symbol 
ne l(m=n) 
mn == O(m x n) 


(cf. Dirac 1947). 
For powers of a* and a the terms become accordingly 


. paar Weed mrs se 
(naa) =(olat[ny = Psat | SEE T. . . 8.5) 


(n—})! 
Furthermore 


(nja*als) =45,,,(n—}), (nlaa*|s)=48,,(n +)... (3.6) 


Considerable simplification could therefore be achieved if the expression 
V =exp (—kq*)=exp [—k(a+o*)?] . . . . (8 


could be replaced by a product involving exponentials of a2, a*2, and 
aa* or a*a. From (3.1) and (3.3) we find the commutation rules 


[@; aS {=o A eee eee 
leg, | =[—a*9,.97 |=92 ae eee 
[ag,0*|= =a? )[a*g,a*?) =a ae en See 


The squares q°, a* and a** therefore act as shift operators with regard to 
the products aq and a*q, and the formulae derived in § 2 are immediately 
applicable. In particular we find from (2.13) and (3.9) 


exp (— Ka*q) =exp [«(aq—q?)]=exp [(1—e*)q?] exp (kag). . . (3.11) 
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‘Similarly from (2.13) and (3.10) 


exp (—«xa*a)=exp [—x(a*q—a**)]=exp (—«a*q) exp [(e*— 1)a*?] 


and (3.12) 
exp (—xaa*)=exp [—«(aq—a*)]=exp [(e*— 1)a?] exp (—xaq). 


(3.13) 
Equating 


k=e*—1, e=log(1+k) . . . . . . (3,14) 
we obtain from (3.11)-(3.13) 


exp (—kq*)=exp (—xa*a) exp (—ka**) exp (— ka?) exp (—xaa*). (3.15) 


Of the four factors on the right-hand side the first and last are diagonal 
matrices in view of (4.6) with elements 


{njexp (—xa*a)|n)=(1+k)-"?+44, (nlexp (—Kaa*)|n) = (1+ hk) 214, 
(3.16) 


The matrix elements of the other two exponentials are, in view of (3.5), 


<s|exp (—ka*?)|n) = (n| exp (—ka?)|s)= Sarl oe AF (3.17) 


[2(s—n)]! L(n—3)! 

provided s—n is a non-negative even number, and zero otherwise. The 
matrix for exp (—kq?) is thus reducible to two matrices involving states 
with n= 21+ 4 and n= 21+ 3 respectively (/=0, 1, 2...). Use of (3.16) and 
(3.17) leads to the final result 
(=y*[2n!2sT" 

(14 &)rts+12 

ib \r+s Lk r+s—2 
x {eo + Lei Tas a (3.18) 


O!lr!s! 2'(r—1)!(s—1)! 
— )rts[(2r + 1)!(2s + 1)! JV? 
(1+k)rts+3e 
(4h)"ts (sk) 2 
TEE C= VERE mere Ore 3.19 
‘ce ce nieeie he ae 
Written in ascending powers of k, this becomes for r > s 


— )r+s[(2r)!(2s)! V2 
<2r+3|V[2s+3)= a 


(akyr-* (hy? 
seen Oe 3.20 
: ea * T1Q@s—2)\r—s+ 1)! ee) 
— )r+s[(2r + 1)!(28+ Tie 

(27 + $|V|2s+ s)= ( al + k)yrtst3l2 

(aky'-° ee e021 
A ae 1)!(r—s)! as IN(2s—1)"(r—s+1)! ee 
For s>r the expressions are equivalent to (3.20), (3.21) with r and s inter- 
changed. 


{2r+4|V|2s+4)= 


(2r + 3|V|2s+3)= ( 


bo 


502 R. A. Sack on 


The derivation (3.8)-(3.15) applies equally in the case of a Gaussian 
potential energy in several dimensions as any creation and destruction 
operator relating to one direction commutes with all the operators relating 
to a normal direction. In particular, for an isotropic two-dimensional 


oscillator 


Happ +0, +d.) toy). 9 eee ee 

the eigenstates can be classified according to their energy nand momentum m 
H\|n,m)=n|n,m), (U2Py—VyPx)|rm)=m|nzm)y . . (3.23) 
m=0,+1,+2..., n= |m|+1,|m|+ 3,... . . . (3.24) 


(cf. Pauling and Wilson 1935). For a potential energy 
V =exp [—k(q.?+ 9,7) ]=exp {— kl (a, +g")? + (ay +a,*) Tf (3.25) 


the only non-vanishing matrix elements occur between states of constant m. 
The matrix elements for (a, +ia,) and (a,* +7a,*) have been derived by 
Moffitt and Liehr (1957) (cf. also Longuet-Higgins et al. 1957); they are 


{n,m|a,—ia,|n+1,m+1)=(n+1,m+I1|a,* +ia,*|n,m) 


= [B(n-+m-+ 1) 7 
(3.26) 


---- 


{n,m|a,* —ia,*|n—1,m+1)=(n—1,m+1|a,+ia,|n,m) 
= [}(n—m—1)}? 


whereas all the other terms vanish. From this follows 


| 
| 
J 


{n,m|a,?+4,?|n+2,m) = (n+ 2, mla,*?+a,*2|n, m) =$[(n + 1)?—m?]}!? 
(3.27) 
{nm|a 4." + dyay*|nm) —$= (nmla,*a,+a,*a,|nm)+$=n. (3.28) 


The results can be expressed somewhat more simply if we confine our 
attention to non-negative m (for negative m the results are entirely ana- 
logous); then the matrix elements of the relevant exponentials become 
with the use of (3.27) 


(2r+m+1,mlexp [—k(a,*? +a,,**)]|2s+m-+1,m) | 


{2s+m+1,mlexp [—k(a,2+a,?)])2r+m+1,m) { 


JCS Ee ~ a (8.29) 


(r—s)! Lsl(s+m)! 
The final result is then, in view of (3.15), (3.25), (3.28) and (3.29), 


_)r+sfplg! 
(2r+m+1,m|V[2s+m+1,m)= (=P lrlst(r + m)i(s +m) PP 


(1 fe k)yrtstm+1 
| are a pakeee her+s—4 
Olmiris! ” Wm + Ir—He—1)l  Bn+2)r—Dl(e—2)! \ 
(3.30 
* (—k)rts (r+m)!(s +m)! 12 
(+ kytstetim! | rtsl oft = Si Lk) ee ae 
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where ,/’, denotes the hypergeometric function (cf. Whittaker and Watson 
1927). ‘The occurrence of non-positive integers as the first two parameters 
ensures the breaking off of the series, the resulting polynomials are known 
as Jacobi polynomials. The expressions (3.18) and (3.19) correspond 
to (3.30), (3.31) if m is taken to be — 4 and 4 respectively. 

These formulae could also have been derived by the use of the generating 
function of Laguerre polynomials (cf. Pauling and Wilson 1935). 
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ABSTRACT 


It is shown that the method of moments can be applied to many problems 
of random walks on lattices with inequivalent sites. The method is illus- 
trated by two simple examples. 


Tue theory of the random walk in one or more dimensions when there is 
no correlation between successive steps is well established and makes use 
of Laplace’s method of characteristic functions for a complete specification 
of the results. The same method applies if the possible positions after 
each step are confined to a lattice, i.e. a discrete set of points, provided 
all the lattice points are equivalent. 

The method of characteristic functions is not applicable in its simple 
form whenever an explicit correlation exists between two or more successive 
steps or for walks on lattices in which the points are not all equivalent. 
In such cases recourse must be made to the more refined theory of Markov 
chains (cf. Montroll 1947, Daniels 1952, Hammersley 1953). 

In many cases a less complete description of the results is required, in 
particular if the only quantities of interest are the mean displacement 
and the mean square deviation therefrom or the parameters specifying 
the asymptotic Gaussian distribution which is obtained after a large 
number of steps. In such cases the use of moments provides a more 
elementary, and usually simpler, method of derivation of the results for 
uncorrelated steps than the method of characteristic functions. For 
motion in one dimension the nth moment m, relating to a single step is 
defined as 


Min = > Def i ae?) eee 


where p; is the probability that the displacement x will be of magnitude 2, 
and H(f) denotes the expectation value for the function f. The nth 
moments [J,,], of the distribution for the algebraic sum 


Xi) = Uy) + LQ) +... Ly) 
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of NV successive displacements a,) satisfy the well-known recurrence rela- 
tions 


[Ma lyas a S (”) [M1 ,-,]ym, ; [Mah =™Mry. yp rte fo (2) 


v= 


This formula is most simply derived from the definition 


Xwiy=A wt My» 


the binomial theorem and the equation 
E(X yy". 4.1)") = E(X yy”). E(t yn") = eee) 


which expresses the independence of the distributions for X(y) and ay). 

The zero’th moments m, and [My], are unity in most applications; 
a smaller value of m, corresponds to the possibility that the moving particle 
is eliminated from the system, the quantity 1—m, is then called the 
probability of absorption. A somewhat different case is found if the p, are 
interpreted as transition probabilities in a fixed time interval; in this 
case there may be a non-zero probability that the particle remains on the 
same lattice point, but m, will be unity if there is no absorption. 

Extension of (2) to random walks in several dimensions is straightforward. 
Thus for a two-dimensional walk we may define 


Mn 2 DRE Yip es ME Te (4) 
1 


and the recurrence formula analogous to (2) becomes 


n 


a w ING = > 


, 


MS elie eee cae) 
v=0 v=0 
the generalization to more dimensions is obvious. 

The purpose of the present note is to point out that the method of 
moments is often applicable even for random walks on lattices with non- 
equivalent sites, provided that steps are uncorrelated in the sense that 
given the starting point of a step the transition probabilities do not depend 
on how this point has been reached. The relations (3) and hence (2) 
remain true for all values of » for which the moments m,(v <n) defined 
in (1) are constants of the lattice, independent of the starting point of the 
step. In particular if the moments up to second order are the same for all 
lattice points the asymptotic distribution after a large number of steps 
can be calculated by the usual methods. Similar considerations apply 
to (4) and (5) for walks in several dimensions. 

The usefulness of the method will be illustrated by two examples: 


(i) A one-dimensional example 
A particle can move on a linear lattice, the distance between neigh- 
bouring points being alternately aand 2a. Withina standard time interval 
the particle can move to either adjoining site or remain where it is. The 


P.M. ee 
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probability of its performing the longer or shorter step are p and 2p: 
respectively, and 1—3p is the probability that it will not move. The 


moments for a single time interval are thus 


mpy=1, m= +(p.2a—2p.a)=0, my=G6pam . . . (6 a) 


for all sites, whereas 


M,= + [p.(2a)? — 2p.a*]= + 6pa* a» 22 (oem 


depends on the starting point. Formula (2) yields immediately 
[Molv=1; [Myly=9; Ly]y =6N pe’. 2 ae 


(ii) Random walk on a hexagonal close-packed lattice 


In the lattice structure each site has 12 nearest neighbours at a distance a. 
The sites are arranged in parallel layers composed of nets of equilateral 


Sites of a hexagonal close-packed lattice. The intersections of the full lines 
and of the broken lines are the normal projections of the lattice sites 
in alternate layers on to a plane parallel to the layers. 


triangles ; the sites of one layer lie above the geometric centres of alternate 
triangles of the adjoining layers, but the lines connecting corresponding 
sites of next nearest layers are normal to the layers. E 

Assume that a particle can carry out an uncorrelated random walk in 
such a way that in a standard time interval it can move to any nearest 
neighbour site with a probability p/6 or q/6 according to whether the 
motion is within a layer or to an adjoining layer; the probability that it 
remains on the same site is 1—p—g. The sites are not equivalent ; in 
rectangular Cartesian coordinates (x, y, z) with the zy-plane parallel to the 
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layers and the x-axis parallel to one set of sides of the triangles the compo- 
nents of the possible steps within the same layer are 


(+a, 0,0) (+4a, +44/3a, 0) 


for all lattice sites, and for steps to adjoining layers (0, a/4/3, +a./%) and 
(+34, —}a/1/3, +a4/3) for sites in even planes and their negatives for 
sites in odd planes. The moments for a single interval are thus in a self- 
explanatory notation 


M,=M,=My,=My,=My,=0 . . . . . (Ta) 
for reasons of symmetry, also 
My= +3q(2.a/4/3—4.4a/1/3)=0, . . . . . (7b) 
Mee = bp(2a?+ 4.40%) + 4g(4.4a%)= (hp + Age, ] 
Myy = EP(4-20*) + 6q(2.50? + 4.50") = (bp + qe, (7c) 
NM = 6q(6.3a*) - 3qa°. 


The inequivalence of sites in alternate layers shows itself only in the third 
and higher moments. For N intervals the non-vanishing second moments 
are thus in view of (4) 


[M n2lw=[Myyly=($9+39)Na*; [M,.ly=3qNo®. . . (8) 


As these results depend only on the products pN and qJ, they are equally 
valid for a continuous process in time in which p and qg denote transition 
rates and N the timet. In view of the Einstein (1905) relation 


[ea 2D fete: ee eee eg ne) 


LL? 


the diffusion constants D,,, D,,,, D,, are given by the moments in (7 c) 
divided by two. The results (8) have been derived by Hammersley (1953) 
by means of a complicated application of the theory of Markov chains. 
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ABSTRACT 


Whiskers of sodium chloride, potassium chloride, copper sulphate and 
sodium nitrate, grown from aqueous solution, have been studied by optical 
and electron microscopy. ‘The results lead to the conclusion that growth is in 
two stages: an initial linear growth and subsequent whisker thickening by 
two-dimensional nucleation. This picture of whisker growth is in accordance 
with a theory which has been proposed by Sears. Certain direct electron 
microscope results suggest that the mechanism of whisker thickening is 
responsible for the known decrease in strength with increasing diameter. 


§ 1. INTRODUCTION 


THERE are two distinct methods by which whiskers can be grown from 
solution. In one method, carried out by Newkirk and Sears (1955) and 
Gordon (1957), the whisker grows totally immersed in the solution. The 
second method, as developed by Gyulai (1948) and Kato (1955), produces 
whiskers connected to the parent solution only at their bases: the solution 
is retained in a porous plate or exists as a film of solution on a large 
crystal. 

One of the most interesting properties of whiskers, their high mechanical 
strength, was demonstrated for the case of sodium chloride whiskers by 
Gyulai (1954). Some of the whiskers of about 1. diameter had breaking 
strengths of about 100 kg/mm? but thicker whiskers had progressively 
decreasing strengths, as their diameters increased. 

In the present work the final shapes and surface structures of whiskers 
of several inorganic salts, grown by the method Gyulai, have been studied 
using the electron microscope. The surface structures can in some 
respects be explained in terms of the growth mechanism. 


§ 2. METHOD orf GROWTH AND OpTicaAL OBSERVATIONS 


Pieces of a refractory boat were soaked in a concentrated solution of 
the salt under investigation, either sodium chloride, potassium chloride, 
sodium nitrate or copper sulphate, and then allowed to dry out very 
slowly at room temperature (20-25°c). The whiskers varied in diameter 


between 0-1 4 and 20 4 and grew several millimetres in length. 
a 


+ Communicated by Dr. G. A. Geach. 
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Observations of sodium chloride whiskers under a low-powered micro- 
scope, showed that growth was preceded by an induction period of about 
a day after which whiskers grew to a length of about 1 mm in 15 hours. 
There was some evidence that, towards the end of growth in length, 
thickening occurred. In many cases the whiskers were seen to grow out 
from small sodium chloride crystals on the surface of the refractory boat. 
Sodium chloride whiskers were usually quite straight and unbranched 
and kinks were rare; occasionally whiskers in the form of right-angled 
crosses were observed. Whiskers of the three other substances investi- 
gated were more frequently kinked and branched. 


§ 3. ELEcTRON DIFFRACTION 


The orientations of sodium and potassium chloride whiskers were 
determined by electron diffraction. Out of 20 sodium chloride whiskers 
examined 19 had [100] orientations and one had a [110] orientation. <A 
similar number of potassium chloride whiskers examined all had a [110] 
orientation. 

Although the electron diffraction specimens were up to 1, thick, 
extended cross-grating patterns were observed and Kikuchi lines were 
absent. It was concluded that the whiskers were not perfect crystals 
but were distorted. 


§ 4. CARBON REPLICA OBSERVATIONS 


Single stage carbon replicas were made of sodium and potassium 
chloride and copper sulphate whiskers, shadowed tan-!3 with gold- 
palladium, in order to study their shapes and surface structures. Most 
of this work was done on sodium chloride whiskers with a smaller amount 
on potassium chloride and copper sulphate. 


4.1. Sodium Chloride 


It was found that the whiskers could be broadly classified into three 
groups on the basis of their size and surface structure. In the first group 
are those whiskers, or sections of whiskers, which have diameters of about 
1 or less. The whisker cross sections were irregular and varied along 
the length; the surfaces were smooth and rounded. Two examples of 
this type are shown in figs. 1 and 2, Pl. 20; in each case the thin section 
is more rounded than the thicker section which is faceted and regular. 

The thick section of the whisker in fig. 2 is typical of the second group 
in which the whisker cross sections were regular and constant for con- 
siderable lengths. These whiskers were typically a few microns in dia- 
meter. A square cross section bounded by {100} planes was most 
common but occasionally cross sections which were apparently triangular 
were seen as in fig. 3, Pl. 20. It is possible that the triangular whisker 
had a [111] orientation. 

The most commonly observed whiskers, in contrast to the groups 
already discussed, showed many growth layers. These whiskers were of 
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a lath-shaped cross section between 1 and 10, in width and with width 
to thickness ratios up to 3:1. Two parallel faces were {100} planes 
(developed to various extents) and the whisker cross section was com- 
pleted by the irregular termination of growth layers. Typical examples 
of this type are shown in figs. 4 and 5, Pl. 20; in fig. 5 many fine lines can 
be seen on the side of the central ridge indicating that this part of the 
whisker has been built up by layer growth. The growth layers observed 
were very variable in height and, where it could be decided, they were 
initiated at the centre of a face, as in fig. 6, Pl. 20. 


4.2. Potassium Chloride 


A large proportion of these whiskers were smooth thin laths of ir- 
regular width a few microns wide and a fraction of a micron thick, as in 
fig. 7, Pl. 21; they only exhibited faint growth markings. They were 
frequently kinked; in some cases at angles of 60° corresponding to the 
angle between (110) growth directions. Another type was a more 
regular, straight, trapeziform whisker as shown in fig. 8, PI. 21. 


4.3. Copper Sulphate 


The most commonly observed whiskers of copper sulphate were only 
a few tenths of a micron in diameter and appeared to be square in cross 
section as in fig. 9, Pl. 21. Only short lengths of such whiskers were 
observed, probably as a result of fracture in the preparation of the 
replicas. The tangled appearance of copper sulphate whiskers, as seen 
in the optical microscope, is probably a result of many thin fibres which 
coalesce at various points; some support for this idea is given by the 
replica shown in fig. 10, Pl. 21. Since the whiskers were always sensibly 
straight, as seen on the electron micrographs, they must break in a 
brittle masiner and do not kink like sodium chloride whiskers (Gyulai 
1954). 

Occasionally thick copper sulphate whiskers were seen and in these 
cases they were irregular and rounded and exhibited growth layers. 


§ 5. WHISKER DISINTEGRATION IN THE ELECTRON BEAM 


When sodium and potassium chloride whiskers were subjected to an 
intense electron beam in the microscope interesting effects were observed. 
In fig. 11, Pl. 21 is shown the result of irradiating a sodium chloride 
whisker about 7 wide; the whisker has almost disintegrated and a 
pattern of rectangles of different opaqueness has been revealed. Effects 
similar to fig. 11 were always seen during the disintegration of sodium or 
potassium chloride whiskers providing they were greater than about 1 pu 
in width. Although the whisker axes were different for these two sub- 
stances (§ 3) the rectangles were always oriented parallel to [100]. 

Electron irradiation of sodium chloride whiskers a few tenths of a 
micron in width never resulted in a regular pattern during disintegration. 
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Figure 12, Pl. 21 is typical of our observations on this class of whisker; 
although the transparency is not uniform there are no clear lines of 
_ demarcation and there is no pattern. 


§ 6. Discusston 


For all three kinds of whiskers investigated by electron microscopy a 
general picture of growth in two stages seems satisfactory. Initially there 
is a rapid linear growth of a thin whisker and later growth results in a 
thickening of the whisker by layer growth. In the case of sodium 
chloride, for which a more detailed picture has been obtained, the initial 
thin whisker is irregular and rounded in cross section and layer growth 
develops {100} planes to varying extents. The difference between 
whiskers of the second and third groups can be attributed to the different 
thickness of the growth layers resulting from small differences in con- 
ditions of growth. 

Gyulai (1954) and Kato (1955) reported that whiskers of this type grow 
from the tip, so that it is probable that the initial growth in length is a 
result of one or more axial screw dislocations. Some support for this 
may be found in the absence of Kikuchi lines in the electron diffraction 
patterns; although misorientation as a result of layer growth could also 
extinguish Kikuchi lines it is probable that most of the electron diffraction 
specimens had not thickened to any extent. 

The observations of the way in which whiskers thicken support the 
theory of growth proposed by Sears (1957). On this theory solute is 
transported to the whisker tip by the hydrodynamic flow of a thin ad- 
herent film of solution; the driving force for the flow is the gradient of 
surface tension resulting from a gradient in the concentration of the 
adherent film. As the solution travels up the whisker from the base its 
supersaturation steadily increases, falling again at the tip where the 
supersaturation would remain small. The maximum supersaturation 
along the whisker increases as the whisker grows in length until two- 
dimensional nucleation can occur on the whisker surface, resulting in a 
thickening of the whisker by layer growth. The onset of layer growth 
will inhibit the further growth in length of the whisker. In contrast to 
the growth of lead nitrate crystals by layer growth (Williams 1957) two- 
dimensional nucleation on sodium chloride crystals occurs at the centre 
of a face and not at the edges. This difference, however, may merely 
be a result of considerable differences in growth conditions. 

The possibility that the observed growth layers arise from dislocations 
which have a screw component normal to the whisker axis, has been 
rejected for two reasons: firstly, the number of centres of nucleation is 
sometimes improbably high (> 108/cm?); secondly, the observation of 
whisker crosses, in the case of sodium chloride, can be interpreted as due 
to screw dislocations formed at the growing tip giving rise to two collinear 
branches. Thus screw dislocations transverse to the axis do not give 
rise to layers but branches during the growth of sodium chloride whiskers. 
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The thickness of the adherent film of solution, necessary for the whiskers. 
to grow at the observed rate, can be calculated using the equation derived 
by Sears}. Taking the diameter of the growing whiskers to be 1 yu and 
supposing that the whiskers grow 1 mm in 15 hours, the film thickness is 
about 60 A, which is not an unreasonable figure. 

It is not proposed to attempt a complete explanation of the mechanism 
by which whiskers disintegrate in an intense electron beam; it is probable 
that evaporation takes place from the crystal surface which the electrons 
enter. Nevertheless, the size effect observed in these experiments, which 
corresponds to that whisker diameter below which no layer growth could 
be expected, is significant. The correspondence implies that layer growth 
creates a structure in the crystal which is not produced by the initial 
linear growth. Since sodium chloride whiskers only have a high strength 
when they are of about 1 in diameter it is suggested that the structure 
created during layer growth is responsible for a lowering of the mechanical 
strength. 
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+ The equation as given by Sears (1957) is somewhat incorrect. TI i 
for the film thickness 7 should be: as 
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where / is the length of whisker of radius 7 which grows in a time f:; p is the 
density of the whisker; 7 is the viscosity of the saturated solution: M is the 
molecular weight of the solute; Cy is the equilibrium concentration of the 
solution; o* is the critical supersaturation for two-dimensional nucleation and 
dy/dc is the rate of change of surface tension with concentration. 
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ABSTRACT 


The number of slip bands formed during the fatigue of O.F.H.C. copper 
at 90°K and 293°K have been compared. For equivalent stresses to cause 
failure in 2 x 10° cycles, more slip bands are formed at 90°x than at 293°K. 
The number of slip bands formed increases with the length of the test. 
When specimens are fatigued at 90°K, and heated to room temperature 
at intermediate stages, the rate of increase of slip bands is much smaller 
than for specimens kept entirely at 90°K. The intermediate heating 
effectively makes the slip bands softer than the surrounding regions so that 
on further fatiguing slip is localized near or within the slip band. The 
critical temperature for the intermediate heating, above which the annealing 
occurs, is approximately 180°x. The results are attributed to point 
defect hardening during the annealing treatment. 


§ 1. INTRODUCTION 


McCammon AND ROSENBERG (1957) have found a marked temperature- 
dependence of the fatigue life of metals. They reported that the stresses 
required to fracture copper in 2 x 108 cycles at 293, 90, 20 and 4-2°K were 
1-2, 1-75, 2-4 and 2-55 x 10° dynes cm-? respectively. It is well known 
that in fatigue large slip movements to and fro are localized in slip bands 
and that eventually fatigue cracks are formed in these bands. The 
temperature-dependence of fatigue life is probably related to the changes 
with temperature of the ease of formation of slip bands and the amount 
of slip in each band. In an investigation into the appearance of slip bands 
formed at different temperatures (Hull 1958) it was observed that slip 
bands became finer and more closely spaced as the temperature of 
fatiguing was decreased. Closer study of the difference in the number of 
slip bands in specimens fatigued at 293 and 90°x, after different lengths 
of time, indicated that some changes occurred in the slip bands when 
specimens fatigued at 90°K were heated to 293°K for examination. It 
was found that fewer slip bands formed on subsequent fatiguing at 90°K 
than were expected. This effect has now been examined in detail and 
the results are described here. 


§ 2. EXPERIMENTAL 


The fatigue tests were made in the apparatus used by McCammon 
and Rosenberg (1955). Specimens of O.F.H.C. copper, approximately 
rent ces, ree 
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1/32 in. thick, were shaped as shown in the figure (a). The specimens were 
annealed at 600°c for 4 hours in vacuo and then electropolished before 
fatiguing. 

The number of slip bands formed on the specimen after fatiguing was 
measured by determining the distribution of slip bands over the whole 
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: : 2 : , 
(a2) Shape of specimen used illustrating the way in which the surface was 
scanned to determine the slip band distribution. The number of slip 
bands was counted in the microscope at all the points marked (X). 


(b) A typical slip band distribution curve of a fatigued specimen. 


specimen. Counts were made of the number of slip bands in one quadrant 
of the eyepiece of an optical microscope (magnification x 500, area of one 
quadrant 3x 10~4cm?). This was scanned over the specimen in the way 
illustrated in the figure (a). The average of the three results across the 


———— 
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specimen was plotted against the length of the specimen. The figure (b) is 
a typical distribution of slip bands. The value Smax, marked in the 
figure (b), was taken as a measure of the number of slip bands formed. 

The specimens were fatigued at 293 and 90°x. After fatiguing at 90°K 
for 104 cycles some specimens were annealed at temperatures between 293 
and 140°K and then refatigued for a further 2 x 104 cycles at 90°x. The 
annealing was carried out by immersing the specimens in a bath of pre- 
cooled acetone. The temperature could be kept constant to within 3 
degrees for over an hour. The specimens had to be removed from the 
fatigue apparatus for examination, but could be replaced without dis- 
tortion. All the annealing treatments, except those at 293°K, were made 
im situ. 


§ 3. RESULTS 


The O.F.H.C. copper used in the present work showed the same 
dependence of fatigue life on temperature as the H.C. copper used by 
McCammon and Rosenberg, but the actual fatigue curves were slightly 
lower. For a life of 2 x 10% cycles the stresses required for failure at 293, 
90, 20 and 4-2°K were 1-0, 1-6, 2-2 and 2-4 10° dynes em~?; when these 
same stresses were applied in unidirectional tension, then 1-0 x 10° dynes 
cm? produced an extension of 1-9°% at 293°K, 1-6 x 109 dynes cm-? an 
extension of 5-3°% at 90°K and 2-2 x 10° dynes cm~? gave an extension of 
8% at 20°. 

When a specimen was fatigued at 293°K at a stress of 1-0 x 10° dynes 
em? for 104, 3x 104 and 10° cycles the value of Smax increased pro- 
gressively as shown in the table. This was repeated at 90°K with a stress 
of 1-6 x 10° dynes cem-?, but in this case the specimen had to be warmed 
to room temperature to be examined after each fatigue (i.e. after 104, 
3x 104 and 10° cycles) and then recooled. The table shows that after 
104 cycles at 90°K many more slip bands have been formed than with an 
equivalent stress at 293°K. After 3 x 104 and 10° cycles there were more 
slip bands, but the differences are less pronounced. To avoid the effect 
of heating the specimens to room temperature at intermediate stages of 
the fatigue test for examination, three separate specimens were fatigued 
at 90°K for 104, 3 x 104 and 10° cycles respectively. The increase in Smax 
(see table) are quite different from the increases that occur with inter- 
mediate heating; the interruption of the fatiguing for room temperature 
examination has resulted in a marked decrease in the number of slip bands 
formed subsequently. The changes are due probably to a thermally 
activated process. To determine the temperature range in which this 
occurs a number of specimens were fatigued at 90°K for 10* cycles and 
then annealed at temperatures between 190 and 140°K for one hour, 
followed by a further fatigue of 2 x 104 cycles at 90°K. The results which 
are summarized in the table show that for annealing temperatures below 
167°xK, the value of Smax is approximately the same as for specimens 
given 3x 104 cycles without an anneal. For annealing temperatures 
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above 190°K Smax is approximately the same as for specimens heated to- 
293°k. The annealing effect occurs, therefore, between 167 and 190°K. 
Annealing at 393°K has the same effect as an anneal at 293°K. 

The intermediate annealing of specimens fatigued at 90°K also affects 
thé thickness of the slip band formed. At room temperature, fatigue 


Values of S,,,,, After Different Fatigue and Annealing ‘Treatments 
Number of cycles. 104 3x 104 10° 
Fatigue at 293°K 5 13 26 


Fatigue at 90°K on one specimen 
warming to 293°K for examination 
between each fatigue period 19 26 = 


Fatigue at 90°K using three different 
specimens, one for each fatigue 
period 19 66 ~ 90 


Fatigue at 90°K to 10* cycles. Anneal 
at 148°xK and then fatigue to 3 x 10# 


cycles —- 65 = 
Asabove. Anneal at 157°K _ 68 — 
As above. Anneal at 167°K = 69 — 
As above. Anneal at 190°K “— 31 —- 


Fatigue at 90°K on one specimen 
heating to 393°K between each 
fatigue period and then examining 20 


bo 
(or) 


36 


produces a few broad slip bands which gradually thicken by formation of 
further slip steps adjacent to the band. A number of new bands are 
formed between these broad bands. At 90°K a large number of fine 
bands are formed and there is more emphasis on formation of new bands 
than on thickening of existing bands. However, if the specimen is 
warmed to 293°K after fatiguing at 90°K further fatigue at 90°K causes. 
the existing bands to thicken. 


§ 4. Discussion 
Comparing the effect of increasing the length of the fatigue test on the 
number of slip bands formed, we see that at 293°K a relatively small 
number of bands form and that these become broader and increase in 
numbers slowly. At 90°K a larger number of bands form which do not 
broaden appreciably, and new bands form between the existing bands. 
The essential difference appears to be that at 293°x further slip can occur 
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close to existing bands throughout the test, but at 90°x the bands rapidly 
become resistant to further slip and new bands form between them. 
The larger number of slip bands formed at 90°K than at 293°K at the 
beginning of the test (10* cycles) is due partially to the larger amount of 
plastic strain, as shown by the tensile tests, in specimens fatigued at 
1-6 x 10° dynes cm=* at 90°K than in specimens fatigued at 1-0 x 10° dynes 
em™ at 293°K. The results resemble closely those of Brown (1951) 
for aluminium. He found that in direct tension many more slip bands 
were formed at 90°K than at 293°K and 543°K. At the higher tempera- 
tures further slip occurred in planes adjacent to the existing bands far 
more readily than at 90°k. This is most probably due to cross slip 
which is a thermally activated process occurring when specimens are 
deformed in the parabolic hardening or Stage III region of the stress— 
strain curve. More cross slip occurs in the existing slip bands at 
higher temperatures and so produces broader bands. In the fatigue 
experiments a constant cyclic stress was applied in a range where cross 
slip is expected and this may be responsible also for the formation of 
broader bands at 293°K than at 90°K in this case. 

The annealing treatment after fatigue at 90°K effectively makes the 
existing slip bands softer relative to the surrounding material so that on 
further fatigue at 90°K slip remains localized in these slip bands. This 
can be achieved either by hardening the regions between the slip bands, 
which would produce an overall increase in hardness of the specimen, or 
by softening the regions close to and within the bands, which would 
result in softening of the specimen. Broom and Ham (private communi- 
cation) have recently observed that a specimen fatigued at 90°K and then 
heated to room temperature had a higher yield stress when extended in 
tension at 90°K than a specimen which was fatigued at 90°K and extended 
at 90°K without intermediate heating. The observation indicates that 
hardening of regions between the slip bands occurs on annealing. The 
hardening between the slip bands will be an additional effect to cross slip 
in determining the width of the bands formed at different temperatures. 
Broom and Ham (1957) reported that more fatigue hardening occurs at 
293°K than at 90°K which suggests that slip will be more localized into 
bands at the higher temperatures. 

Electrical resistivity’ measurements at low temperatures by MacCrone, 
McCammon and Rosenberg (private communication) show that con- 
siderable annealing of the extra resistivity produced by fatiguing copper 
at 90°K occurs in the region of 170°K. This is in the same temperature 
range at which the present effects are first observed. Since electrical 
resistivity is relatively insensitive to dislocation regrouping compared 
with the rearrangement of point defects (Nowick 1955), it is probable that 
the annealing effects are due to point defects. 

Hardening of regions between the slip bands could result from diffusion 
of defects in the lattice to free dislocations, thus pinning them. It is 
unlikely that at temperatures as low as 167—-190°K vacancies will move 
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more than a few atomic spaces, but this may be sufficient to take them to 
dislocations. More complex defects, such as divacancies (Henderson and 
Koehler 1956), are expected to be more mobile than single vacancies. 

In conclusion, the difference in the number and width of slip bands 
formed at 293 and 90°K is due probably to the larger amounts of cross 
slip and greater hardening at the higher temperature. The effect of 
heating to above 167 to 190°K, after fatigue at 90°K, on the number of 
slip bands formed on subsequent fatigue at 90°K is thought to result from 
point defect hardening particularly between the existing slip bands. 
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ABSTRACT 


Neutron irradiation, gamma irradiation and quenching all change the 
internal friction of pure polycrystalline rods of copper. Only if the rods are 
initially given an anneal near the melting point for many hours do these 
three treatments give consistent and similar results, and then the logarithmic 
decrement is reduced. It is concluded from the results that it is the vacancies 
introduced by each of the three treatments which pin the dislocation lines 
by diffusing to them, and that unless the purity of the copper is high their 
effect can be obscured. 


§ 1. INTRODUCTION 


The internal friction of pure metals in the kilocycle per second frequency 
range is drastically reduced, and the Young’s Modulus correspondingly 
increased, by the introduction of small quantities of impurity (Weertman 
and Salkovitz 1955, Takahashi 1956, Weinig and Machlin 1956). Recently 
irradiation and quenching have been shown individually to produce a 
similar reduction in the internal friction and an increase in the Young’s 
Modulus. Electron bombardment of polycrystalline copper at — 180°C 
increases the Young’s Modulus, a dose of 4:5 x 101° 1 Mev electrons/cm? 
producing almost the maximum (i.e. complete) change obtainable for 
large doses, if there is an intervening room temperature anneal (Dieckamp 
and Sosin 1956). Neutron or gamma irradiation of copper single crystals 
produces comparable changes in the Young’s Modulus with 10 neutrons/ 
em? and 1015 gamma photons/em? respectively (Thompson ef al. 1955, 
Thompson and Holmes 1956, 1957). 'The internal friction of gold decreases 
when it is cooled in air from 900°C (Roswell and Nowick 1957) and the 
amplitude dependence is reduced. Levy and Metzger (1955) report that 
after quenching from 640°C, the decrease in the internal friction of alum1- 
nium occurs slowly at room temperature. 

This paper describes some experiments designed to compare the effect 
of three different ‘ treatments ’"—neutron irradiation, gamma irradiation 
and quenching upon the internal friction of polycrystalline copper. 
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§ 2, EXPERIMENTAL METHOD 


The apparatus used was essentially the same as that described by 
Forster (1937). Each sample, in the form of a rod 0-5 cm diameter and 
7-5.em long, was hung horizontally in two cotton stirrups placed 1 cm 
from each end, i.e. about 0-6cm outside the nodes when the rod was 
vibrating in its fundamental mode. One of the cotton stirrups was used 
to induce the rod to vibrate transversely and the other to detect the 
oscillations. Calibration of the detection system (consisting of the cotton 
stirrup, the piezo-electric pick-up, the head amplifier and the main ampli- 
fier) enabled the strain amplitude of the rod to be estimated. 

The advantage of this apparatus over other types lies in the ease of 
-demounting and replacing rods (sometimes radioactive) with the minimum 
of deformation. This was an essential requirement as the rods had to be 
removed from the apparatus either for irradiation or quenching, and 
carelessness in handling the rods was found to alter their internal friction. 

The rods normally have two resonances, up to 100 ¢/sec apart, which 
have been ascribed to the maximum and minimum moments of inertia 
of cross sections of the rod (Russell and Boag 1953). Each rod was always 
rotated in its stirrups until the amplitude of vibration at one of these 
frequencies was dominant. The rod was then always measured while 
so oriented in the stirrups so that this resonance was used. The logarith- 
mic decrement (A) was determined by measuring the frequency width 
of the resonance peak (dv) at half the amplitude at the resonant frequency 
(v= 3300 c/s approximately) A= (z/+/3)(dv/v). 


§ 3. EXPERIMENTS 


The material used throughout the investigation was Johnson Matthey 
spectroscopically standardized copper which has a stated impurity content 
of about 8 parts per million of lead, silver, iron and nickel together, and 
less than 5 parts per million of other elements spectroscopically detectable. 
This material was received in the worked state, when it had the values of 
internal friction almost independent of strain shown as curve (a) in fig. 1. 

Several rods were annealed in vacuo successively for one hour at the 
temperatures 200, 400, 600, 800 and 1000°c. No appreciable change in the 
logarithmic decrement occurred until the last of these anneals when the 
values shown in fig. 1 (b) were obtained. These values compare with 
those for polycrystalline copper used by other workers. 

Some preliminary experiments were performed with rods initially 
annealed for one hour in vacuo at 1000°c. It was found that whereas 
their decrements decreased slightly upon neutron irradiation, both gamma 
irradiation and quenching produced slight increases in decrement. The 
results, however, were not consistent from one rod to another. It was 
thought that the low initial decrement of these rods and the inconsistent 
results were due to impurities in the copper. 

Eventually an annealing procedure was evolved which consistently 
gave rods with logarithmic decrements higher than that shown as curve (c) 
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infig.1. This consisted of initially acid etching the rods and then annealing 
them for 48 hours at 1050°c ina silica tube vacuum furnace. In an attempt 
to avoid the introduction of impurities the rods were placed upon a rack 
made of spectrographically pure graphite rods. During the anneal 
the pressure was maintained at less than 10-5 mm of mercury. The rods 
were allowed to cool slowly with the furnace. The grain size of the copper 
after this treatment was about 0-3cem. By performing several high 
temperature anneals on a given rod it was possible to obtain logarithmic 
decrements as high as that shown as curve (d) in fig. 1. 
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Logarithmic decrement versus strain amplitude for the copper rods before and 
after various heat treatments. 


One rod which had been annealed in this way was etched in stages in 
dilute nitric acid to two-fifths of its original diameter. After each reduction 
in diameter the logarithmic decrement was measured. There was no 
change in the values, suggesting that the internal friction was homogeneous 
throughout the rod. 

Some rods were annealed at temperatures up to 500°C in vacuo and 
slowly cooled. Their values of internal friction were retained. 

Rods which had received identical initial anneals giving them high 
values of logarithmic decrement were irradiated with neutrons or gamma 
rays or quenched and the effect upon their internal friction compared. 
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§ 4. REesuLtTS 


Figure 2 shows graphs for three of the polyecrystal rods which were used 
to investigate the effect of neutron irradiation (number 77) gamma irradia- 


tion (number 75) and quenching (number 81). For each rod the graphs: 


show, firstly the initial values of logarithmic decrement versus strain, 
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Logarithmic decrement versus strain amplitude for three copper rods number 77 


irradiated with neutrons, number 75 irradiated with gamma rays and 
number 82 quenched. 


amplitude, secondly the final values after the treatment and lastly the 
values after an anneal for 2 hours at 400°c in vacuo. As soon as possible 
after the irradiation or quenching treatment the rods were replaced in the 
cotton stirrups and the logarithmic decrement at a low strain amplitude: 
measured at intervals. These values for rod numbers 77 and 75 are plotted 
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versus the time after treatment in fig. 3. A similar graph for rod number 55 
(quenched to room temperature in silicone fluid from 530°C) is given for 
comparison. 
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Logarithmic decrement measured at small strain amplitude versus time after 
(a) neutron irradiation (6) gamma irradiation and (c) quenching. 


The individual treatments given the rods were as follows : 


Rod Number 77. Neutron Irradiation 
The rod was gently placed in an aluminium container packed with 
cotton wool which was then lowered into a vertical hole in BEPO where it 
was held for 15 minutes and then removed. The rod received a dose of 
4 x 104 thermal neutrons/cm? while still at room temperature (it is assumed 
that it received an equal dose of fast neutrons); it was not in BEPO long 
enough to take up the temperature of the reactor (about 60°C). 


Rod Number 75. Gamma Irradiation 
This rod, also carefully packed, was placed 2 cm from a 300 curie ®Co 
source and left there for 168 hours. It is estimated to have received a 
dose of 6 x 1016 gamma photons/cm? while at room temperature. 


Rod Number 81. Quenched 

This rod was heated to 1044°c for 8 minutes while held in a molybdenum 
cage. It was then gently lowered from the furnace to cool for 3 minutes 
in argon and then in silicone fluid at room temperature. | 

The critical condition at which each of the three treatments produced a 
complete decrease in the logarithmic decrement was estimated to be 
approximately (a) a dose of 3 x 10!4 fast neutrons/cm?; (b) a dose of 101” 
gamma photons/cm?} and (c) a quench in silicone fluid from a temperature 
of 550°C. 

+ It should be noted that the gamma photon dose given rod 75 (see fig. 2) 
was less than this critical value and that the minimum value of logarithmic 


decrement has not been reached. 
2N2 
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The temperature at which the logarithmic decrement recovered was 
determined by performing a series of anneals. Those rods which had 
received the above critical treatments all recovered at about 400°c. 
For conditions less than critical, the values recovered at temperatures 
as low as 100°c, while for conditions above critical the recovery temperature 
was above 400°C. 


§ 5. Discusston 


At frequencies above a kilocycle per second the estimated damping due 
to reorientation effects, thermoelasticity and grain boundary sliding is 
much smaller than the value which is observed and, since grain boundary 
sliding is of little importance, single crystal and polycrystalline material 
of the same purity will have very similar values of logarithmic decrement. 
The internal friction of metals measured at these frequencies is considered 
to be due to the motion of dislocation lines within the crystals. The 
oscillatory motion of these dislocation lines (and hence the internal friction) 
will be very sensitive to the free length of the dislocation loops, which will 
normally be less than the distance between dislocation nodes due to the 
presence of pinning agents along the dislocation line. The breakaway 
of the dislocations from these pinning agents at large stress amplitudes 
has been used to explain the marked amplitude dependence of the internal 
friction which appears above a critical stress amplitude (Nowick 1950, 
Marx and Koehler 1950, Granato and Liicke 1956). Any discontinuity 
in the lattice which will attract a dislocation line to it can act as a pinning 
agent. Thus impurities, interstitial atoms and vacancies can act in this 
way. For the free length to be greater than 10-> cm the atomic concentra- 
tion of impurities on dislocation lines must be less than 10-§. Thus, in 
even the most pure materials, the free length of the dislocation loops will be 
limited. 

The impurities in the copper used initially obscured the effects due to 
irradiation and quenching presumably because of the heavy pinning by 
impurity atoms and by their trapping of point defects. 

Experiments have shown that an increase in the impurity content 
increases the critical stress at which the internal friction becomes amplitude 
dependent and reduces the value of the amplitude independent part 
(Weertman and Salkovitz 1955, Takahashi 1956). Thus the changes 
brought about by the high temperature vacuum anneal, i.e. a decrease 
in the critical stress and raising of the amplitude independent part (see 
fig. 1) correspond to a reduction in the impurity content of the copper. 

The curve (d) of fig. 1 shows values of the logarithmic decrement at a 
given strain amplitude comparable with those given for single crystals. 
The initial values of the internal friction of copper used by various workers 
have been listed by Kamentsky (1956) who showed that the values of the 
logarithmic decrement at low strain amplitude range from 7x 10-® to 
4-8 x 10°, the higher values being generally associated with single crystals 
and with the higher purity material. The higher internal friction of 
single crystal copper is probably due to the removal of some of the impurities 
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from the material by the single pass zone-refining which occurs during the 
Bridgman method of single crystal growth. 

Removal of the surface layers (even to a depth of 0-15 em), does not 
alter the logarithmic decrement of a rod. This suggests that any impurity 
extracted during the high temperature anneal had been removed homo- 
geneously. Thus its diffusion coefficient in copper at 1050°c must be 
greater than 107‘ cm?/sec. However, spectroscopic analysis before and 
after the anneal at 1050°c did not detect any change in the metallic 
impurities. This indicates that the anneal is not reducing the impurity 
content of the copper but it is possible that the impurities aggregate so 
that their spacing along dislocation lines is greatly increased. 

Retention of the high value of damping during the subsequent anneals 
at temperatures up to 500°C suggests that the high damping does not result 
from some unstable configuration of the dislocation lines. 

Gamma irradiation and quenching can hardly introduce impurities into 
the copper rods, and the atomic concentration of zine introduced by the 
thermal neutrons is less than 10-§. The internal friction after both neutron 
and gamma irradiation and subsequent annealing at 400°c is higher than 
it was initially, which suggests no increase in the impurity content. 

The similarity in behaviour of rods which have been quenched to those 
which have been irradiated indicates that the same process is producing 
the effect with all three treatments. Dislocation lines could be produced 
by internal stress during the quench and it has been suggested that neutron 
bombardment can lead to the introduction of dislocations ; however, it is 
most unlikely that gamma photons could do this. Whereas interstitial 
atoms and vacancies must be produced in equal numbers by both the 
irradiation treatments the number of interstitial atoms produced by the 
quench will be 20 orders of magnitude less. The only fundamental process 
which is common to all three is the introduction of vacancies. 

The internal friction of copper single crystals measured at —195°c 
was not altered by 15 minutes neutron irradiation while at —195°c. 
Only after pulse annealing for several hours at 23°c did the internal friction 
decrease (Barnes and Hancock 1958). This experiment suggested that the 
number of interstitial atoms produced (which might be expected to be 
mobile at — 195°c) was not modifying the value of internal friction and that 
whatever agent lowered it could only do so while the temperature was 
above 0°c. 

The movement of vacancies to the dislocation lines can explain the 
decrease of damping. Vacancies introduced either by quenching or 
irradiation will diffuse to dislocation lines at room temperature and pin 
them. The decay of the logarithmic decrement observed at room tem- 
perature after the treatments suggests that the vacancies introduced 
diffuse slowly to the dislocation lines at room temperature. 

The atomic concentration of vacancies produced by the critical conditions 
for each treatment can be estimated. These values are approximately 
(a) 10-7 for neutron irradiation, assuming that each primary knock-on 
produces 100 permanently displaced atoms () 10-7 for gamma irradiation 
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assuming a cross section for the displacement of copper atoms of 1 barn for 
the gamma photons} and (c) 2 x 1077 for the quench assuming an activation 
energy of formation of vacancies of 1-1 ev and that all these vacancies are 
retained during the quench. There is remarkable agreement between 
these figures, which adds to the evidence that vacancies are responsible 
for the pinning. These figures suggest that the dislocation damping is 
greatly reduced only when the atomic concentration of vacancies in the 
material is about 10-7, i.e. when there is one vacancy for each atomic site 
along the dislocation lines in the crystal if the density of dislocation lines 
is 10°/cm2, 
§ 6. CONCLUSIONS 

The main source of variation in the logarithmic decrement of different 
rods appears to be their purity, very high values being associated with pure 
material. It is essential that very pure material be used if the dislocation 
pinning behaviour of agents other than impurities is to be studied. 

The similarity in the internal friction behaviour, after either neutron 
irradiation, gamma irradiation, or quenching suggests that the pinning 
agent produced by all three is vacancies. This view is upheld by the 
suppression of an internal friction change when irradiated at — 195°C 
and the slow decrease in the decrement observed at room temperature 
after each treatment. 
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ABSTRACT 


The internal friction of copper single crystals remains unaffected by a 
short neutron bombardment at —195°c if no warning is allowed before 
measurement. Successive pulse anneals produce no marked change until 
23°c is reached, when the internal friction is greatly reduced. 


§ 1. [INTRODUCTION 


‘THE internal friction of metals normally decreases rapidly during neutron 
bombardment at room temperature, the fast neutron dose associated with 
10 thermal neutrons/cm? in BEPO being sufficient to complete the change 
in polycrystalline copper (Barnes e¢ al. 1958). Thompson et al. (1957) 
have recently reported that the internal friction of a copper single crystal 
is reduced during irradiation, even at about 21°K, but at this temperature 
the change was not completed after a dose of 2 x 10! fast neutrons/cm?. 
An experiment designed to investigate the annealing behaviour of copper 
single crystals after neutron irradiation at liquid nitrogen temperature, is 
‘described here. 
§ 2. EXPERIMENTAL DETAILS 


The two single crystals (referred to as crystals 1 and 2) were grown by 
the Bridgman technique from Johnson Matthey spectroscopically 
standardized copper and were 8-5 cm long and 0-5 cm in diameter. Their 
internal friction was measured using an apparatus essentially the same as 
that described by Forster (1937) and shown in block diagram form in 
fig. 1. The single crystal rod was suspended horizontally by two cotton 
threads tied permanently near to each end. One of these cottons was 
used to induce the rod to vibrate transversely and the other to detect the 
oscillations. Measurements at —195°c were made by surrounding the 
crystal and its supporting threads in a metal jacket, that part of the 
jacket around the crystal itself being immersed in liquid nitrogen. The 
jacket was evacuated to avoid ice formation on the rod and threads, and to 
prevent the condensation of atmospheric oxygen. The logarithmic decre- 
ment (A) was determined by oscillating the crystal in its fundamental mode 
and measuring the frequency width of the resonance peak dv at half the 
amplitude at the resonant frequency v, A=(7/V/ 3)(Sv/v). The strain 
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amplitude of the crystal during oscillation was estimated by calibrating 
the detection system (consisting of the thread, the piezo-electric crystal 
pick-up, the head amplifier, and the main amplifier). 


Fig. 1 


PIEZO- 
VIBRATION ELECTRIC MAIN 
TOR PICK-UP 
ph don AMPLIFIER 


DECADE WITH 
OSCILLATOR ATTENUATOR 
GAIN 102-108 


OSCILLOSCOPE 


Block diagram of the apparatus used to measure the internal friction and 
Young’s modulus of the single crystal rods. 


The logarithmic decrements of both crystals were measured at room 
temperature at various strain amplitudes up to about 10~7. Crystal 1 was 
then deformed in torsion (maximum strain 2 x 10~-?) and remeasured after 
several days when the values had stabilized. After measuring the 
logarithmic decrement of both crystals at —195°c they were gently 
transferred by their attached threads into a Dewar flask where they were 
completely immersed in liquid nitrogen. This flask was gently lowered into 
a vertical hole in BEPO (measured thermal neutron flux 1-2 x 10!*/cm?/sec) 
and held there for 15 minutes. It was left for several hours after 
removal to allow the radioactivity to decay. Crystal 1 was then quickly 
transferred from the flask back into the measuring apparatus which had 
already been cooled with liquid nitrogen. To prevent the crystal heating 
during this transfer, a loosely fitting aluminium tube had been placed 
around it before the first measurement; the suspending threads passing 
through two holes in the tube. The logarithmic decrement at —195°c 
was measured again and remeasured after each of several pulse anneals, 
the final of these being at room temperature (23°c) when values were also 
measured during the anneal. Crystal 2, which had meanwhile been kept 
immersed in liquid nitrogen, was then similarly treated. 


§ 3. RESULTS 
The results of the measurements on crystals 1 and 2 are given in figs. 
2 and 3 respectively. 
The final values measured at room temperature after irradiation at 


—195°o (see figs. 2 and 3) are comparable with those obtained after: 


neutron irradiation of a crystal to the same dose at room temperature. 
Tables 1 and 2 show that the logarithmic decrement measured at — 195°C 
decreases appreciably only after the anneal at room temperature. 
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Logarithmic decrement at various strain amplitudes measured at room 
temperature for crystal 1. 


é Fig. 3 
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Logarithmic decrement at various strain amplitudes measured at room 
temperature for crystal 2. 


Table 1. Table of the Logarithmic Decrements Extrapolated to Zero 
Strain Amplitude Measured after the Sequence of Anneals for Crystal 1 


T Logarithmic decrement x 10° 
emperature 
of 
After After After After 
eet ee me Ave é . es: 1 hour 1 hour 1} hours | 20 hours 
c grown | torsion | irradiation St 11 4°o | at —60°o at 0°C at 23°c 
23 9-5 5:5 = = — — 1-5 
—19 2-5 2-0 2-0 2-0 2-0 2-0 1-0 
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Table 2. Table of the Logarithmic Decrements Extrapolated to Zero 
Strain Amplitude Measured after the Sequence of Anneals for Crystal 2 
ee ee eee 


Logarithmic decrement x 10* 


Temperature: |—-&@$- << —<—$<$<—<— <_< rrr |e 
of After After After After 
mecenent tee” _ After 1 hour 22 min 20 min 113 min 
26 grown | irradiation | 2+ _g9g°0 | at —60°c | at 0°0 at 23°C 
23 22:5 ae: = a = 3-5 
Sra5 3:5 3-1 34 4-2 4-2 95 


§ 3. Discussion 


These results are similar to that of Dieckamp and Sosin (1956) who 
found that a room temperature anneal could complete a Young’s Modulus 
change in copper that was begun by electron bombardment at — 180°c. 

It is generally accepted that the internal friction measured in the kilo- 
cycle per second frequency range is due to the oscillatory motion of dis- 
locations lines and that it is drastically reduced by the introduction of 
pinning agents along the dislocation lines. That the damping is not 
appreciably affected by the fast neutron dose associated with 10% thermal 
neutrons/cm? at —195°c, suggests that insufficient extra pinning points 
are introduced. As interstitial atoms are thought to be mobile at — 195°C 
they must either have little ability to pin dislocation lines which seems 
unlikely, or they very rarely reach the dislocation lines. The latter would 
occur if they precipitate as platelets in the crystal, or are trapped at 
impurities; as crowdions perhaps (Lomer and Cottrell 1955). However, 
on warming to a temperature greater than about 0°c the internal friction 
is reduced, suggesting that the dislocation lines become heavily pinned. 
This temperature range is that in which vacancies are thought to be mobile 
in copper (Broom and Ham 1957) which suggests that vacancies move to the 
dislocation lines at these temperatures and pin them. The experiment 
thus supplies evidence that vacancies are the important pinning agents 
introduced during the irradiation. 
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ABSTRACT 


Transmission electron micrographs of aged aluminium 4%-copper foils 
have shown unusual contrast near precipitated particles. This is interpreted 
in terms of Bragg diffraction from regions distorted by elastic coherency 
strains between the precipitate and the matrix. 


THE first phase precipitated when a supersaturated solid solution de- 
composes is usually coherent with the matrix, i.e. there is complete 
continuity between the atom planes of the precipitate and the matrix. 
If the normal lattice parameters of the two phases are not equal, high 
elastic strains may be generated near the interface and localized lattice 
distortion occurs in both the precipitate and the matrix (Nabarro 1940). 
The strains in the matrix are thought to account for some of the hardening 
produced by this type of precipitation process. 

When a supersaturated solid solution of aluminium 4°%-copper is aged 
at low temperatures the first formed precipitates are discs, one atom thick, 
of copper aligned on the {100} planes of the aluminium matrix. These 
discs are known as G.P. [1] zones and their formation brings about harden - 
ing of the alloy. At higher ageing temperatures or after longer ageing 
times G.P. [2] zones form and a further increase in hardness occurs. It is 
thought that the G.P. [2] zones consist of alternate copper and aluminium 
rich layers of atoms parallel to the {100} planes of the matrix. The crystal 
structure of G.P. [2] has been shown to be tetragonal with a=b=4-04 A 
whilst the ¢ parameter varies between 8-0 A and 7-6 A depending on the 
precise heat treatment. When first formed the c parameter of the G.P. [2] 
zones is 8-0 A and this changes to 7-6 A as the zones become larger (Silcock 
et al. 1953). The structure may be considered as two collapsed cubic 
aluminium unit cells. Both the G.P. [1] and G.P. [2] zones are therefore 
coherent precipitates. 

We have recently observed some unusual contrast effects in the electron 
optical images from thin foils of an aged Al 4% -Cu alloy, which we believe 
are a result of coherency strains. 

Samples of an Al 4%-Cu alloy, in the form of thin sheets, were solution 
heat treated and water-quenched to produce a supersaturated solid solution. 
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They were then aged for 16 hr at 130°c to give a structure consisting 
entirely of G.P. [1] zones and for 5 hr at 160°c to give a structure consisting 
chiefly of G.P. [2] zones. These are the structures predicted from x-ray 
diffraction studies on identical specimens by Silcock et al. (1953). Thin 
foils 200 A to 500 A thick were prepared from the sample using the technique 
described by Nicholson et al. (1958) and were examined by transmission 
in the Siemens Elmiskop I operating at 80 kv. 

The microstructure observed after the ageing treatment at 160°C is 
shown in fig. 1, Pl. 22. The G.P. [2] zones form as discs along the {100} 
planes of the matrix and as the surface of the foil is about (101) the discs 
parallel to (010) are seen as needles about 20 A thick and 300 4 long. In 
this case the precipitate diameter is about equal to the thickness of the 
foil. Under these circumstances the contrast may arise from two sources: 
the preferential scattering of the electrons due to the greater mass thickness 
of the copper rich zones, and Bragg diffraction from the crystal planes of 
the zone. For this orientation no diffracted beams are collected by the 
20 micron aperture of the objective lens and a dark image results. The 
dises parallel to (100) and (001) are not observed because their resolved 
thickness in the plane of the incident beam is insufficient to produce 
either type of contrast. If the foil is appreciably thicker than the pre- 
cipitate diameter and if the orientation of the foil surface is nearly (001) 
unusual contrast effects can be obtained by slightly tilting the specimen 
using the stereo holder. At a critical angle of tilt some of the G.P. [2] 
zones appear lighter than the matrix, whilst there are dark regions in the 
matrix adjacent to some of the precipitates. A typical micrograph is 
shown in fig. 2, Pl. 22. 

In other parts of the same grain the regions adjacent to the precipitates 
appear lighter than the matrix on one side of the precipitate and darker 
on the other side. This effect is shown in fig. 3, Pl. 28. The extreme 
sensitivity to orientation and the possible reversal of the contrast indicate 
that these effects are diffraction phenomena similar to those observed 
arcund dislocations in thin metal foils by Hirsch et al. (1956). It follows 
that the aluminium lattice adjacent to the G.P. [2] zone is displaced so 
that the crystal is no longer perfect. 

In view of the results of Sileock et al. (1953) on the crystallography of 
the G.P. [2] zones it may be assumed that the distortion of the matrix 
lattice planes adjacent to the precipitate is as shown in fig. 4, This 
requires an induced tetragonality of the face-centred cubic lattice of the 
matrix in a direction normal to the faces of the G.P. [2] zones. 

Using this model, a possible explanation of the unusual contrast effects 
can be made in terms of double diffraction. Dark field micrographs have 
shown that in certain orientations the G.P. [2] precipitates give strong 
diffracted beams. The use of a 20 micron objective aperture makes it 
impossible for any of the primary diffracted beams to contribute to the 
bright field image (Menter 1956). However, it has been shown by Pashley 
et al. (1958) that double diffraction by two crystal lattices is possible even 
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when the depth of one of them is only 50 4. If the depth of the metal foil 
is greater than the diameter of the precipitate a primary diffracted beam 
from the aluminium lattice may be diffracted a second time by the G.P. [2] 
lattice. Alternatively the order of diffraction may be reversed, as this 
depends solely on whether the precipitate is at the top or bottom of the 
foil. A simple calculation shows that doubly diffracted beams from any 
important sets of crystallographic planes will fall inside the objective 
aperture and will contribute to the bright field image of the precipitate. 


Fig. 4 


G.P[2] zone , . a4 
Q=b=4:04A crr78A . 


Aluminium matrix 
a=b=c=4.04 4 | 


TILT 


Diagram showing distortion of (100) planes near G.P. [2] zone. The dotted 
line indicates the approximate extent of the induced tetragonality (not 
to scale). 


It has been found by Pashley et al. (1957, 1958) that the doubly diffracted 
beam may interfere with the undeviated beam to give unusual contrast 
effects in the bright fieldimage. Ifthere are two infinite diffracting lattices 
a moiré pattern is produced but it is difficult to predict the precise effect 
when one of the lattices is in the form of a thin vertical disc as in the present 
experiments. In fact the precipitates appear black, white or as in a few 
cases in fig. 2, Pl. 22, a double line of white and black. Thus the contrast 
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probably depends on the exact diffracting conditions arising from the 
size and relative position of the precipitate. 

A similar argument may be advanced to explain the abnormal contrast 
in the regions of the matrix adjacent to the precipitate. In this case the 
lattice displacement extends over a comparatively large area, therefore 
the difference in lattice parameter is very small. If the two lattices were 
infinite a very large moiré spacing would result and this may explain why 
the small regions do not show the banded contrast characteristic of moiré 
patterns. As before these regions may appear dark or light but for reasons 
which are not fully understood there is sometimes no contrast near a 
precipitate or contrast on one side only. 

So far it has not been possible to obtain direct proof of double diffraction 
in these foils. The difficulties arise from the weakness of the doubly 
diffracted beam (due to the very small volume of the foil which is contri- 
buting to the beam) and its proximity to the very intense undeviated beam. 
The angular separation of the two beams is less than 10~* radians. How- 
ever, indirect proof is provided by two phenomena. No contrast effects 
are observed when there is only one diffracting lattice, i.e. when the pre- 
cipitate diameter exceeds the maximum thickness of the foil (about 500 4A), 
yet at this stage the precipitates are still fully coherent and large local 
strains would be expected. Secondly the contrast effects always disappear 
towards the edge of the specimen. An example of this is shown in fig. 5, 
Pl. 23 where the edge is just beyond the right-hand side of the picture. 
The mode of preparation of these foils necessarily involves the production 
of a thin wedge rather than a uniform thin film. This disappearance of the 
contrast effects in and near the precipitates is again associated with a reduc- 
tion in foil thickness and the consequent elimination of the two superposed 
lattices. 

If this explanation is correct the contrast effects observed near the 
precipitates arise from coherency strains induced by the tetragonality of 
the G.P. [2] structure. The nature of these strains is of considerable 
importance due to their influence on dislocation movement. With 
increasing size of the G.P. [2] zones and greater density of precipitation 
the strain to be accommodated in the matrix becomes larger and eventually 
the whole of the regions between the precipitates will be strained. This 
effect is shown in fig. 5, Pl. 23. 

The formation of G.P. [1] zones which is associated with the first increase 
in hardness during the ageing of a supersaturated solid solution should 
also give rise to small coherency strains. Figure 6, Pl. 24 shows a typical 
micrograph from the sample aged for 16 hr at 130°c. The G.P. [1] discs 
can be seen in only one orientation and appear as needles 4 to 6 A thick 
and about 80A in diameter. The very great density of precipitation 
leads to strains between the precipitates. The contrast effects produced 
by these strains can be seen in the micrograph. As in fig. 3, Pl. 23 and 
fig. 4 the thickness of the foil exceeds the precipitate diameter so that 
there is also abnormal contrast in some of the G.P. [1] precipitates. 
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Further experiments to determine the magnitude of the coherency 
strains associated with the formation of the G.P. zones are in progress 


and it is hoped to correlate the results with the age-hardening behaviour 
of this alloy. 
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IN a previous paper (Wohlfarth and Tonge 1957) some theoretical results 
were obtained for the remanent magnetization of single domain particles 
having mixed uniaxial anisotropies. This work has now been extended 
to the problem of the remanence of particles whose anisotropy is made up of 
a uniaxial and a cubic magnetocrystalline component. This situation 
arises, for example, in iron particles which have an additional shape or 
strain anisotropy. 

Let the direction cosines of the magnetization vector J, be x, 8, yand of the 
direction of uniaxial anisotropy J, m, n, both referred to the cube axes. 
Then the total anisotropy energy is 


H= K(a2B? + B*y? + y2a?) — s(al+ Bm + yn), evar t: 


where K, which may be positive or negative, is the usual cubic anisotropy 
constant, and s the uniaxial anisotropy constant which is here taken 
positive (negative values have also been considered). The remanence can 
be calculated from (1) by using the method of Gans (1932). This would 
necessitate locating the stationary directions of the magnetization vector 
as dependent on/, m,n and px, where 


pee SLT on Rah cy eo ae ee 
followed by an averaging process, to give the reduced remanence 
p= Lat oe) on Se ee ee 


as a function of J, m,n, and |u|. It has been found that the problem of 
locating these stationary directions is so laborious that its complete solution 
in the absence of a high speed computer is impracticable. Instead, only 
special cases of (1) have been considered, and the general case has been 
discussed qualitatively, leading to a fairly clear understanding of what a 
complete solution of the problem would give. A summary of the results 
is given in this note. 

When jz = 0, (1) leads to 6 or 8 minimum directions for K + or — respec- 
tively. As |, |increases these minima move towards the direction [/, m, n], 
and in general they reach this direction only when |u| oo. It has been 
found, however, that for values of | «| as low as 2 or 3, the angles which the 
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minimum directions make with [J, m, n] are quite small. For the special 
cases for which the directions [1, m, n] are given by 


(A) [lll], K>0; (A’‘) [111], K<0; 
(B) [110], Ks 0; (B’) 1110], K<0; 
(C) [100], Ks 0; (C’) [100], K <0; 


the number of minimum directions reduces to 2 at finite values of | )|, 1, 
say. These values of yu, are: 


(A) (A’) (B) (B’) (C) (C’) 
2 0-22 1 2 1 1 


The case (C’) was considered in one instance (magnetic field applied parallel 
to [100]) by Akulov (1933), who could not, of course, carry out an averaging 
process. 

The value of j, is, in general, above 4, but the deviations from 4 will be 
small when the minimum directions are close to [/, m,n], which, as already 
stated, happens for |u|aslowas2or3. For the special cases (A), (A’) ... 
(C’) just discussed, j,=4 for | u |S p,. 

Complete calculations of a curve relating j, with |u| are possible without 
excessive labour only for (C) and (C’), and the results are shown in the figure 
which also shows 1, for (A), (A’), (B), (B’). 


0-9) 
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To complete the discussion it would be necessary to average j,, || 
curves over all values of l, m,n. Although this is not at present possible, 
the above results show that the averaged j,, || curve will not deviate 
appreciably from 4 for values of |u| appreciably above 1, which is the point 
where the coercive forces arising from the cubic and uniaxial anisotropies 
separately are equal. 
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Iw one of the many papers (Simon 1951) in which the late Sir Francis Simon 
explored and clarified the Third Law of Thermodynamics he raised the 
following question : under what conditions can one prove that the vanishing 
of the specific heat as 7 approaches zero implies the Third Law? Simon 
himself pointed out that the theorem is true for the ideal quantum gas, and 
this case has since been discussed in greater detail by Landsberg (1954). 
The purposes of this note are to show that the theorem is true for another 
class of systems, those satisfying the Griineisen equation of state, and to 
discuss some of the implications of this result. 
The Griineisen equation of state can be written in the form 
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where pis the pressure, V is the volume, and U is the internal energy of the 
system. U,(V) is the internal energy at absolute zero and y(V) is 
Griineisen’s ‘ constant ’ (see Slater 1939 or Born and Huang 1954). 

From (1) it follows immediately that (dp/d7’),, is equal to (y/V)(dU/0T),, 
i.e. to yC,,/V where C,, is the heat capacity at constant volume. If we 
express (0p/07’),,in terms of the standard quantities 8, the volume coefficient 
of thermal expansion, and x, the isothermal compressibility, we obtain the 
result 


(8p/8T), = B/e=yC,|/V ee ee 


which is the usual Griineisen law. Since (0p/07’),, is equal to .(0S8/0V),, 
where Sis the entropy, it follows thatif, as 7’ approaches zero, C, approaches 
zero, then also (0S/0V),, must approach zero. In other words for a system 
obeying Griineisen’s law the vanishing of C,, as 7’ approaches zero implies 
the Third Law. 

The class of systems satisfying the Griineisen equation of state is certainly 
non-trivial. It will be satisfied by any collection of harmonic oscillators 
all of whose frequencies v depend on volume in the same way, and y is then 
—(dInv/din V) (Slater 1939). More generally Griineisen’s equation will 
be satisfied by any collection of harmonic oscillators in a crystal regardless 
of how their frequencies depend on volume, at sufficiently low temperatures 
and at sufficiently high temperatures (Borns 1923). At intermediate 
temperatures Barron’s theory (1955, 1957) indicates that y will appear to 
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vary with 7’ from its low temperature limit to its high temperature limit. 
(The experimental situation here is not clear: Bijl and Pullan (1955) 
find a pronounced variation in y for several metals, but their experimental 
method is questioned by Simmons and Balluffi (1957) who find much less 
variation.) In any case at sufficiently low temperatures the Griineisen 
equation should apply to the crystal lattice, and that is sufficient for. 
Third Law considerations. 

Tn conclusion we may point out that the ideal gas, for which the theorem 
was previously discussed, also satisfies the Griineisen relation (2) with a 
y equal to 2/3 (Simon 1951). It would be interesting to be able to 
establish Simon’s theorem for systems of more complex character than the 
ideal gas or ideal crystal lattice. 
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In their paper ‘On the Identification of Radio Sources’, Baade and 
Minkowski (1954) referring to the ‘ jet’ in the globular nebula NGC 4486, 
remark that ‘“‘ no possibility exists at this time of forming any hypothesis 
on the formation of the jet, the physical state of its material, and the 
mechanism which connects the existence of the jet with the observed radio 
emission”’. May I point out that many years ago—some years before 
the recognition of radio astronomy indeed—I put forward a hypothesis 
which led to the expectation that such ‘ jets ’ would exist in some Type E 
nebulae (Bruce 1944). The hypothesis was that such extensive atmo- 
spheres as exist in globular and elliptical nebulae, and in late type stars, 
would become the site of radial electric fields, which would increase until 


breakdown occurs. 
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Since then I have shown (Bruce 1955 a, b) that a mechanism does exist 
which will lead to the generation of these electric fields in gaseous atmo- 
spheres, namely the charging of dust particles or grains by asymmetrical 
impacts and the separation of the charged particles in a gravitational field. 
The necessary grains are known to be present in our own and other galaxies, 
and to be orientated in a particular direction in any one place as they would 
be by an electric field. 

Applied to the long-period variable stars the theory of electric fields 
and discharges accounts satisfactorily for the general nature of the out- 
bursts in these stars, the times involved, the nature of the light changes 
and of the gas movements (Bruce 1956, 1957 a). 

It is shown in an E.R.A. Report (Bruce 1957b) that the discharge 
hypothesis similarly accounts for many phenomena observed in the extra- 
galactic nebulae, in addition to the ‘jet’ or discharge channel in 
NGC 4486, and that there would appear to be sufficient reason to accept 
the occurrence of electrical discharges to resolve the nebular theoretical 
impasse, just as it has successively resolved the corresponding impasses on 
the terrestrial and stellar scales. 
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We have measured differential cross sections for the scattering of photons 
of mean energy 90 Mev in uranium through angles ranging from 1-2 to 
4-5 milliradians. The x-ray beam from the Oxford synchrotron opera- 
ting at 110 Mev peak energy was collimated to 2mm diameter at a distance 
of 4m from the machine target. The scattering angle was defined by an 
annular aperture 8m from the target. A total absorption Cerenkov 
counter and biased discriminator were used as a threshold detector for 
high energy photons. The uranium scatterer, placed immediately 
behind the 2 mm collimator, was 11-0 g/cm? thick, and electrons produced 
in the scatterer were removed by a magnetic field. A large background 
due to scattering by the air near the annular collimator was eliminated 
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by evacuating the space between this collimator and the scatterer, and 
cosmic ray background was reduced by gating the counter on during 
the x-ray pulse. The experimental results are shown in the figure. 
Inelastic scattering of photons can occur by the Compton effect and by 
secondary processes such as pair production followed by Bremsstrahlung. 
The differential Compton cross section is 7:3 barn/steradian and_ is 
substantially constant over the small angles investigated. The energy 
loss of Compton scattered photons is so small that the counter did not 
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Bethe and Rohrlich’s formula increased by 18 barn/steradian. Probable 
errors are shown for the experimental data. 


discriminate between Compton and elastic scattering. The intrinsic 
angular distributions of pair production and Bremsstrahlung are narrow 
enough to give a significant variation in the yield of secondary photons 
with angle, but, by using a thick scatterer, this variation is reduced to a 
few per cent over the range of angles of interest in this experiment due 
to broadening of the distribution by multiple scattering of the pair 
electrons. The number of secondary photons counted varies very rapidly 
with discriminator bias, but a rough theoretical estimate gave an effective 
cross section only slightly greater than the Compton cross section for 
the bias used. 
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The forward peak in the experimental distribution is therefore pre- 
sumably due to elastic scattering processes. Of these, nuclear Thomson 
scattering is negligible, probably of the order of 3 x 10~° barn/steradian. 
There is at present no complete theory of Rayleigh scattering at high 
energies. However, simple form factor calculations (Franz 1935) lead 
to a differential cross section of 6 barn/steradian at the smallest angle for 
which measurements were made. This is more than an order of magni- 
tude smaller than the observed elastic scattering, and it is extremely 
unlikely that corrections for relativistic effects and binding in the inter- 
mediate state could change the cross section to this extent. We conclude 
therefore that the main contribution to elastic scattering is Delbriick or 
potential scattering in the Coulomb field of the nucleus. 

Although it thus appears that Rayleigh scattering is small compared 
with Delbriick scattering at 90 Mev, it may still modify the total scattering 
significantly by interference. However, in the absence of a satisfactory 
theory of Rayleigh scattering, we have compared the experimental data 
with calculations of Delbriick scattering alone made by Bethe and 
Rohrlich (1952). The full curve in the figure shows the calculated 
cross sections superimposed on a uniform (inelastic) cross section of 
18 barn/steradian, chosen to give the best fit to the experimental points. 
The agreement is within experimental error. 

Previous evidence for the existence of Delbriick scattering (e.g. 
Goldzahl and Eberhard 1957) has been sought by comparing measured 
scattering cross sections for lead at energies of the order of 1 Mev with 
calculations of Rayleigh and Thomson scattering. At these energies, 
Delbriick scattering is small and the uncertainties in the theory of Ray- 
leigh scattering make unambiguous separation of the different elastic 
processes very difficult. The present results seem to give definite evidence 
for the existence of Delbriick scattering, and that the calculations of 
Bethe and Rohrlich are at least approximately correct. 
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REVIEWS OF BOOKS 


The Exploration of Space by Radio. By R. Hanspury Brown and A. ©. B. 
Lovett. (Chapman and Hall, 1957.) [Pp. 207.] 35s. 


THIS is a very good book. It is clearly written, attractively produced, and 
moderate in price. All branches of radio astronomy are described, including 
the scintillations of radio stars produced by the terrestrial ionosphere, and 
radio observations of the aurora. The emphasis is on work carried out at the 
Jodrell Bank Experimental Station. 

The value of the book would be increased by the addition of a full biblio- 
graphy. In some of the chapters full references are given to the original 
papers, but in others the author’s name only is given. There are a few mis- 
prints, and fig. 33 (a) is wrongly drawn. These are small points which could 
easily be put right in a later edition. B. H. B. 


Die Photographische Empfindlichkeitt. By J. W. Mircuetn. (Helwich : 
Darmstadt.) [Pp.35.] D.M. 8.40. 


THE silver halides belong to that class of silver salts which decompose when 
irradiated with light of the appropriate wavelength. The importance of the 
silver halides les in the fact that they are used in photographic emulsions. 
A theory for the formation of the silver latent image was first proposed by 
Mott and Gurney and in this article Dr. Mitchell discusses how the Mott- 
Gurney theory must be modified to bring it into line with the present day 
experimental evidence. It is now evident that crystal imperfections play a 
decisive role in the formation of the latent image and this is discussed in 
relation to crystal growth and optical and chemical sensitization. The present 
theory of silver formation differs from that proposed by Mott and Gurney in 
the positions assigned to the energy levels of silver sulphide, latent image 
specks and visible silver particles relative to the energy bands of the silver 
halide crystal. On this basis the latent image speck grows by first adsorbing a 
silver ion which traps a conduction electron during exposure and accepts 
electrons from reducing molecules during development. 

The new theory presented in this excellent review owes much to the work of 
Dr. Mitchell and his colleagues at Bristol. It describes not only the formation 
of the latent image but most of the physical and chemical aspects of photographic 
sensitivity. Bele Ey 


The Scientific Papers of Sir Geoffrey Ingram Taylor. Vol. I. Mechanics of 
Solids. Edited by G. K. Batcnetor. (Cambridge: University Press.) 
[Pp. x+593.] 75s. 


Tuts fine volume containing the scientific papers of Sir Geoffrey Ingram Taylor 
on the Mechanics of Solids is a fitting record of his outstanding work on this 
subject from 1917 until the present day. The Editor, Dr. G. K. Batchelor, 
and the Cambridge University Press are to be congratulated on its Ee 
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BOOK NOTICES 


Games and Decisions. By R. D. Luce and H. Ratrra. (London: Chapman 
& Hall.) [Pp. xi+509.] Price 70s. 


Proceedings of the Third Midwestern Conference on Solid Mechanics held at 
the University of Michigan, April 1 and 2, 1957. Edited by A. M. Kurrus. 
(Ann Arbor: University of Michigan Press.) [Pp. vi+250.] Price $5.50. 


Proceedings of the Fifth Midwestern Conference on Fluid Mechanics held at the 
University of Michigan, April 1 and 2, 1957. Edited by A. M. Kurrue, 
(Ann Arbor: University of Michigan Press.) [Pp. viii+388.] Price $8.00. 


Combination of Observations. By W. M. Smarr. (Cambridge University 
Press.) [Pp. xiv+253.] Price 35s. 


Nuclear Masses and their Determination. Proceedings of the Conference 
held in the Max Planck Institut fiir Chemie, Mainz, 10-12 July 1956. 
Edited by H. HintenBERGER. (London: Pergamon Press.) [Pp. ix+267.] 
Price 84s. 


Year Book of the Physical Society 1957. (London: Physical Society.) 
[Pp. 132.] Free to Members ; price to non-members 12s. 6d., postage 7d. 


Integral Equations. By F. G. Tricomr. (Pure and Applied Mathematics, 
Volume V.) (London and New York : Interscience Publishers.) [Pp. viii+ 
238.] Price $7.00. 
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Applied Mathematics, Volume VI.) (London and New York: Inter- 
science Publishers.) [Pp. x+364.] Price $9.50. 
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[Pp. xiv+263.] Price 3,900 frs. 
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Clarendon Press.) [Pp. vii+224.] Price 30s. 
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JSAPTIONS TO PLATES 


si 
Fig. 2 

The small angle scattering from cold rolled commercial copper. The round 

beam stop, which obscures angles less than 14°, was removed for 4 


Q> 
second to expose the direct beam. The beam stop is supported by a 
wire, the shadow of which can be seen. The large rectangular region 
has been blackened by parasitic radiation. The arrows indicate the 
positions of weak streaks due to reflections of the type 222, 311. 
Fig. 8 
Small angle scattering pattern from cold rolled commercial copper predicted 
entirely from the pole figures, figs. 6 and 7 (continuous lines) and from 
fig. 5 (broken lines). The directly measured pattern is shown in fig. 2. 
The concentric ‘reles represent scattering at 14° (size of beam stop) 
5°, 10°, 15° (largest angle observed). 
Fig. 4 
Large angle transmission photograph taken in the point focusing instrument 
showing the strong preferred orientation in the cold rolled commercial 


copper. The central blackening results from the direct beam, which 
also causes some spreading of the reflections. 


Fig 5. 
Bragg reflections from cold rolled commercial copper. The geometry of the 
photograph is described in § 2.2.1. 
Fig. 9 
The small angle scattering from a single grain of aluminium deformed 18° in 
tension. 


Fig. 10 


Large angle transmission photograph of aluminium deformed in tension (the 
same specimen as in fig. 9). 


Fig 
Small angle scattering from fatigued pure copper. The letters ‘ b’ indicate 
scattering of the type 200, 111. 
Fig. 12 
Large angle transmission photograph of fatigued pure copper (the same specimen 
as in fig. 11). 
Fig. 13 


Large angle transmission photograph of fatigued commercial copper. The 
specimen produced no observable small angle scattering. 


Fig. 14 
Small angle scattering from a single grain of fatigued aluminium. 
Fig. 15 
Large angle photograph of fatigued aluminium (same specimen as 1n fig. 14). 
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Carbon replicas of sodium chloride whiskers grown from solution 


P, CHARSLEY and P. E. RUSH Phil. Mag. Ser. 8, Vol. 3, PI. 21. 
Fig. 7 
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Carbon replicas of potassium chloride whiskers (figs. 7 and 8) and copper 
sulphate whiskers (figs. 9 and 10). 


Fig. 11 Fig. 12 
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Figures 11 and 12. Direct electron micrographs of sodium chloride whiskers 
after exposure to the electron beam. 
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Thin foil Al 4%-Cu aged 5 hr at 160°c. (x 160 000.) 


Fig. 2 


1000 A 
Thin foil Al4%-Cu aged 5 hr at 160°c. (x 160 000.) 
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Fig. 3 
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1000 A 


Thin foil Al4°%-Cu aged 5 hr at 160°C. (x 160 000.) 


Thin foil Al 4°%-Cu aged 5 hr at 160°C. (x 80 000.) 


R. B. NICHOLSON and J. NUTTING Phil. Mag. Ser. 8, Vol. 3, Pl. 24. 
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1000 A 


Thin foil Al 4%-Cu aged 16 hrs at 130°c to produce G.P. [1] zones. ( x 500 000 ) 


